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Soft leptogenesis is a scenario in which the cosmic baryon asymmetry is produced from a 
lepton asymmetry generated in the decays of heavy sneutrinos (the partners of the singlet 
neutrinos of the seesaw) and where the relevant sources of CP violation are the com- 
plex phases of soft supersymmetry-breaking terms. We explain the motivations for soft 
leptogenesis, and review its basic ingredients: the different CP-violating contributions, 
the crucial role played by thermal corrections, and the enhancement of the efficiency 
from lepton flavour effects. We also discuss the high temperature regime T > 10 7 GeV in 
which the cosmic baryon asymmetry originates from an initial asymmetry of an anoma- 
lous ij-charge, and soft leptogenesis reembodies in ij-gencsis. 

PACS numbers: 13.30.Fs, 14.60.St, 12.60.Jv, 14.80.Ly 
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1. The Baryon Asymmetry of the Universe 

1.1. Observations 

Up to date no traces of cosmological antimatter have been observed. The presence 
of a small amount of antiprotons and positrons in cosmic rays can be consistently 
explained by their secondary origin in energetic cosmic particles collisions or in 
highly energetic astrophysical processes, but no antinuclei, even as light as anti- 
deuterium or as tightly bounded as anti-a particles, has ever been detected. 

The absence of annihilation radiation pp — > . . . ir° — > ... 27 excludes signifi- 
cant matter-antimatter admixtures in objects up to the size of galactic clusters-^ 
~ 20 Mpc, while observational limits on anomalous contributions to the cosmic 
diffuse 7-ray background and the absence of distortions in the cosmic microwave 
background allows to conclude that little antimatter is to be found within ~ 1 Gpc, 
and that within our horizon an equal amount of matter and antimatter is em- 
pirically excluded!^ Of course, at larger super-horizon scales the vanishing of the 
average asymmetry cannot be excluded, and this would indeed be the case if the 
fundamental Lagrangian is C and CP symmetric and charge invariance is broken 
spontaneously!^ 

Quantitatively, the value of baryon asymmetry of the Universe is inferred from 
observations in two independent ways. The first way is by confronting the abun- 
dances of the light ele ments, D, 3 H e, 4 He, and 7 Li, with the predictions of Big Bang 
nucleosynthesis (BBN) p l k l " ! ' ! " ! ^ ! The crucial time for primordial nucleosynthesis is 
when the thermal bath temperature falls below T < 1 MeV. With the assumption 
of only three light neutrinos, these predictions depend on essentially a single pa- 
rameter, that is the difference between the number of baryons and anti-baryons 
normalized to the number of photons: 



n B - n s 

n = 



(i) 



where the subscript means "at present time" . By using only the abundance of 
deuterium, that is particularly sensitive to r\. Re f. 0] quotes: 

10 10 77 = 5.7 ±0.6 (95%c.L). (2) 

In this same range there is also an acceptable agreement among the various abun- 
dances, once theoretical uncertainties as well as statistical and systematic errors are 
accounted forP 

The second way is from measurements of the cosmic microwave background 
(CMB) anisotropics (for pedagogical reviews, see Refs.fTO l rTT j) . The crucial time for 
CMB is that of recombination, when the temperature dropped low enough that, 
in spite of the extremely large entropy, protons and electrons could form neutral 
hydrogen, which happened at T ^ 1 eV. CMB observations measure the relative 
baryon contribution to the energy density of the Universe multiplied by the square 
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of the (reduced) Hubble constant h = iJo/(100kmsec _1 Mpc" 1 ): 

n B h 2 = h 2 ^, (3) 

Pcrit 

that is related to r\ through 10 10 7? = 274 £Ib h 2 . The physical effect of the baryons 
at the onset of matter domination, which occurs quite close to the recombination 
epoch, is to provide extra gravity which enhances the compression into potential 
wells. The consequence is enhancement of the compressional phases which translates 
into enhancement of the odd peaks in the spectrum. Thus, a measurement of the 
odd/even peak disparity constrains the baryon energy density. A fit to the most 
recent observations (WMAP7 data only, assuming a ACDM model with a scale-free 
power spectrum for the primordial density fluctuations) gives at 68% c.lM^ 

10 2 ft B ft 2 = 2.258±°;°56- (4) 

There is a third way to express the baryon asymmetry of the Universe, that is by 
normalizing the baryon asymmetry to the entropy density s = 5*(27r 2 /45)T 3 , where 
5* is the number of degrees of freedom in the plasma, and T is the temperature: 



Y AB ^ 



(5) 

The relation with the previous definitions is given by the conversion factor so/n 7 o = 
7.04. Yab is a convenient quantity in theoretical studies of the generation of the 
baryon asymmetry from very early times, because it is conserved throughout the 
thermal evolution of the Universe. 

In terms of Fab the BBN results © and the CMB measurement (g|) (at 95%c.Z.) 
read: 

y ab N = ( 8 - 10 ± °- 85 ) x 1(rll > y ab B = (8-79 ± 0.44) x 10 -11 . (6) 

The impressive consistency between the determinations of the baryon density of 
the Universe from BBN and CMB that, besides being completely independent, also 
refer to epochs with a six orders of magnitude difference in temperature, provides 
a striking confirmation of the hot Big Bang cosmology. 



1.2. Theory 

From the theoretical point of view, the question is where the Universe baryon asym- 
metry comes from. Could it simply be the result of a fine tuned initial condition, 
one that would require just one quark in excess over 6,000,000 antiquarks and an 
exactly conserved baryon (or more appropriately B — L) number? The inflationary 
cosmological model excludes this possibility, and since we do not know any other 
way to construct a consistent cosmology without inflation, this veto is a very strong 
one. The argument goes as follows: the inflationary stage, that is the epoch in which 
the volume of the Universe undergoes exponential expansion, can only be successful 
if it lasts at least 65 Hubble times H]t] > 65. During this epoch the energy density 
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of relativistic baryons would drop exponentially as exp(— AH it). However, exponen- 
tial expansion requires that the total energy density is (approximately) constant. 
From Eq. (|6]) we see that just about seven Hubble times backward from the end of 
inflation ps would become 0(1) and dominate the (non-constant) Universe energy 
density, and this would destroy inflation. This simple argument implies that baryon 
number cannot be conserved, which opens the way to the possibility of generating 
the Universe baryon asymmetry dyna mica lly, a scenario that is known as baryoge- 
nesis. In fact, as Sakharov pointed out J^ the ingredients required for baryogenesis 
are three: 

(1) Baryon number violation: This condition is required in order to evolve from an 
initial state with Yab — to a state with Yab t^ 0. 

(2) C and CP violation: If either C or CP were conserved, then processes involving 
baryons would proceed at precisely the same rate as the C- or CP-conjugate pro- 
cesses involving antibaryons, with the overall effects that no baryon asymmetry 
is generated. 

(3) Out of equilibrium dynamics: Equilibrium distribution functions are determined 
solely by the particle energy E and chemical potential jjl 



(e-h)/t 



±1 , (7) 



and when charges (such as B) are not conserved, the corresponding chemical 
potentials vanish. On the other hand, because of the CPT theorem masses of 
particles and antiparticles are the same, and thus their equilibrium distributions 
must also be the same, which yields: 

d 3 p 

(2^3)' 

Although these ingredients are all present in the Standard Model (SM), so far all 
attempts to reproduce quantitatively the observed baryon asymmetry have failed. 



n B =n B = / — 3 , n C q- (8) 



(1) Baryon number is violated in the SM, and baryon number violating processes 
(sphalerons) are fast in the early UniverseM^J B violation is due to the triangle 
anomaly, and leads to processes that involve nine left-handed quarks (three 
from each generation) and three left-handed leptons (one from each generation). 
Sphalcron processes cannot mediate proton decay because of the selection rule 

AB = AL = ±3. (9) 

At zero tempera ture , the amplitude of the baryon number violating processes 
is proportional to^e~ S7T I 9 , which is too small to have any observable effect. 
At high temperatures, however, these transitions become unsuppressedf^ the 
first condition is then quantitatively realized, and would not impede successful 
baryogenesis. 

(2) The weak interactions of the SM violate C maximally while CP is violated by the 
Kobayashi-Maskawa complex phase of the Yukawa couplings!^ CP violation in 
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-20 



the SM can be parametrized by the Jarlskog invarianlM^ which is of order 10" 
Sinc e there a re practically no kinematic enhancement factors in the thermal 
bath ] 1 " ! l" l ^^ l it is then impossible to generate Yab ~ 10~ 10 . 
(3) Departur es from thermal equilibrium occur in the SM at the electroweak phase 
transition! I I Here, the non-equilibrium condition is provided by the interac- 
tions of particles with the bubble wall, as it sweeps through the plasma. The 
experimental lower bound on the Higgs mass implies, however, that this tran- 
sition is not strongly first order, as required for successful baryogenesis!^ 

This shows that baryogencsis requires new physics that extends the SM in at least 
two ways: It must introduce new sources of CP violation and it must either provide a 
departure from thermal equilibrium in addition to the electroweak phase transition 
(EWPT) or modify the EWPT itself. Some possible new physics mechanisms for 
baryogenesis are the foll owing: 

GUT baryogenesi J 24 ! 25 ! 26 ! 27 ! 28 ! 29 ! 30 ! 31 ! 32 ! 33 ! generates the baryon asym- 
metry in the out-of-equilibrium decays of heavy bosons in Grand Unified Theories 
(GUTs). The GUT baryogenesis scenario has difficulties with the non-observation 
of proton decay, which puts a lower bound on the mass of the decaying boson, and 
therefore on the reheat temperature after inflation. Simple inflation models do not 
give such a high reheat temperature. Furthermore, in the simplest GUTs, B + L is 
violated but B — L is not. Consequently, the B + L violating SM sphalerons, which 
are in equilibrium at T < 10 12 GeV, would destroy this asymmetry. 

Electroweak baryogenesi a ^4 | oo| j g a scenar i m which the departure from 
thermal equilibrium is provided by the EWPT. Models for electroweak baryogen- 
esis need a modification of the SM scalar potential such that the EWPT becomes 
first order, as well as new sources of CP violation. One example^ is the 2HDM 
(two Higgs doublet model), where the Higgs potential has more parameters and, 
unlike the SM potential, violates CP. Another well known example is the Minimal 
Supersymm etric SM (MSSM), where a light stop modifies the Higgs potential in the 
required wa y ' ' and where there are new, flavour-diagonal, CP-violating phases. 
Electroweak baryogenesis and, in particular, MSSM baryogenesis, might soon be 
subject to experimental tests at the C ERN Large Hadron Collider (LHC). 

Affleck-Dine mechanism! *^ ! u l The asymmetry arises in a classical scalar 
field, which later decays to particles. The field starts with a large expectation value, 
and rolls towards the origin. At the initial configuration displaced from the origin 
there can be contributions to the potential from baryon or lepton number violating 
interactions, that impart a net asymmetry to the rolling field. Generically, this 
mechanism could produce an asy mmetry in any combination of B and L. 

Spontaneous Baryogenesis! I I In this scenario, baryon number is an ap- 
proximate symmetry spontaneously broken at some large scale. A baryon asymme- 
try can develop while baryon violating interactions arc still in thermal equilibrium 
by using the effective breaking of CPT invariance caused by the Universe expansion, 
which breaks time-invariancc. Furthermore, both the ground state and fundamental 
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interactions in these theories can be CP conserving: the Universe as a whole is CP 
symmetric, but a period of exponential expansion blew domains of antimatter well 
outside our horizon. No sacred principles are violated, and although at first sight 
the mechanism could seem quite exotic, it is in fact rather natural. 

Leptogenesis. This scenario was first proposed by Fukugita and Yanagida in 
Ref. |43l and in its simplest and theoretically best motivat ed realization i s intrinsi- 
cally related to the seesaw mechanism for neutrino masses! I I I ' ' To imple- 
ment the seesaw, new Majorana SU(2)l singlet neutrinos with a large mass scale M 
are added to the SM particle spectrum. The complex Yukawa couplings of these new 
particles provide new sources of CP violation, departure from thermal equilibrium 
can occur if their lifetime is not much shorter than the age of the Universe when 
T ~ M, and their Majorana masses imply that lepton number is not conserved. 
A lepton asymmetry can then be generated dynamically, and SM sphalerons will 
partially convert it into a baryon asymmetry!^ A popular and well studied possibil- 
ity is "thermal leptogenesis" where the heavy Majorana neutrinos are produced by 
scatterings in the thermal bath starting from a vanishing initial abundance, so that 
their number density can be calculated solely in terms of the seesaw parameters 
and of the reheat temperature of the Universe. 

1.3. Prerequisites 

This review focuses on a particular realization of thermal leptogenesis, that was first 
proposed in Refs. 1501 151[ in which the lepton asymmetry is generated in the decays 
of heavy sneutrinos (the supersymmetric partners of the Majorana neutrinos of the 
seesaw) and where the relevant sources of CP violation are the complex phases of 
soft supersymmetry-breaking terms.* It is then clear that for reading this review 
some acquaintance with standard leptogenesis as well as with its supersymmetric 
version is necessary. Thermal leptogenesis has been studied in detail by many peo- 
ple, and many general papers and pedagogical reviews are available. Early studies 
that mainly focused on hierarchical singlet neutrinos include Refs. [53] [54j [55] [56] 
The importance of including the wave function renormalization of the decaying sin- 
glet neutrinos in calculating the CP asymmetry was recognized in Ref. 1571 Various 
reviews were written at this stage, and a pedagogical presentation that introduces 
the Boltzmann equations for thermal leptogenesis can be found in Ref. 1581 A partial 
set of thermal corrections to leptogenesis processes were first given in Ref. [59] while 
more complete and detailed calculations can be found in Refs. [60] 

All these studies did not include flavour effects that were first discussed in 
Refs. EU [62] but whose importance was fully recognized only later in Refs. [63] 
164] 1651 They can play an even more important role in soft leptogenesis (see Sec- 
tion [5]) than in standard leptogenesis. A pedagogical introduction to flavour effects 



"The idea of utilizing soft supersymmetry-breaking terms to realize low scale leptogenesis was 
first put forth in Ref. 1521 
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can be found in the review Ref. [66] together with all technical details. Short but 
self-contain ed resum es are also given in TASI lectures^ as well as in conference 
proceedings! ' ' ' Finally, a comprehensive study of supcrsymmctric lcptogcnc- 
sis in which the effects of non-supcrequilibration (see Section [6]) have been included 
for the first time can be found in Ref. [71] 



1.4. Reading this review 

This review is organized as follows: in Section [2] the basis of soft leptogenesis (SL) 
are reviewed and the main results are recapped. The relevant Lagrangian for SL is 
introduced in Section |2~T1 The CP asymmetries are derived in Section l2~2l bv using 
two different approaches. In Section 12.2.11 a held theoretical approach is followed, 
while in Section [2. 2.21 the same quantities are evaluated with a quantum mechanical 
approach. Beyond this section only the results of the field theoretical approach are 
used, and thus the reader can skip the details of the quantum mechanical approach, 
without affecting the understanding of the rest of the review. In SL thermal effects 
are needed to prevent a vanishing total CP asymmetry. This is a fundamental issue 
and is reviewed in detail in Section [231 

Section [3] begins with a general discussion fSection l3.1[) of how the appropriate 
effective theory to study dynamical processes in the early Universe can be for- 
mulated. Its aim is to render clear the different steps taken in studying SL with 
increasing degree of precision. The first step is discussed in Section 13.21 where the 
dynamics of SL in the so-called 'one flavour approximation' is addressed, and an 
initial set of Boltzmann equations is derived, in which flavour as well as other im- 
portant effects are left out. This Section is crucial to understand the dynamics of 
SL and to follow the qualitative discussion presented in Section 13.41 although the 
quantitative results, that are given in Scction [3.5[ can give at best a rough estimate 
of the baryon asymmetry yield of SL. 

The resonant enhancement of the CP asymmetries from self-energy contribu- 
tions is an important ingredient of SL, and for this type of contributions quantum 
corrections to the dynamical equations can be important. This issue is reviewed in 
Section [4] that, however, being a bit technical can be skipped at a first reading. 

The inclusion of lepton flavour effects in SL studies is mandatory, because SL 
always occurs in the flavoured regime. The role of lepton flavours is reviewed in 
Section [5] The flavoured CP asymmetries are introduced in Section 15.11 and two 
flavour structures representative of different soft supcrsymmctry breaking patterns 
are discussed in Section 15.21 Lepton flavour violation from soft breaking slepton 
masses is part of the phenomenology of supcrsymmctry, and if the related processes 
are sufficiently rapid all flavour effects would be efficiently damped. This issue is 
discussed in Section [5.31 and it is addressed again in relation with low energy data in 
Section [5.71 The network of flavoured Boltzmann equations, including also Higgs and 
other spectator effects, is presented in Section [5 .41 and in Scction [5.5l thc numerical 
results obtained with these equations are discussed. Finally, the impact that flavour 
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enhancements of the final baryon asymmetry can have on the soft supcrsymmetry- 
breaking parameter space is discussed in Section [5.61 

In the high temperature regime (T ^ 10 7 GeV) SL, as described in the previous 
sections, is no more the appropriate theory. Important modifications take place, 
that are related with the fact that reactions that depend on the soft gaugino masses 
and on the higgsino mixing paramete r a become irrelevant, and a new effective 
theory, that has been named i?-genesisj-^ should be considered instead. This is the 
topic of Section HI Various details of the construction of _R-gcncsis are discussed 
in Sections 16.11 and 16. 2[ and the corresponding Boltzmann equations are given in 
Section 16.31 A simplified scenario that illustrates rather clearly what is new in R- 
genesis with respect to SL is presented in Section 16.41 and numerical results are 
discussed in Section I6T51 

The prospect of (not) being able to experimentally verify the standard SL sce- 
nario is briefly discussed in Section [7] The variations of SL with their possible 
experimental signatures are reviewed in Section 17.11 

The main topics discussed in the review are resumed in the conclusions in Sec- 
tion while the more technical details are collected in two appendices. 
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2. Soft Leptogenesis: the Basic Ingredients 



The basic ingredients for generating a lepton asymmetry in SL are the CP asymme- 
tries induced in the decays of the right-handed sneutrinos (RHSN) by the complex 
phases of the soft supersymmetry(SUSY)-breaking terms. Starting from the relevant 
soft leptogenesis Lagrangian, in this section we compute the CP asymmetries fol- 
lowing first a field theoretical approach, and then a quantum mechanical approach. 
In spite of minor differences between the results obtained with the two approaches, 
it is found that in both cases to an excellent approximation the total CP asymme- 
tries for decays into scalars and into fermions vanish in the zero temperature limit. 
In fact, a general proof for the vanishing of the one-loop CP asymmetries in decays 
can be given without resorting to explicit computations, and will be presented in 
Section E3l 



2.1. Lagrangian for soft leptogenesis 

The superpotential for the supersymmetric seesaw model is: 

W = I^mssm + Y ia e ah N?ft a H b u + ^M^Nf, (10) 

where a, b = 0, 1 are the SU(2)l indices with eoi = — £ io = 1, a = e, jLt,r is the 
lepton flavour index and i = 1, 2, ... labels the generations of right-handed neutrinos 
(RHN) chiral supcrficlds defined according to usual convention in terms of their 
left-handed Wcyl spinor components (Nf contains scalar component N = v^. and 

fermion component (v^Y)- i a = (j>L Q ,d^) , H u = [H^,H^) are the left chiral 
supernelds of the lepton and up-type Higgs SU(2)l doublets respectively. Without 
loss of generality, one can work in the basis where the Majorana mass matrix M is 
diagonal. Notice that due to the Majorana mass term, one cannot consistently assign 
lepton number to N such that the superpotential (|10p remains invariant under 
global C/(l)z, Q . In other words, both L and L a are broken by the superpotential 

Starting from Eq. (|10p , the interaction Lagrangian density involving N = v^ + 
{vRi) an d Ni can be written as follows: 



- £i„t = Y la e ab [MtN*l a a H b u + H^Pl^N +Hl'°P L Nd a a + NiP L t a H b u ) + h.c.,(ll) 

where Pl,r = |(lT7s) are respectively the left and right chiral projectors. In 
Eq. dIT|) the SU(2) L doublets are I a = (£ La ,5^) T , H u = (H+ , H°) T , and H c u = 

I H+' C ,H®' C j . Notice that since H+ = H^ L is the left-handed positively charged 

Weyl higgsino, H+' c = H~ R is the right-handed negatively charged Weyl higgsino. 

The relevant soft SUSY-breaking terms involving Ni, the SU(2)l gauginos A 2 ' , 

the U(l)y gauginos Ai and the three sleptons £ a in the basis in which the charged 



September 7, 2011 0:24 
breaking 



Review: Lcptogencsis from soft supersymmetry 



10 Chee Sheng Fong, M.C. Gonzalez- Garcia and Enrico Nardi 



lepton Yukawa couplings are diagonal, are given by 

- L soft = MfjNfftj + ( AY ia e ab NiCH b u + IbM.N.N, + h. 



- (?7i 2 A^ P L A^ + miAiP^Ai + h.c. 



(12) 



where for simplicity proportionality of the bilinear and trilinear soft breaking terms 
to the corresponding SUSY invariant couplings has been assumed: Bi = BMi and 
Ai a = AYi a . In Section [5] this assumption will be dropped in favour of a more 
general flavour structure for the trilinear couplings Ai a = AZi a and, as we will see, 
this can result in important qualitative and quantitative differences. 

Even if the off-diagonal terms in the soft breaking mass matrix My are assumed 
to be negligible M^j <C Ma, the presence of the B term implies that the RHSN 
and anti-RHSN states mix in the mass matrix with mass eigenstates 



N_i 



-1= (e^' 2 N t 
V2 V 



i*,/2jy> 



V2 



(e** i/2 Ni - e'^ i/2 Nt ), 



where $,■ 



arg(B-Mj). The corresponding mass eigenvalues are 
M? ± = M? + M? i ±\BM i \. 



(13) 
(14) 

(15) 



In the following we will set, without loss of generality, $^ = 0, which is equivalent 
to assigning the phases only to A and Yi a . Including the soft terms from Eq. (|12j) . 
the Lagrangian involving the interactions of the (s)leptons and Higgs(inos) with 
the RHSN mass eigenstates N ±i , the RHN N t , and with the SU{2) L and U(1) Y 
gauginos, is given by 

Y, 



— &SL = 



V2 
+Y ia e 



-|e a6 {N +l [H C u h P L t a + (A + A/,) £ a a H b u 

} 



m- b P L £ a a + (A- Mi) £ a a H b 
N % P L £ a a H b 



ab [H u ' 



+.92 (<7±u (xjp L e a a ?: + m' a PL\tH b : 

+^= (a 3 ) ab (IPlC^* +WP L \ a 2 H b ; 

+ ^ab \^l{yZLPL-y? R PR)Z a a t +W*P L \ X HZ*j +h.C, (16) 

where g 2 and gy are respectively the SU(2)l and U(l)y gauge couplings, yn, = — 1 
and yiR = 2 denote respectively the hypcrcharges of the left- and right-handed 
(s)leptons, and a± = (01 ± ia 2 ) /2 with ui the Pauli matrices. 

All the parameters appearing in the supcrpotcntial (|10j) and in the Lagrangian 
()12|) (or cquivalently in the first three lines of Eq. (|16jl) are in principle complex 
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quantities. However, superficial phase redefinition allows to remove several phases. 
In the following, for simplicity, we concentrate on SL arising from a single RHSN 
generation i = 1 and to simplify notations we will drop that index (Y a = Yi Q , Z a = 
Z\ a , B = -Bn, etc.).* After supcrficld phase rotations, the relevant Lagrangian 
terms restricted to i = 1 arc characterized by only three independent physical 
phases: 

<t> A = aig(AB*), (17) 

^Tg(Bm* 2 ), (18) 

iarg(Sm*), (19) 

which can be assigned to A, and to the gaugino coupling operators 172, 9y respec- 
tively. Thus, in the calculation of the CP asymmetry described below M, B, 7712, mi 
and Y a correspond to real and positive parameters, while A, 52 and gy are complex 
quantities with respective phases cf>A, 4> g2 , and 4>g Y - 
The tree-level RHSN decay width is given by 



■'il2 



'ilY 






Rc(A) ( B\ UP B 2 „ /r , % 



(20) 



4tt 
where 

fc^A^A (21) 

and 5s <C 1 is assumed. Neglecting SUSY-breaking effects in the RHSN masses and 
in the vertex, we have 

r tf + * r *L* rs ^E# (22) 

a 

2.2. CP asymmetries 

The total CP asymmetry in the decays of N± is defined as: 

Y^ n(Ni -> a a ) - 7(JVi -> a Q ; 



t-±,a D 



(23) 



where 7(iV,; — > a Q ) is the thermally averaged decay rate''' for the decay of Ni into 
final state a Q (a a = s a , f a with s a = ^-ff^ and / a = £^H^' h ). 

* This simplification does not imply any crucial loss of generality. As it is explained in detail in 
Ref. 1731 the dynamics of the heavier leptogenesis states can become important only in temperature 
regimes in which the flavours of the leptons are not completely resolved by their Yukawa mediated 
interactions with the Higgs. The relevant temperature range falls in any case above T ~ 2 X 10 9 GcV 
(see Section [SJ, while SL can proceed successfully only at lower temperatures. 

'The thermally averaged reaction density is defined in Eq. (|A.2H >. 
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H„ 



H„ 



N>. 



Nj N H 



■ N< 



Ni. 



H'„ 




Nj N { 



H' 




N 3 



Fig. 1. Feynman diagrams contributing to the RHSN self-energies at one-loop. 



A^ 




H' 




I" 



Fig. 2. Feynman diagrams contributing to the RHSN decay vertex at one-loop. 



Ignoring thermal effects and taking into account only the mass splitting in the 
decay width and amplitudes, Eq. (|23|) becomes 



£ (|iri 2 -K1 2 )M 



e a (T = 0) 



i= ±,a Q 



i— ±,ap,fi 



(24) 



where A®" is the amplitude for the decay of Ni into a a . 

To fully account for finite temperature corrections several different effects must 
be considered: 

(A) Thermal corrections to (s)leptons and Higgs(inos) propagators, 

(B) Final state statistical factors, 

(C) Thermal masses of (s)leptons and Higgs(inos), 
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(D) Thermal corrections to gauge and Yukawa couplings, 

(E) Particle motion in the thermal bath. 

In the two pioneering paper d ' ' it was shown that the most relevant thermal 
effects in SL are those of type (B) that arise from final state Bosc-enhanccment and 
Fermi-blocking for RHSN decays respectively into scalars and fermions. These effects 
spoil the cancellation between the decay asymmetries into scalars and fermions, and 
are large enough to render SL viable. In Ref . [50l only effects of type (B) were taken 
into account. In Ref. [ST] effects of types (C) and (D) were also included, but it was 
found that they did not change significantly the overall picture. However, in all 
these studies, effect (E) was always ignored. Later, the authors of Ref. [SO! studied 
the full-fledged thermal effects (A)-(D), and concluded that all the effects previously 
neglected did not introduce significant changes. As regards specifically the effects of 
type (E), Refs. l59l l60l showed that in the case of SM type I leptogenesis the related 
corrections are at most ~ 20% with respect to the T = case, which suggests that 
they can be neglected also in SL. 

Including only the main thermal effects (B), (C) and (D) the total CP asymmetry 
(|23| simplifies to: 

e Q = e +Q + e_ Q + e +Q + e_ Q , (25) 

where 



c ±a 



j2 (\Ar\ 2 + \W\ 2 )c/M/ 



(26) 



i=±,ap, 

i/c|2 _ I Afa\2 \ Ja 



|^| 2 -|41 2 4°M 

J — V L (27) 

a £ (iiri 2 +i^i 2 )c-7M/ 

In these equations the finite temperature corrections from thermal phase-space, final 
state Bose-enhancement for decays into scalars and Fermi-blocking for decays into 
fermions have been factored out in the thermal coefficients c^_° , c± , so that A" a and 

A1" are the zero temperature amplitudes. Note that as long as the zero temperature 
lepton and slepton masses and small neutrino Yukawa couplings are neglected, the 
thermal coefficients are flavour independent, and if the mass splitting between N + 
and N- is also ignored, they are the same also for i = ±. In the approximation in 
which N± decay at rest the thermal coefficients are given by: 

c{ = cJ{T) = (l-x t - x s JX(l,x t ,x s J [1 - f! g ] [l - /£] , (28) 

4 = c s (T) = \{l,x Hu ,Xi) [1 + fZ] fl + ff\ , (29) 



September 7, 2011 0:24 Review: Lcptogcncsis from soft supcrsymmctry 
breaking 



14 Ghee Sheng Fong, M.C. Gonzalez- Garcia and Enrico Nardi 

where 

A(l, x, y) = V(l + a:-y)2-4a:, x a = "^f- , (30) 

and 

f e i = \ f e J = - CS11 

J ^J oME Hu .j/T]-r J h u i exp[£ gu /T] + l' l ' 

are the Bose-Einstein and Fermi-Dirac equilibrium distributions, respectively, with 

_ M _M 

^i,h u ~ ~2~' ~*~ x i,S u ~ x h u ,i» H u ,i ~ ~2~' ~*~ x H u ,e ~ x i,hJ- \y2) 

The CP asymmetry is generated at the loop level from the interference between 
the tree- level and the one- loop diagrams shown in Figs. [1] and [21 that correspond to 
different sources of CP violation: the first one arises from the self-energy corrections 
(Fig.QJ while the second arises from vertex corrections (Fig. [2]). In the following we 
describe how the decay asymmetries are computed within two different approaches: 
the first one relies on field theory, the second one on quantum mechanics. 

2.2.1. Field theoretical approach 

When T > AM± = M + - M_, the two RHSN states are not well-separated 
particles!^ In this case, the result for the asymmetry depends on how the initial 
state is prepared. * In what follows we assume that the RHSN are in a thermal bath 
with a thermalization time T" 1 shorter than the typical oscillation time AM_jT . In 
this case coherence is lost, and it is appropriate to compute the CP asymmetries 
in terms of the mass eigenstates fjl4[) . The relevant decay amplitudes can be ob- 
tained following the effective field-theoretical approach described in Refs. [75J [76l 
[77l \78\ l79j which takes into account CP violation due to mixing and decay (as well 
as their interference) of nearly degenerate states, by using resummed propagators 
for unstable mass eigenstate particles. The decay amplitude A" a of the unstable 
external state Ni into final state a a (a a = s a , f a with s a = i a a ii\ and f a = £^H^' b ) 
is described by a superposition of amplitudes with stable final states: 

At = (^ + *vr b V)) - (a? + *n- b v)) x M2 _ M r* bs > (33) 



Ml- M^+ iE*™ ' 

.— abs 

iF=(^r+*n Qabs V)) - (<°*+*n° abs v))x - l f _ abs -(34) 



M|-M4+iE T=F 



In these equations A±" are the tree-level amplitudes: 



A 3 ? =^(A*+ M)e ab , A s _° = ~i^(A* - M)e ab , (35) 

Y Y 

A+ = -j=[u{pi,)P R v{p^ c )]e ab , At = -i-^[u(pi)P R v(p S c)]e ab . (36) 

* The effects of initial conditions in SL have been studied in Ref . 1741 
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Ef bs are the absorptive parts of the Ni — > Nj self-energies (see Fig. [T]) which can 
be obtain by directly evaluati ng t he imaginary part of the Feynman integral or by 
using Cutkosky's cutting rules!^ 

"1 M 2 
_2 + 2M 2 

-r lm{A). 



-,(l)abs __ 

-,(l)abs _ 
J =F± ~ 



rM 



\A\ 2 t Re(A) 



2M 2 



M 



(37) 
(38) 



V 



a a abs 



are the absorptive parts of the vertex corrections (see Fig. [5]) : 



V 



s a abs 



v 



v£ 



s a abs 



(p 2 ) 
(p 2 ) 



Y a 3ra 2 / N2 

e a,b — ^7~ (92) 



abs 



Cab 



V2 32tt 

Y a 3m 2 
V2 32T 
y Q 3m 2 



In 



p 2 + to?. 



-Itab 



( 52 ) 2 ln- 



mo 



jr + m 2 

2i 



v^ abs ( P 2 ) 



V2 32ttp 2 
Yq, 3to 



(A*+M)( 9 *)^ln 



(39) 

(40) 

T X [u(pi)P R v(p 8e )}, (41) 



-«e G &- 



r (A*-A/)( 52 ") 2 ln- 



'V2 32V- " p2+ ^|X[^)P^fe £ )] ) (42) 

where only the contribution from SU(2)l gauginos has been included. The contribu- 

I |2 

tion from C/(l)y gaugino can be obtained by simply substituting a 2 — > ay = ? 
and 3 ->• 1 in Eqs. (f39 |) - ((42]) . 



Substituting Eqs. <(33j) and (J34J) into Eqs. flU and ([27]) and using that S^ = 



El* and E^ bs * 



TT 



J T± 



—abs 

: E^, one gets: 



|^°| 2 - |^| 2 ~ -A\ -Im [A^*A^E^ b |] 



M 2 



M 2 



(Af 2 -Af 2 ) 2 + |E^| 2 



-Im 

him 



AJ»*VJ» abs (M; 



^^a a abs* / ^ j 2 



(M|)A" a £ 



abs 



A£ 



¥ V a ^ hs (Ml)^l 



^abs 



(A4-Af|) 2 + |E^| 2 



(43) 



where the ~ sign means that terms of order <5| and higher are ignored, with 8$ 
defined in Eq. (f2Tj) . The three terms inside curly brackets in Eq. (|43|) correspond 
respectively to CP violation in N mixing from the off-diagonal one-loop self-energies 
that will be denoted below with S (='self-encrgy'),§ CP violation due to the gaugino- 
mediated one-loop vertex corrections to the N decay^J denoted by V (='vertex'), 
and CP violation in the interference of vertex and self-energies denoted by / (=' in- 
terference'). On the other hand, the amplitudes appearing in the denominators 
in Eqs. PS ]) - ([27 |) verify the tree-level relations \A a ^\ 2 + \A^\ 2 = 2\A a £\ 2 , with 



|^ |2 



Y 2 [|A| 2 + M 2 ± 2MRe(A)] and \A f ± a 



Y 2 Ml. 



§This corresponds to the effects originally considered in Refs. 1501 1511 
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Using the explicit forms in Eqs. (j55"j) - (|4"2"j) one can verify that the the three 
contributions S, V and I to the CP asymmetry from scalar and fcrmion decays 
satisfy : 



ef a = A s (T)e s ±a , 
ef a = A s (T)el a , 
4 J a = A s (T)e£ Q , 



eH = -A>(T)e| a , 
e f Z = -Af(T)el a , 
e{' a = -Af(T)ei 



±a' 



with 

e ±a — 
e ±a = 
e ±a = 

and 



~ M Sm( -' 

3P a a 2 m 2 

8 M~ 

3P a a 2 m 2 A 

4 m"m 







In- 



2BT 



ab 2 + r 2 : 



M 2 



.4 



■ sin 



(0A + ^4>g) 



In 



+ M 2 



M 

sin(</>A)cos(2(/)g) 



M 
r 2 



sin(2</> ff ) ± sin(20 9 ) 



AB 2 + r 2 ' 



A S ,/(T) 



„V 



(T) 



5 (T) + c/(T)' 



(44) 

(45) 

,(46) 

(47) 

(48) 



I |2 

In the expressions Eqs. (|4"5j) -([4"T f we have introduced a 2 = ^-, and the physical 



complex phases 4>a and 



have been explicitly written, so that all the 



parameters A and Y a etc. are understood to be real and positive. We will adopt 
this convention also in the following, unless explicitly stated in the text. The flavour 
projectors are defined as 

Y 2 



P a = 



and satisfy the conditions 



Z^ 



P n 



zyr 



1 and < P a < 1 . 



(49) 



(50) 



Summing up the contributions from the decays of N + and iV_ into scalars and 
fcrmions, one obtains the three contributions to the total CP asymmetry Eq. (|23| . 



where 



^v _ 



sin ( 

M V 

3a 2 m 2 

A~JI 

3a 2 m 2 A 

~Y~~MM 







In- 



4(T) 

4 CO 

ABT 

AB 2 + r 2 ' 



P a e s A BF (T), 
P a e v A BF {T), 
P Q e- 7 A BF (T), 



M 2 



M 



hi- 



M 2 



sin {<S>a + 1<Pg ) 
sin(0^)cos(2(/)g) 



B 

Yi 
r 2 



sin (2(j) g ) 



AB 2 + r 2 ' 



(51) 
(52) 
(53) 

(54) 
(55) 
(56) 
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Fig. 3. The thermal factors Abf (black solid curve), A s (blue dashed curve) and A-* (red dotted 
curve) as a function of z = M/T. 



and the thermal factor Abf (T) is given by 

A BF (T) = A S (T)-A'(T). (57) 

Eq. (|51l) contains the contribution to th e asym metry due to CP violation in 
RHSN mixing discussed in the original works! I I Eqs. (|5"2"j) and (|55|) give respec- 
tively the contribution to the asymmetry from CP violation in decay and in the 
interference between mixing and decay. These last two contributions have paramet- 
ric dependence similar to the ones obtained in Ref. 1811 However, as it is explicitly 

shown in Eqs. (|44|). t he scalar and fermionic CP asymmetries cancel each other at 

fool f 

zero temperature)^ because as T — > both c s (T), c' (T) — » 1. Consequently up to 

second order in the soft parameters, all contributions to the SL CP lepton asym- 
metry require thermal effects in order to be significant. More precisely, e^ (T) and 
e^(T) vanish exactly in the T = limit, in agreement with a general proof that will 
be presented in Section [2~3l As regards e~ (T), it does not vanish exactly; however, 
the surviving terms arc of order 0(<5|) and thus completely negligible. 

Fig. [3] displays the thermal factors Abf (black solid curve), A s (blue dashed 
curve) and A-^ (red dotted curve) as a function of z = M/T. For z < 0.8, the decays 
of RHSN to scalars and fcrmions are kincmatically forbidden. In the small interval 
0-8 ;$ z < 1.2 the fermionic channel becomes accessible although the scalar channel 
is still closed; this is because the thermal masses for the fermions are half than 
the ones for the scalars. For z > 1.2, the scalar channel opens up as well, however 
because of thermal effects the cancellation between A s and A-* is not very effective, 
and for relatively small values of z a sizable total asymmetry survives. For z > 10 
thermal effects arc strongly suppressed and the cancellation becomes almost exact. 

As a final remark, let us note that in this derivation thermal corrections to the 
loop diagrams responsible for the CP asymmetries have been neglected. That is, the 
imaginary part of the one-loop graphs has been obtained by directly evaluating the 
imaginary part of the Feynman integrals or by Cutkosky's cutting rules at T = 0!^ 
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2.2.2. Quantum mechanical approach 

In this section we describe the computation of the CP asymmetry using a 
quantum m echanical (QM) approach, based on an effective (non Hcrmitian) 
Hamiltonian P u l a I " I In this language an analogy can be drawn between the N- 
N^ system and the system of neutral mesons such as K°-K , for which the time 
evolution is determined, in the non-relativistic limit, by the Hamiltonian: 

-(!!H(4¥). 

with T given in Eq. (|22]l . 

In Refs. [50j [51] [81] the QM formalism was applied for weak initial states N and 
N^ . In practice, the formalism can be applied to study the evolution of initial states 
that are either weak or mass eigenstates. In order to illustrate the dependence of 
the results on the choice of initial conditions, we compute the asymmetry for both 
types of initial states. Let us define the basis: 

iVi = (gN + hitA ) n 2 = e j s fhN - gfiA . (59) 

The mass basis introduced in Eq. (fT4")) , corresponds to (g,h,6) = (-4=, -4=, — \ ). 

Pure N and W initial states correspond instead to (g, h 7 9) — (1, 0, 7r). 

Including the one-loop contribution from gaugino exchange, the decay ampli- 
tudes of N\ and N 2 into fermions f a = ^a-ff^' 6 are: 

At = ^Y a h-^(gM + hAn(g;f^If\[u(p £ )P R v(p SCu )]e ab , 

M" = [Y a g - |^ (hM + gA) {g 2 f ^/>| [u(p Sc JP L v ( Pe )]e ab , 

At = -e~ ie \Y a g - |^ {hM - gA*) (g* 2 f ^//j [u(jpt) P R v(p Sg )]e ab , 

At = e- ie {Y a h -^(hA- gM) (g 2 ) 2 ^//j [u(p^)P L v ( Pe )]e ab , (60) 

where A denotes the decay amplitudes into antifcrmions. The corresponding decay 
amplitudes into scalar s a = i^H^ are: 

At" = |r Q (gM + hA*) - ^h (g 2 f 4^/ s } e ab , 
A^ = {Y a (hM + gA) ^g (g*f ^-I s \ e ab: 



A? = e-* W \ Y a (hM -gA*) + -^g (<? 2 ) 2 £±I. \ e ab , 



3Y a 2 m 2 

— g{g2) 16V 

~^T h{9 * 2)2 jt Is 
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Re(//) = f R = 



1 / m 

2- 



ml + M 2 



Li 2 



Re(/ fl ) = s R = - 



h 



m 2 + M 2 



Li 2 



mi + M 2 



ml + M 2 
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C(2) 



-C(2) 



B (M 2 , m 2 , 0) + B (M 2 , 0, m 2 



lm(I f ) = f I = Im(I s ) = s I 



In 



ml + M 2 ' 



In terms of Ni and N 2 the eigenvectors of the Hamiltonian are: 



N L ) = (gp + hq) N, ) + e- W (hp - gq) 



N H ) = (gp - hq) 



Nx 



(hp + gq) 



N 2 
N 2 



where 



Q_ 

P 



•1 



TA 



2 A2 



T Z A 



■sm.{(p A ) 



BM " M 2 B 2 

At time t the states N\ and N 2 evolve into 



cos 2 ((j> A ) 



i T 2 A 2 
2M 2 B 2 



sin (2(f)/ 



#i, 2 (t)) = \ {[e L (t) + e H (t) ± Co (e L (t) ~ e H (t))\ 
+e* ie C lt 2(e L (t) -ejr(*))|jV- 2>1 )} , 



^i, 



where 



Go = gh I - + - , Ci = h g -, C 2 = h g 

q p 



and 



e H , L (t) = &' i(MH ' L ~ iTH,L)t - 
The total time integrated CP asymmetry is 

J2 T{Ni ->• a a ) - T(Ni -> a a ) 



C QM 



J2 r(N l ^a p ) + r(N l ^a f3 ) ' 

i=l,2,a«,/3 



where the time integrated rates r(iVj — > a a ) can be obtained from Eq. ([65 



T(Ni -> a Q ) 



1 c a ° 



4 16ttM 



KTGW 



40a 



Gj- 



Re(^^.)Gg-Im(^*^i)^ 



(62) 



(63) 
(64) 



(65) 
(66) 
(67) 

(68) 



(69) 
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and the rates T(Ni — > a a ) for antiparticles are obtained by replacing A^ a — > A\° 
In Eq. (|69|) we have introduced the time integrated projections 



<3i(2)+ — 2 ( 2 + T~ 

' \ 1 — v 1 + : 



±8 



y 

Re (Co) 



2|G> 



1 



1 



.<; 



Gi(2)- — 2 |Ci,2| 



1-2/2 

1 1 



- Im (Co 



1 - y 2 1 + x' 
x 



1 + x 



..2 



1 — y 2 1 + x 2 

y 



Gf 1(22) = 2 JRe [e^C m ] ^- 2 Im [e^C 1(2) ] ^ } 



y 2 1 + :r' z 



±2Re [e^C *C 1(2) ] (^y 
G{ 1(22) = 2 jlm [e^C 1(2) ] ^ + Re [e^C 1(2) ] ^ 



±2Im [e^C *C 1(2) ] (- 



1 



y 2 1 + x' 

written in terms of the mass and width differences^: 

M H - M L B 1 TA 2 



T 2 BAP 



sm 10^ 



y 



r H -r L a b 

= — cos ( m a) . 

2r M K ' 2M 



Using Eqs. (J69J)— (J73J) one can write the numerator in Eq. (|68|) as 

^r(iV ?; ->■ a a ) - T(N t -> a a ) = Ar a - R + AT a "' NR + Ar Q - 7 , 



(70) 
(71) 

(72) 

(73) 



(74) 



(75) 



with 



Ar 



a a ,R 



1 c a ° 



(|c 1 | 2 -|c 2 | / 

, ,C | Ji+ -3"i- 



a; 2 + y 2 



2 167rM(l-y 2 )(l + x 2 ) 1'""' " iT 2 

+2 [Rc5- 2 -Rc(e- i0 C o *Ci) -Rey 3 +Re(e ie C *C 2 )] 

-2 [Imy 2 _lm (e- ie C^Ci) -ImJ 3+ Im (e ie C *C 2 )] I, 



(76) 



^Wc use the expression of Th — ^L from Ref. 1811 Notice that with this definition Fh — Tl ^ 
^N ~ ^N wnore Tjv is defined in Eq. J20)l . 
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Ar 



a a ,NR 



2yRe(C )Ji-+Ji+ 



16ttM (1-y 2 ) , 

+y [ReJ 2 „Re (e _ie Ci) +ReJ 3 _Re (e i6 C 2 )] 

-y [lmJ 2 _Im (e- i0 d) +ImJ 3 _Im (e ie C 2 )] L 



(77) 



where 



Ar a Q ,/ 



-2Im(C* )Ji 



16ttM(1 + x 2 ) . 

- [ReJ 2 -Re (e~ ie Ci) +ReJ 3 _Re (e ie C 2 )] 

- [lm?2-Im(e" ie Ci) +Im3 r 3 _Im (e ie C 2 )] 



(78) 



J 



±L4£° 



T 



^° 



(79) 

(80) 
(81) 



1± = |AM 2 - A»- 

J 3± = A^'Af" ilf 7 *^ . 
In writing the above equations we have classified the contributions as resonant (R) 

2,2 1 

if they include an overall factor \T^i and non-resonant (NR) if no factor of 1 .^ is 
present, while the remainder has been labeled as interference (I). After substituting 
the explicit values for the amplitudes and the coefficients, and neglecting all the 
terms that cancel in both bases, the following relations are obtained: 



AT U,R 



Ar /c 



,NR 



Ar /= 



/ AT«, 


AP- fl = c s Arf , 


\ r NR 


Ar s a ,JVil = c « Ar Wfi ) 


^r 7 

^" l - L ft! 


Ar W = c ^ATi, 



with 



ATi 



-in 

1 
x- 



Y l [(9 2 ~ h2 ) 2 + ( 2 9 h ) 2 cos(26»)] Asin( 



x 2 + y 2 



x(l-y 2 )(l + x 2 Y 



Ar 



NR. 



16tt 



Y" Q a 2 In ■ 



rn 2 



M 2 M 1 - y 2 



-Asin(</u+20 9 ) 



-yM (2(2 ff /i) 2 + (g 2 - h 2 ) 2 cos(20)) sin(20 fl )] , 



Ar 7 = — F Q 2 a 2 In 



m 2 



1 



16tt 



i 2 + M 2 M 1 



rj4 



(82) 

(83) 
(84) 

(85) 
(86) 



x sin(^) 008(20) cos(20 g ). 

Eqs. (j82"|) explicitly show that the T = cancellation of the CP asymmetries occurs 
also in the QM formalism in both cases of RHSN as initial mass or weak eigenstates. 
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Given that the dependence on the thermal factor Abf{T) Eq. ([57)1 is the same as in 
the field-theoretical approach and, after normalizing to the total decay width, the 
same projectors P a Eq. (|49)) multiply the CP asymmetries, the results can again be 
recast in terms of flavour and temperature independent quantities e defined as: 



JC)QM(is) 



(T) = P a e 



(C)QM(is) 



Abf(T), 



(87) 



where the superscript (C) = R, NR, I refers to the resonant, non-resonant, and 
interference contributions, while (is) — w,m refers to the case of weak or mass 
RHSN initial states. Substituting the values for the coefficients for initial weak 
RHSN, together with the expressions for x and y in Eqs. (|74[) and expanding at 
order 6 S , one gets 



-R.QM.w 



F NR,QM,w 



F I,QM,w 



sin ( 

M V 

3a 2 m 2 

Fi 

3a 2 ?Tl2 A 







BT 



b 2 + t 2 



In. 



M 2 

,,2 



A 
M 



sin((j)A) cos(2(j)g 



In- 



M 2 



6^)008(2^) 



B 

2M 

r 2 

s 2 + r 2 



sin (2cf) g ) 



Correspondingly, for initial N± states one gets: 



-B„QM,m 



-NR,QM,m 



-I.QAI.m 



A_ 
M 



sin ( 







BT 



B 2 +T 2 '' 



3a 2 m 2 
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(91) 
(92) 
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Comparing Eqs. (|9TJ) (J93J) with Eqs. (|88|)-([90]) and Eqs. fl54])-(|56D one sees that 
the parametric dependence is very similar, although there are some differences in 
the numerical coefficients. In particular in either the weak or mass basis e R,( ^ M , 
£i,QM anc j |- ne ^-dependent (second term) in e NR <Q M coincide with e s , e 1 and the 
B-depcndent term in e v derived in the previous section, modulo the redefinition 
A —> 2A, B — » 2B and sm(<pA) —* ± sm((f>A)- There are however, some differences in 
the phase combination which appears in the B independent term in the asymmetries 
^_r,qj\/ and ^v as seen in Eqs ^ (J52j 7 ([gg]) anc | ([92]). I n other words, the choice of 

initial state only leads to minor differences. But the crucial role of thermal effects 
to avoid exact cancellations and to allow for a non- vanishing CP asymmetry is the 
same in both the QM and field-theoretical approaches, and is independent of the 
particular basis chosen for the initial RHSN states. 



2.3. The vanishing of the CP asymmetry in decays at T = 

As we have seen in the previous two sections, the original claim that the sources of 
direct CP violation from vertex corrections involving the gauginos do not require 
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Fig. 4. Soft leptogenesis diagrams for RHSN decays into scalars (la), (lb), (lc) and into fermions 
(2a), (26), (2c). 



thermal effects to produce a sizable lepton asymmetry in the plasmaN^ is incor- 
rect, and after including vertex corrections the CP asymmetries for decays into 
scalars and into fermions still cancel in the T = limit. This issue is of some inter- 
est, because if thermal corrections are necessary for SL to work, then non-thermal 
scenarios, like the ones in which RHSN are produced by infiaton decays and the 
thermal bath remains at a temperature T <C M during the following leptogenesis 
epoch, would be completely excluded. In the following, we present a simple but 
general argument proving that at T = the direct leptonic CP violation in RHSN 
decays vanishes at one loop, due to an exact cancellation between the scalar and 
fcrmion contributions. 

Let us take for simplicity $ = in Eq. (fl4")l (this amounts to assign the phases <pA 
and 4>g hi Eqs. (fTTf and (fl"8)) respectively to A and m^ ■ Since the lepton flavour a 
will not play a role in this proof, we will suppress in this section the corresponding 
label. Let us introduce for the various amplitudes the shorthand notation A. = 



N(N*) 



A(N ± -> IH$, A 

other final states. From Eq. ([T4|) we can write 



2 Aj 



2 \Ai 



A[N(N*) — > £H^) with similar expressions for the 
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II Here we only consider the contributions from SU{2)j J gauginos since for U(1)y gaugino the 
proof proceeds in exactly the same way. 
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where the complex conjugate amplitudes in the last terms of both these equations 
have been rewritten as follows: (Af )* = A £ ~ = Af and (A 



N' 



if: 



A\ 



N* 



Af by 

using CPT invariance in the second step. The direct CP asymmetry for N± decays 
into fermions is given by the difference between Eqs. (j94| and (l95]l : 
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(96) 



With the replacements t — > i and £—>£*, a completely equivalent expression holds 
also for the decays into scalars. 

The tree-level and one-loop diagrams for the various decay amplitudes into 
scalars and fermions arc given in Fig. SJ We note at this point that A~ has no 
one-loop amplitude to interfere with (see diagram (1 a)) and thus, up to one-loop, 
the full amplitude coincides with the tree-level result, and is CP conserving. A^ 1 is 
a pure one- loop amplitude (see diagram (2 c)) and therefore is also CP conserving. 
It follows that: 
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(97) 



We can thus change simultaneously the signs of Af and Af in eq. (|M|) 
without affecting the equality, and the same we can do in the analogous equation 
for the scalars. This gives: 
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Using CPT invariance 
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and unitarity 
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(102) 



we can readily see that the sum of zero temperature fermionic CP asymmetry 
Eq. (|98|) and scalar CP asymmetry Eq. (|99|) vanishes. We have thus proved that for 
N + and N_ independently, at one loop there is an exact cancellation between the 
scalar and fcrmion final state contributions, and thus at T — the direct decay CP 
asymmetries vanish. 
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3. One-flavour Approximation and Superequilibration Regime 

3.1. Effective theories in the early Universe 

In the expanding early Universe, at each temperature T is associated a characteristic 
time scale given by the Universe age tu(T) ~ H~ l (T) (H(T) being the Hubble 
parameter at T). Particle reactions must be treated in a different way depending if 
their characteristic time scale r (given by inverse of their their thermally averaged 
rates) is: 

(i) Much shorter than the age of the Universe: r <C tjj{T); 
(ii) Much larger than the age of the Universe: r S> tjj{T)\ 
(iii) Comparable with the Universe age: r ~ tjj{T). 

The first type of reactions (i) occur very frequently during one expansion time and 
their effects can be simply 'resummed' by imposing on the thermodynamic system 
the chemical equilibrium condition appropriate for each specific reaction, that is 
y^j \xi = ^2 F Hf, where /xj denotes the chemical potential of an initial state particle, 
and hf that of a final state particle. The numerical values of the parameters that 
are responsible for these reactions only determine the precise temperature T when 
chemical equilibrium is attained and the resummation of all effects into chemical 
equilibrium conditions holds but, apart from this, have no other relevance, and do 
not appear explicitly in the effective formulation of the problem. 

Reactions of the second type (ii) cannot have any effect on the system, since they 
basically do not occur. Then all physical processes are blind to the corresponding 
parameters, that can be set to zero in the effective Lagrangian. By doing this, it 
is then easy to read out if new global symmetries appear and, if no anomalies 
are involved, these symmetries correspond to exactly conserved quantities. The 
corresponding conservation laws must be respected by the equations describing the 
dynamics of the system. 

Reactions of the third type (iii) in general violate some symmetries, and thus 
spoil the corresponding conservation conditions, but are not fast enough to enforce 
chemical equilibrium conditions. Only reactions of this type appear explicitly in the 
formulation of the problem (they generally enter into a set of Boltzmann equations 
for the evolution of the system) and only the corresponding parameters represent 
fundamental quantities in the specific effective theory. 

Several examples of the importance of using the appropriate early Universe effec- 
tive theory can be found in leptogenesis studies. Leptogenesis was first formulated 
in the so-called 'one flavour approximation' in which a single SU(2)l lepton doublet 
of an unspecified flavour is assumed to couple to the lightest singlet seesaw neutrino, 
and it is thus responsible for the generation of the lepton asymmetry. Indeed, until 
the works in Refs. I63[ I64| most leptogenesis studies were carried out within this 
framework. Nowadays, it is well understood that the 'one flavour approximation' 
gives a rather rough and often unreliable description of leptogenesis dynamics in 
the regime when the flavours of the leptons are identified by in-cquilibrium charged 
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leptons Yukawa reactions. This is because such an 'approximation' has no control 
over the effects that are neglected, and thus the related uncertainty cannot be es- 
timated. On the other hand, if leptogenesis occurs above T ~ 10 12 GcV, when all 
the charged leptons Yukawa reactions have characteristic time scales much larger 
than tu, the 'one flavour approximation' is not at all an approximation. Rather, it 
is the correct high temperature effective theory that must be used to compute the 
baryon asymmetry. The corresponding effective Lagrangian is obtained by setting 
to zero, in the first place, all the charged lepton Yukawa couplings, so that the only 
remaining flavour structure is determined by the Yukawa couplings of the heavy 
Majorana neutrinos. 

In supersymmetric leptogenesis instead, the effective theory that was generally 
used was in fact only appropriate for temperatures much lower than the typical 
temperatures T ^> 10 8 GcV in which leptogenesis can be successful, and only quite 
recently it was clarified that in t he re levant temperature range a completely differ- 
ent effective theory holds insteadM-^ More specifically, it was always assumed that 
lcpton-slcpton reactions like e.g. ££ <h> ££ that are induced by soft SUSY-breaking 
gaugino masses arc in thermal equilibrium, and this implies equilibration between 
the leptons and slcptons density asymmetries (superequilihration). Supcrequilibra- 
tion (SE) instead, only occurs below T ~ 10 7 GeV, and thus supersymmetric lepto- 
genesis always proceeds in the non-superequilibration (NSE) regime. 

As regards SL, it always occurs in a temperature regime in which the charged 
lepton Yukawa couplings cannot be set to zero, and thus flavour effects must be taken 
into account, while, since SL can be successful from T ~ 10 8 GeV downwards, the 
two possibilities that it will occur in the SE or in the NS E regimes remain open. 

Here, as it was done in the original formulation ] ^" ' we first describe SL taking 
into account only reactions of type (iii). That is, we will neglect all considerations 
about flavour effects, that are related to reactions of type (i), as well as NSE effects, 
that are related to reactions of type (ii). These two issues are addressed respectively 
in Section [5] and in Section [6] 

3.2. Boltzmann equations in the unflavoured approximation 

In order to quantify the parameter ranges in which SL is successful one needs to 
solve the relevant set of Boltzmann equations (BE). All technical details about the 
BE for SL are given in |Appcndix A| 

To eliminate the dependence on the expansion of the Universe it is customary 
to recast the BE in terms of the variables Yx — n x/s, that is in terms of particle 
number densities nx normalized to the entropy density s = -jr-g^T 3 where g* is 
the total number of relativistic degrees of freedom. To account for the sources of 
violation of lepton number one then needs to follow the evolution of Y^ and, since 
RHN decays are also AL = 1 processes, the evolution of Y/v must also be considered. 

To simplify the understanding of how a sizable density asymmetry is dynamically 
generated it is convenient to adopt a certain number of approximations. 
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The first approximation is to neglect lepton flavours, and work in the so-called 
'one flavour approximation'. The relevant quantities one wants to estimate in this 
case are the fermionic Y&e and scalar Y^j lepton asymmetries generated in the 
leptonic states coupled to the RHSN (that in general correspond to a superposition 
of the different lepton flavours) . They are defined respectively as Y A i = (Y/ — Yjt) /gt 
and Y A j = \Yj—Yj^)/gj, that is, we define the density asymmetries for single 
SU(2)l degree of freedom, with gg = gj = 2. 

The second approximation is to neglect all "spectator effects" IgalSa Of course, 
besides the lepton density asymmetries, many other asymmetries related to the fi- 
nite chemical potentials of the Higgs, higgsinos, quarks and squarks, SU(2)l singlet 
leptons and sleptons, are also present in the plas ma, an d affect indirectly the out- 
come of SL through the so-called spectator effects! I I In this section all effects of 
this type will be neglected, which amounts to assume that all particles except the 
heavy (s)neutrinos and the SU(2)l doublet (s)lcptons follow either Bose-Einstein 
or Fermi-Dirac distribution with vanishing chemical potential / = (e E ' T =F 1) _1 - 

A third simplification arises from the fact that at relatively low temperatures 
(T < 10 7 GeV) reactions that transform leptons into sleptons and vice versa are 
much faster than the Universe expansion rate. Consequently, the chemical potentials 
of lepton and slcpton equilibrate m = \xj or equivalcntly -y% = tM-, a condition 

known as SE. In the NSE regime T J^ 10 7 GeV interesting new effects arise that, 
however, introduce highly non-trivial modifications in the description of SL. For this 
reason in this section SE is assumed even when the relevant temperature regimes 
fall above T ~ 10 7 GeV. 

Neglecting SUSY-breaking effects in the RHSN masses and in the vertices, all 
the amplitudes for N + and N_ decays are equal, as well as their corresponding 
equilibrium number densities, n~ = n e ~ = n e ~. Thus, in this approximation, a 
unique BE for Y^ = Y^ + Yjj suffices to account for the RHSN densities that, 
together with the BE for Yjv, give two equations for the out-of-equilibrium heavy 
neutral states. Using the SE condition 2Yg = Yj one can combine the BE for the 
unflavoured asymmetries Yai and Y A j into a single equation by defining a global 
density asymmetry in the SU(2)l lepton doublets 

Y Aitat = 2 (Y Ae + Y Al ) , (103) 

where the factor of 2 comes from summing over the SU(2)l degrees of freedom. We 
can also define a total CP asymmetry 

e(T) = e s (T) + e f (T) = I ■ A BF (T), (104) 

where 

e *J(T)=± J2 eOA^iT), e= ^ e", (105) 

q=S,VJ q=S,V,I 

and the thermal factor Abf{T) is given in Eq. f|5T[) . Given that in the one flavour 
approximation all lepton flavours are treated on an equal footing, it is left under- 
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stood that in the previous equations the various components of the CP asymmetry 
have been simply summed over lepton flavour e q = ^ a e q a . The relevant parameters 
that appear in the CP asymmetries then are A, rri2, B, M and the two CP-violating 
phases 4>a and (j> g . The BE for the unflavoured case read: 

(106) 
,(107) 

(108) 

where the time derivative is defined as Yx = s H z ^jf~ with z = M/T, Y~ q 
and Yf q = , 4 a 5 . The washout term W in the equation for Yat,, . reads: 

W = \ ( 7 tf + 7.) + jf + ^§f ( 7t (5) + lr) + ^ (2 7t (0) + 7 f 

^y \ / N 

+ 2 ( 7t (1) + 7t (2) + 7t (4) + 7t (6) + iP + 7 t (9) ) + (i + \yW-) 722 ■ (109) 

Assuming Maxwell-Boltzmann equilibrium distribution, the RHSN and RHN equi- 
librium abundances can be written as: 

n q = ^-/ % ^ Y * = ^i/ % ^- ( n °) 

(See Ref. |85]for a discussion of the validity of the use of integrated BE.) 

The derivation of the factorization of the relevant CP asymmetries including the 
thermal effects is somewhat lengthy but straightforward (see | Appendix POl ). The 




different 7's are the thermally averaged reaction densities for the different processes 
(they are defined in |Appcndix A.2.T ). In all cases a sum over the CP conjugate final 



states and lepton flavours is left implicit. 

Eqs. (|106[) - (|108[) include the N± and N decay and inverse decay processes as well 
as all the AL = 1 scattering processes induced by the top-quark Yukawa coupling. 
AL = 2 processes involving the on-shell exchange of N or N± are already accounted 
for by the decay and inverse decay processes. The AL = 2 off-shell scatterings 
involving the pole-subtracted s-channel and the u- and t-channcls, as well as the 
the .L-conserving processes from N and TV pair creation and annihilation, have not 
been included. The reaction rates for these processes are quartic in the neutrino 
Yukawa couplings and therefore can be safely neglected as long as these couplings 
are much smaller than one, as it is the case for the relevant mass range M < 
10 9 GcV required for successful SL (see next section). The non- resonant AL = 2 
processes only become important (strongly suppressing the final asymmetry) when 
the neutrino Yukawa couplings become of order of one which implies M > 10 GeV 
(see e.g. Ref. [60j). Note that in Eqs. ()106|) - ()108p only the CP asymmetry in the 
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N± two body decays has been included, while CP violating effects in three body 
decays and in scattering s] ' " I po I od I s ' I have been left out. Strictly speaking, when 
the effects of washout from scatterings are included, for consistency one should 
include also the corresponding CP asymmetries. However, in the case of standard 
leptogenesis it has been found that CP asymmetries in scatterings are important 
(and dominant) only at high temperatures z < 0.5!^ Hence they are only releva nt 
in the weak washout regime, and in the case of zero initial RHN abundance! I I In 
this case, the inclusion of the scattering CP asymmetries suppresses the final lepton 
asymmetry because it results in a balance between the two opposite sign lepton 
asymmetries respectively generated during the RHN production phase and when 
the RHN eventually decay away, giving rise to a strong cancellation which, in the 
limit of vanishing washout, is actually exact !^ In SL, however, the inclusion of the 
CP asymmetries in scatterings is not straightforward, because scattering thermal 
factors constitute a new set of non trivial quantities. Nevertheless, it is reasonable 
to expect that at least for the strong washout regime the effects of scattering CP 
asymmetries are negligible also in SL. Having said that, a careful quantitative study 
in this direction is still lacking. 

3.3. Leptogenesis efficiency 

The effectiveness of leptogenesis for producing a final lepton asymmetry Y^ = 
^A^tot ( z ~^ °°) ( or ^a° _ if -^ violation from sphalerons is accounted for) could 
be conveniently parametrized in terms of the fractional amount of the maximum 
available asymmetry that is eventually converted into Y£% . However, such a 
parametrization can be consistently introduced only for standard thermal lepto- 
genesis when it occurs at temperatures above the onset of flavour effects, and this 
is because only in this case the maximum available asymmetry can be reliably es- 
timated. In this case Y^ q = Y^(z — > 0) corresponds to the maximum possible 
density of decaying RHN, and an amount of L-asymmetry equals to e is produced 
in each decay. Then the maximum available asymmetry is e ■ Y^ q , and one can write 

Y£ tot =V-eY3 fi , (111) 

where r\ is a non-negative parameter satisfying < r\ < 1 that represents the 
leptogenesis efficiency. 

However, in other scenarios different from unflavoured thermal leptogenesis it 
is more difficult to determine the maximum am ount of available asymmetry. For 
example, if the RHN are produced non thermallyp^ it can easily happen that Y$ 
is much larger than Y^ q , and in this case, if Y^ is still expressed in units of 
eY^ 9 , values of r\ > 1 will result. Needless to say, this does not mean that the 
efficiency of the leptogenesis dynamics is higher than 100%, but it simply follows 
from underestimating the maximum amount of available asymmetry. 

In the presence of flavour effects, the available amount of CP violation is no 
more described by the total CP asymmetry summed over lepton flavours e = ^ Q e Q 
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(a = e, (j,, t) but rather by the three flavoured CP asymmetries e a , and it can easily 
occur that the absolute value of some (or even of all) flavoured CP asymmetries are 
larger than the absolute value of the total CP asymmetry, with some e a having a 
sign opposite to the one of el^l Clearly, also in this case e • Y^ 9 does not account 
for the maximum available asymmetry, and since particular flavour configurations 
can produce Y^ e with a sign opposite to the one of the total CP asymmetry e, 
using Eq. (|llll) could even result in negative values of the 'efficiency' r\. 

In SL, estimating the maximum amount of available asymmetry is basically 
an impossible task. This is because, besides the effects of lepton flavours (that 
is mandatory to include in reliable SL numerical studies, see Section [5]) the CP 
asymmetries for RHSN decays into scalars and fermions e s 'f(T) depend on the 
temperature and, as it is depicted in Fig. [3l the total CP asymmetry obtained from 
their sum Eq. (|104|) can have different signs, depending on the temperature interval 
considered. Nevertheless, it became customary, and it is often convenient, to express 
the effectiveness of SL in generating a lepton asymmetry in terms of the fractional 
amount r\ of a large reference asymmetry 2iY~ q , that is: 



yoo 
I A£ tat 



2eY e ~ q0 

N 



(112) 



with a similar definition if Y^° _ is considered instead. In the denominator of the 
r.h.s. of Eq. (|1 12[) the factor 2 has been included because there are two RHSN states, 
while Y"~ 9 = 45 /(27r 4 <7*) is defined for one degree of freedom. Solving the BE for SL 
then effectively means finding the value of r\ for the specific SL setup. The value of n 
takes into account the possible inefficiency in the production of RHSN in the weak 
washout regime, the erasure of the asymmetry by L-violating washout processes, 
and the temperature dependence of the CP asymmetry through the thermal factor 
Abf{T). In the more complete treatment of Section [5l r\ will also include the effects 
of flavours and of spectator processes, and in Section[6]of the non-supercquilibration 
of the particles and sparticles density asymmetries. 

Note that, although as discussed above in general cases r\ does not correspond to 
an efficiency, often in comparing different SL setups with equal initial Y~ and equal 
total (or flavoured) CP asymmetries, the ratios of the different 77's do correspond 
to the ratios of the corresponding efficiencies, and for this reason we will follow 
the general convention of referring to r\ as to the SL efficiency. Note also that 
the relative sign between e and Y^ can sometimes be important to understand 
the details of the SL dynamics, however, as defined in Eq. (|112p . 77 is always a 
positive quantity. Nevertheless, since the sign of e is determined by soft SUSY- 
brcaking phases whose values are presently unknown, and unlikely to be measured in 
foreseeable experiments (see Scction[7|), from the practical point of view no relevant 
information is lost in characterizing the results through the 'efficiency' r\. 
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3.4. Boltzmann equations: qualitative discussion 

Eqs. (|106|) — ()108[) constitute a rather nontransparent set of differential equations. In 
order to illustrate their physical content let us discuss an oversimplified example 
that, although it refers to A?" decays, it still captures the most relevant features of 
the general mechanism of leptogenesis. Let us write down simplified BE under the 
assumption that only the decays of N are relevant to generate the lepton asymmetry 
^A^tot = 2(Ya^ + ^a?) ( wnere the factor of 2 takes into account the two SU(2)l 
degrees of freedom) and let us describe the evolution of Y^ and Y&i tot by including 
only decays and inverse decays: 

^= d { y n- n q ) . ( ii3 ) 



eD{Yx-YZ?)-W ID Y AiM , (114) 



dz \ N N 

where e is the CP asymmetry parameter, and the decay and washout (inverse decay) 
terms are respectively given by 

V t \ Y~ q 

with % n the modified Bcssel function of the second kind of order n. From Eq. f|l 14[) 
we see that in thermal equilibrium, when Y^ = Y~ q the source term vanishes and 
no asymmetry can be generated. Let us define the decay parameter K as the ratio 
between the RHSN decay width F and the Universe expansion rate at T = M 

H(M) = V^ft£ : 

K = H(M) = ^- (U6) 

In this equation we have introduced the effective neutrino mass parametc 



= ^I»- ( 117 ) 



m e ff — 

M 

a 

with v u = v sin/3 (with u=174 GeV) the vacuum expectation value (VEV) of the 
up- type Higgs doublet, and tan/3 = v u /vd with Vd the VEV of the down- type Higgs 
doublet. Note that although m e ff is related to the light neutrino mass matrix, it has 
no direct connection with its eigenvalues, and therefore it is generally treated as a 
free parameter. The equilibrium mass appearing in the denominator of the second 

equality in Eq. (fTTBj) is defined as to* = ^v/^, where M pl = 1.22 x 10 19 GeV 
is the Planck mass. In the MSSM g* = 228.75, yielding m* = 7.8 x 10~ 4 cV. 

Clearly m s, or equivalently K, characterizes the condition for the RHSN decays 
to be in equilibrium or out of equilibrium at z = 1: the strong washout regime 
corresponds to K 3> 1, the weak washout regime to K <C 1, while the intermediate 
washout regime corresponds to K ~ 1. Another factor that concurs to determine 
the final result (in the weak washout regime) is the assumed initial abundance of 
RHSN Yft(z —> 0). Two possibilities are generally considered: 
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1. Vanishing initial abundance Y^{z — > 0) = Y~ = 0. This case relies on the 
assumption that the population of RHSN is generated only through neutrino 
Yukawa interactions in the thermal bath. 

2. Thermal initial abundance Y~ = Y~ q (z -> 0) = Y~ q . This possibility can 
be realized if the RHSN have additional interactions with the particles in the 
plasma that at early times are fast enough to generate a thermal abundance. 

Qualitatively, if K 3> 1 decays occur rapidly and quickly generate a lepton 
asymmetry. However, inverse decays are also fast and efficiently erase it. In this 
case, irrespectively of the initial abundance, Y^ approaches its thermal abundance 
already at z < 1, and any lepton asymmetry generated in the early N production 
phase, as well as any preexisting asymmetry generated through some other mech- 
anisms (e.g. from the decays of the heavier TV) gets washed out completely. The 
final lepton asymmetry can be generated only when z > 1, that is when the N 
decays start occurring out of equilibrium (i.e. YVf > Y~ q ), and leptogenesis pro- 
ceeds until the last of TV's have decayed away. In this regime r), and hence the final 
lepton asymmetry, decreases with increasing values of K because the time at which 
an asymmetry can be generated is shifted towards larger values of z, where the N 
abundance gets exponentially suppressed by the Boltzmann factor. In SL the sup- 
pression effect with increasing K is even much larger than in standard leptogenesis 
because, as discussed above, the CP asymmetry quickly decreases with decreasing 
temperatures. 

When K <C 1 the washout of the lepton asymmetry is negligible, and the initial 
conditions play an important role. Assuming a thermal initial abundance Y~ = Y~ q 
and taking (just for exemplification) a constant CP asymmetry e(T) = e the final 
lepton asymmetry saturates to the maximum possible value Y^ sa e Y~ q that is 
rj = 1. On the other hand, for zero initial abundance Y~ = and K -C 1, basically 

no TV's would decay because none would be produced in first place, and thus no 
asymmetry can be generated. Relaxing the condition to K < 1 a "wrong" sign 
lepton asymmetry is generated as long as inverse decays keep populating the N 
degree of freedom (i.e. Y^ < Y~ q ).** Since the washout is weak, this asymmetry 
only suffers mild washout effects. Eventually, when at z > 1 inverse decays start 
becoming Boltzmann suppressed and slow down, the out of equilibrium N decays 
take over (i.e. Y^ > Y~ q ) producing a "right" sign asymmetry. Because all washout 
processes are now Boltzmann suppressed, this asymmetry suffers an even milder 
erasure than the "wrong" sign one, and the imperfect cancellation between the two 
asymmetries of opposite signs results in a non-vanishing Y^ e . In this regime the 
final asymmetry increases with the value of K because of two reasons: the first one 
is that the total N population is created solely through its Yukawa interactions, and 
thus the larger is K the larger is the N abundance. The second reason is that larger 



** Notice that labeling with "right" or "wrong" sign of the asymmetry is completely arbitrary. 
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K implies stronger washouts processes, and this enhances the imbalance between 
the "wrong" and "right" sign asymmetries. 

Finally, for K ~ 1 and vanishing TV initial abundance a thermal abundance 
is still reached at z ~ 1, while all washout processes remain as small as possible. 
This is the 'optimal' regime for thermal leptogenesis, that mediates between the 
requirement of generating the largest possible N abundance, while at the same 
time minimizing washout effects. 



3.5. Quantitative results 

Reliable quantitative results for 1a4 , can only be obtained by solving numerically 
Eqs. ()106|) - (|108p . Before embarking in the details of the analysis, let us remark 
that VA^tot is n °t the most convenient quantity for writing the BE for SL. This is 
because heavy (s)neutrino decays are not the only source of lepton number viola- 
tion: sphaleron transitions, that are the crucial processes to realize baryogenesis via 
leptogenesis, also violate lepton number, and in the temperature regime in which 
SL can take place they proceed with in-equilibrium rates violating L at a fast pace. 
The quantity that is best suited for numerical studies of leptogenesis is the density 
asymmetry Ya b _ l (or in the flavoured case the asymmetries of the flavour charges 
A Q = B/3 — L a y r ' ' This is because sphalerons conserve B — L, and thus N 
and N related processes are the only ones that can generate such an asymmetry or 
change its value. However, to relate the asymmetry Y&g tat that is generated by N 
and N decays exclusively in the SU(2)l lepton doublets to Ya b _ l that is given by 
a sum over the asymmetries of all the particles with non- vanishing B — L, requires 
also a detailed knowledge of the network of B and L conserving processes that are 
in thermal equilibrium, and this is because through these processes the asymmetry 
generated in the decays gets spread among all types of particle species. We will delay 
the details of the evaluation of the Ya b _ l <->• Y^i tot conversion factors to Section [5] 
and, a s anti cipated, here we will ignore sphalerons as well as all other spectator 
effects! I[ This boils down to take simply 

dY A B - L = dY Aetot 
dz dz 

and the efficiency in producing the B — L asymmetry can be expressed in terms of 
■q Eq. (fTT2|) with the replacement Y^ tot -> Y A X> B _ L . 

After SL is over, the L and B asymmetries keep being converted from one into 
the other by the sphaleron processes until at the EWPT or slightly after it, B + L 
violation gets switched off. How much Yab is generated from a certain amount 
of Ya bl then depends on the number and types of particles that are present 
in the bath with large (thermal) abundances when sphaleron processes drop out 
of equilibrium. Assuming that at the EWPT all supersymmctric particles already 
decayed away or have negligible residual densities, and that the only remaining 
relativistic degrees of freedom are the SM states and the up-type and down-type 
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Higgs doublets, the relation between Y^ B and Y^° _ i: 

Yab = -^Y? b _ l . (119) 

This relati on can change somewhat if the EW sphaleron processes decouple after 
the EWPT ' ' or if threshold effects f or hea vy sparticles or particles like the top 
quark and Higgs are taken into account! I I 

Solving the BE (|106p - (|108[) one can obtain r\ for different choices of the relevant 
parameters m c g and M. Fig. [7] displays r\ as a function of m e g for M = 10 7 GeV 
and for the two initial conditions discussed above, and it shows how in the strong 
washout regime, the efficiency is independent of the initial conditions. This is also 
illustrated by the evolution of Y& B _ L in the strong regime for both thermal and 
zero initial RHSN abundances (bottom panels of Fig. [5] and Fig. |6]). 

Notice that with thermal initial RHSN abundance and in the weak washout 
regime, the efficiency does not flatten out to a maximum value as we would have 
expected if the CP asymmetry were constant, i.e. temperature independent. What 
we observe in Fig. [7J is that in this case, r/ (dashed red curve) decreases with de- 
creasing m e ff due to the temperature dependence of the CP asymmetry. As m c ff 
decreases and Yukawa interactions become correspondingly weaker, the RHSN de- 
cay at a later time (see the top panel of Fig. [5]) when the CP asymmetry is smaller 
(see Fig. |3]), and this explains the smaller efficiency. In the strong washout regime, 
the efficiency decreases with increasing m e ff due to increasing washout (see the 
bottom panel of Fig. [5]). If the CP asymmetry were constant, the efficiency would 
decrease roughly as ~ l/m e g (see e.g. Ref. [91]for a discussion of leptogenesis in the 
strong washout regime). However, larger m e g also shifts towards smaller tempera- 
tures the moment when the B — L asymmetry is generated. Because of the strong 
temperature dependence of the CP asymmetry in SL, this implies that the efficiency 
decreases faster ~ l/m e s as can be seen from Fig. [71 

The solid black curve in Fig. [Jj shows that with zero initial RHSN abundance, 
the efficiency ry quickly drops to zero somewhere around the intermediate washout 
regime, to rise again for larger values of m e s- This corresponds to a change of sign in 
the ratio Y^° _ /e that occurs for the following reason: during the RHSN production 
phase (i.e. Y^ < 2Y~ q ), the "wrong" sign lepton asymmetry is generated. In the 
weak washout regime, a large part of "wrong" sign asymmetry survives because the 
washouts are weak and also because the "right" sign asymmetry is generated at later 
times when the CP asymmetry is smaller. As a result, the "right" sign asymmetry 
cannot overcome the "wrong" sign one (see the top panel of Fig. [6|). In the strong 
washout regime, the washout of the initial "wrong" sign lepton asymmetry is more 
efficient and also the RHSN will decay earlier when the CP asymmetry is larger. 
The combination of these two effects results in a final "right" sign lepton asymmetry 
(see the bottom panel of Fig. [6j. In the intermediate regime, a perfect cancellation 
between the "wrong" and "right" sign lepton asymmetries occurs in the dip observed 
in Fig. [Jj where the efficiency vanishes. 
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Fig. 5. Evolution of Y~ (black solid curve) and Y/± BL (red dashed curve) assuming an initial 

thermal RHSN abundance Y° = 2Yt q0 for m cff = 10" 4 eV (top) and m cff = 10" 2 eV (bottom). 

The equilibrium RHSN abundance Y- 9 is given by the gray dotted curve. 

JVtot 

The upper panels in Fig. [8] show the regions of parameters for which CP violation 
from pure mixing effects (e s ) can produce the observed asymmetry. Due to the 
resonant nature of this contribution, these effects are sufficiently la rge onl y when 
B ~ 0(r) which, as discussed in the the original proposals of SLJ ' ' implies 
unconventionally small values of B and an upper bound M < 10 9 GeV.^ 

The central panels of Fig. [5] give the corresponding regions for which CP vio- 
lation from gaugino- induced vertex effects (e v ) can produce the observed baryon 
asymmetry. Despite being higher order in 5s and proportional to the square of the 
gauge couplings (0:2), this contribution can be relevant because it dominates for 



tt Models that can naturally yield small values of B are considered in Rcfs. 92 93, 94, 95 
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conventional values of the B parameter. However, in order to overcome the 5s and 
«2 suppression this contribution can only be sizable if the RHSN are light M < 10 6 
GeV (with the approximation used in this work: 6s <C 1, A,rri2 ~ O(TeV)). 

The figure corresponds to values of the parameters such that the second term 
in Eg. ([55)1 dominates, so that the allowed region depicts a lower bound on B. Con- 
versely, when the first term in Eq. (j52")l dominates, e v becomes independent of B. 
In this case, for a given value of M and 6s the produced baryon asymmetry can be 
sizable within the range of m G ff values for which r\ is large enough. For example for 
M = 10 5 GeV, and m% = A = 1 TeV and | sin(<^4 -I- 2(j) g ) = 1| with vanishing initial 
conditions 
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or 8 x 10~ 4 < — < 3 x 10~ 2 , (120) 
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Fig. 7. Efficiency factor r) as a function of m e fj for M = 10 GeV and tan/3 = 30, and for 
vanishing (solid black curve) and thermal initial RHSN abundance (dashed red curve). 

where each range corresponds to a sign of the CP phase sm((pA + 2</> ff ) 

Finally we show in the lower panels of Fig. [5] the values of B and m e g for 
which enough baryon asymmetry can be generated from the interference of mixing 
and vertex corrections e 1 in Eq. (|53p . Generically speaking, e 1 is subdominant with 
respect to e s , since both involve the same CP phase sin (^4) but e 1 has additional 
5s and a 2 suppression: 



3 m 2 

—a? In 

8 M 



— — — — 77 cos(26 a ) — 



(121) 



Consequently, as it is also shown by the figure, e 1 can dominate only if B is extremely 
small (B <C T), when it becomes independent of B. Note also that for M < 10 4 GeV 
and ra e ff > 10 -2 eV the baryon asymmetry generated by this contribution becomes 
independent of m c g. This is because in this regime of strong washouts the rn^ s 
dependence from T cancels the approximate 1/m 2 ^ dependence of r\. 



4. The Possible Role of Quantum Effects 

Most studies of thermal leptogenesis (both for the standard seesaw case, as well 
as for the SL scenario) rely on the classical BE approach that was described in 
the previous section. The possibility of using quantum BE (QBE) in leptogene- 
sis was first discussed in Ref. I96[ and more recently analyzed in greater detail in 
Refs. [UlinHlIMl EM EEQ EH EH In Ref. W\ QBE were obtained starting from 
the non-equilibrium quantum field theory based on the Closed Time-Path (CTP) 
formulation, and differ from the classical BE in that they contain integrals over the 
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Fig. 8. B, m c g regions in which successful SL can be achieved (flavour and spectator effects 
neglected). In all cases we take A = rri2 = 10 3 GeV and tan/3 = 30 and different values of M 
and <j>j\ and (f> g as labeled in the figure (see text for details). The left (right) panels correspond to 
vanishing (thermal) initial N abundance . 



past times. In the classical kinetic theory instead the scattering terms do not include 
any integral over the past history of the system, which is equivalent to assuming 
that any collision in the plasma is independent of the previous ones. In the CTP 
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formalism, the energy conservation delta functions appearing in the evaluation of 
the reaction rates are substituted by retarded time integrals of time-dependent ker- 
nels, and cosine functions whose arguments are the energy involved in the reactions. 
In the limit in which the time range of the kernels is shorter than the relaxation 
time of the particle abundances, and the time integrals are taken over an infinite 
time (i.e. neglecting memory effects), the standard time- independent reaction rates 
are recovered. Furthermore, the CP asymmetry also acquires an additional time- 
dependent piece, that at any given instant depends upon the previous history of 
the system. 

In Rcf.[97lit was argued that the additional time dependence of the CP asymme- 
try is quantitatively the most relevant effect for leptogenesis. However, if the time 
variation of the CP asymmetry is shorter than the relaxation time of the particles 
abundances, the solutions to the quantum and the classical Boltzmann equations 
are expected to differ only by terms of the order of the ratio of the time-scale of the 
CP asymmetry to the relaxation time-scale of the distribution. This is typically the 
case in thermal leptogenesis with hierarchical RUN. Conversely in the resonant lep- 
togenesis scenario, (Mj — Mj) is of the order of the decay rate of the RH neutrinos. 
As a consequence the typical time-scale to build up coherently the time-dependent 
CP asymmetry, which is of the order of (Mj — Mi) -1 , can be larger than the time- 
scale for the change of the abundance of the RHN. As shown in Refs. 11041 I105[ in 
the case of resonant leptogenesis this leads to quantitative differences between the 
classical and quantum approach and, in particular, in the weak washout regime this 
can enhance the produced asymmetry. 

Since in SL the CP asymmetry in mixing (Eq. (1511) ) is pro duced resonantly, we 
can expect that this type of effects could be of some relevance! I 

4.1. Modification to the CP asymmetry and quantification 



We have seen in Section [2721 that the relevant CP asymmetry in SL is temperature 
(i.e. time) dependent already in the classical approximation. The inclusion of quan- 
tum effects introduces an additional time dependence. As shown in Refs. [97l 11041 
11051 quantum effects are flavour independent as long as the damping rates of the 
leptons are taken to be flavour independent. Neglecting also the difference in the 
width of the two RHSN, one can show that 

e s (T) =e s x A BF (T) x QC(t), (122) 

where I s is defined in Eq. (|54| and 

QC(t) = 2 sin 2 ( M+ ~ M ~ t) - M+ T _ M _ sin ((M + - M.)t) . (123) 

The factor QC(t) is the one that remains after taking the past time integral to 
large time such that only on-shcll decay processes contribute to the CP asymmetry 
(which is equivalent to neglecting memory effects in decay processes). This factor 
grows for t < 1/AM and starts oscillating for t > 1/AM. The oscillation pattern 
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originates from the CP-violating decays of two mixed states N+ and 7V_ analogous 
to the CP violation in neutral meson systems. If the timescale for the decay t ~ 1/r 
is much larger than 1/AM, the CP asymmetry should average to the classical value. 
However, if the decay timescale t ~ 1/r is shorter than 1/AM, this additional time 
dependence on CP asymmetry may not be negligible. 

As usual, it is convenient to change in Eq. (|123|) from time t to z = M/T. 
For a Universe undergoing adiabatic expansion the entropy per comoving volume 
is constant, i.e. sR 3 = constant, and since s oc z~ 3 then R oc z. Thus the Hubble 
parameter is given by H = R~ 1 dR/dt = z~ 1 dz/dt. After integration, one gets 

t= } * 2 ~ Z ° , (124) 

H(M) 2 ' v ' 

where zq is the temperature at t = 0. Substituting Eq. (|124[) into Eq. (|123l) : 

2 ( 1 M+ - M_ 2 \ T (M+- M_ 2 



QC(z) = 2 sin (j mM) - z ) M+ _ M _ sin y ^ (M) 

= 2sin 2 f^i 1 >^-|smf^i ? ^, (125) 

V to* 8 / R \ m * 4 / y ' 

where zq has been set equal to (corresponding to a very high initial temperature) 
and M + — M_ = B has been used (assuming M -C M (see Eq. ((T5|) V Finally, the 
degeneracy parameter R has been defined as: 

2(M + -M_) 2B , N 

i?=^^ ^ = — . (126) 

In summary, the final CP asymmetry consists of three factors: the first one is 
the temperature independent piece e which can be rewritten as^T 

q ImA 2R 

I s = . (127 

M R+l v ' 

which, as discussed before, it is resonantly enhanced for R = 1. The second one is 
the thermal factor Abf(T) which is non-vanishing only for z > 0.8 (see Fig. [3]). 
The third one is the quantum correction factor QC{z). 

The impact of this additional quantum time-dependence of the CP asymmetry 
on the final baryon asymmetry can be easily quantified by introducing QC(z) in the 
relevant BE ((T06|) . (fi"07| and (|T08| . Fig. U shows the evolution of the asymmetry 
with and without the inclusion of the quantum correction factor for several values 
of the washout parameter m e s both for the resonant case R = 1 and for the very 
degenerate case R — 2 x 10~ 4 . The two upper panels correspond to strong and 
moderate washout regimes, while the lower two correspond to weak and very weak 
washout regimes. 

tTQuantum effects for CP asymmetries in decay or interference of decay and mixing can be 
introduced in similar fashion. However, in the interesting parameter space for e v where R ^> 1, 
these effects are irrelevant, while the effects on CP violation in the interference between decay and 
mixing are expected to be of a similar size. 
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Fig. 9. Absolute value of the lepton asymmetry with the quantum time dependence of the CP 
asymmetry (solid) and without it (dashed) as a function of z for different values of m e ff as labeled 
in the figure. In each panel the two upper curves (black) correspond to the resonant case R = 1 
while the lower two curves (red) correspond to the very degenerate case R = 2 X 10~ 4 . The left 
(right) panels correspond to vanishing (thermal) initial RHSN abundance. The figure corresponds 
to M = 10 7 GeV and tan 8 = 30. 
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The figure illustrates that, as expected, for strong washouts and large degeneracy 
parameter R (see the upper curves in the upper panels), the quantum effects lead 
to an oscillation of the produced asymmetry until it averages out to the classical 
value. Conversely, for very small values of R and strong washouts, quantum effects 
enhance the final asymmetry. For small enough R the arguments in the periodic 
functions in QC(z) are very small for all relevant values of z and m G ff, so that the 
sin term is negligible. By expanding the sin term we get 

W-)-- — y, (128) 

TO* 2 

which, in the strong washout regime, is always larger than 1. 

In the weak washout regime, independently of the initial conditions and of the 
value of the degeneracy parameter R, the quantum effects always lead to a sup- 
pression of the final asymmetry. This is different from what happens in type I 
seesaw resonant leptogenesis in which for weak washouts, R ~ 1, and vanishing 
RHN initial abun dances, quantum effects lead to an enhancement of the asymme- 
try produced! I The origin of the difference is in the additional time dependence 
of the CP asymmetry in SL Abf- In order to understand how this works, we must 
remember that in the weak washout regime for type I seesaw resonant lcptogcnc- 
sis, the resulting asymmetry is the one that survives the cancellation between the 
opposite sign asymmetry generated when RH neutrinos are initially produced, and 
the asymmetry generated when they decay. Including the time-dependent quan- 
tum corr ectio ns spoils this cancellation and as a consequence a larger asymmetry is 
obtainedES 

However, in SL already in the classical approximation the thermal factor Abf 
prevents the opposite sign asymmetries cancellation, and the inclusion of the time 
dependent quantum effects only amounts to an additional multiplicative factor 
which, in this regime, is smaller than one. Consequently, for SL, even in the resonant- 
regime quantum effects do no lead to major quantitative differences, at least in the 
range of parameters for which successful leptogenesis is possible. This is explicitly 
shown in Fig. [10] that depicts the ranges of the parameters B and m c s for which 
enough asymmetry is generated (Y^ B > 8.35 x 10 -11 ), with and without the inclu- 
sion of the quantum corrections. We see that the main effect of including quantum 
corrections is that for a given value of M the allowed regions extend up to larger 
values of m c ff . This is precisely due to the suppression of the asymmetry in the weak 
washout regime just discussed. Because of the enhancement in the very degenerate, 
strong washout regime, for a given value of M the regions also tend to extend to 
lower values of B and larger values of m e g . 

A qualitative difference obtained from the inclusion of quantum effects is that, 
depending on the value of TO c g, r\ can take both signs independently of the initial 
RHSN abundance. Thus it is possible to generate the right sign asymmetry with ei- 
ther sign of Irm4 for both thermal and zero initial RHSN abundance. However, apart 
from this peculiarity, altogether it can be concluded that for a given M the values 
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Fig. 10. -B, m e ff regions in which successful SL can be achieved with (left panels) and without 
(right panels) quantum effects. We take |ImA| = 10 3 GeV and tan /? = 30 and different values 
of M as labeled in the figure. The upper (lower) panels correspond to vanishing (thermal) initial 
RHSN abundance. 



of the L-violating soft bilinear term B required to achieve successful leptogenesis 
are not substantially modified. 



September 7, 2011 0:24 Review: Lcptogcncsis from soft supcrsymmctry 
breaking 



44 Ghee Sheng Fong, M.C. Gonzalez- Garcia and Enrico Nardi 

5. The Role of Lepton Flavours 

The role of lepton flavour in the standard thermal leptogenesis scenario was 

first discussed in Ref. IFT1 However, until the importance of flavour effects 

was fully clari fie d ' _ ' they had been included in leptogenesis studies only 

in a few cases ' ' ' Nowadays flavour effects hav e been investigated in full 
detai i65 | 108 | 109 | 110 | lll | 112 | m | ll4 | 115 | 97 | 104 | 105 | 71 | and are a mandatory ingrc _ 

dicnt of any reliable analysis of lept ogenesi s. 

As regards SL, the original work a I I neglected flavour effects and considered 
the generation of the lepton asymmetry directly in the £ and £ lepton states cou- 
pled to the lightest RHSN, a situation commonly referred to as the 'one flavour 
approximation'. However, RHSN couple in fact to certain lepton combinations, in 
which the different flavours are weighted by the respective RHN Yukawa couplings: 
£ oc ^2 a Y a £ a and £ ex J2 a Y a £u- Only at very high temperatures £ and £ remain 
coherent superpositions, and are the correct states to describe the dynamics of lep- 
togenesis. At lower temperatures scatterings induced by the charged lepton Yukawa 
couplings occur at a sufficiently fast pace to distinguish the different lepton flavour 
components. In this situation, £ and £ cannot be considered anymore as coherent 
flavour superpositions, and the dynamics of leptogenesis must be described instead 
in terms of the flavour eigenstates £ a and £ a . The specific temperature when lepto- 
genesis becomes sensitive to lepton flavour dynamics can be estimated by requiring 
that the rates of processes r a (a — e, /x, r) that are induced by the charged lepton 
Yukawa couplings h a bec ome fast er than the Universe expansion rate H{T). An 
approximate relation give n I I 

r Q (T)~10- 2 ^T, (129) 

which implies that 

T a (T) > H{T) when T < T a (1 + tan j3 2 ) , (130) 

where T e ~ 4 x 10 4 GcV, T^2x 10 9 GeV, and T r -5x 10 11 GcV. Notice that to 
fully distinguish the three flavours it is sufficient that the r and fi Yukawa reactions 
attain thermal equilibrium. Therefore, given that SL can successfully proceed only 
for temperatures below T M , we can conclude that independently of the value of tan /3 
the three flavour regime is always the appropriate one to study SL. The interesting 
issue is if the region of parameter space in which SL can be successful is enlarged 
or reduced when flavour effects are taken into account, and in particular if the 
requirement of unnaturally small values of the soft B term can get relaxed. These 
issues arc addressed in the following sections. 

5.1. Flavored CP asymmetries 

In the flavour regime the CP asymmetries for RHSN decays into the single lepton 
flavours become important. Their flavour structure is determined by the Yukawa 
couplings Y a as well as by the trilinear couplings A a . The soft breaking Lagrangian 
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of the previous section Eq. (|12p assumed for simplicity a universal trilinear soft 
breaking scale A and proportionality of the soft breaking trilinear couplings to the 
Yukawa couplings, that is A a = AY a . In this section a more general flavour structure 
is considered where the ^4-terms have the generic flavour structure 

- L^l = AZ ta e ab N t I a a H b u + h.c. . (131) 

Considering only the lightest RHSN N — Ni and adopting as usual a simplified 
notation Z a = Z\ a etc., with the generalized flavour configuration in Eq. (J131I) we 
have three relevant physical phases: 

cf> Aa = aig{AZ a Y*B*). (132) 

and the CP asymmetries (J5TJ) . (15^1) and (f53")l are now written as: 



where 



e s a ^(T) = P a ^4^A BF (T), (133) 

% = -^*^mri* (134) 

_ v _ 3a 2 m 2 



eX = r- -tt In 



4 M m 2 2 + M 2 



A Y B 

— sin (4> a + 2<i„) sin(2<6„ 



(135) 



ei, = In — s — f — rsm Ui )cos(2(B ) s ^ . 136 

2 MM ml + M 2 KVa ' KV9l AB 2 + T 2 V ; 

In these equations the physical complex phases <f>A a and </> g have been explicitly 
written, so that all the parameters A, Z ai Y a etc. are real and positive. Unless 
explicitly stated in the text, this convention is adopted also in what follows. As 
regards the flavoured reaction rates, they can be simply written in terms of the 
unflavoured rates by means of the flavour projectors Eq. (|49l) : 



5.2. Flavour structures 



Ix = P a jx ■ (137) 



The flavour structure Z a of the A terms Eq. (|131j) is in principle independent from 
the flavour structure of the Yukawa couplings Y a . However, the study of flavour 
effects in a completely general flavour configuration would be rather awkward, be- 
cause of the very large dimensionality of the parameter space. It is thus convenient 
to define less general possibilities, but chosen in such a way to render possible a 
qualitative extrapolation of the impact of flavour effects to the general case. In 
Rcfs. H18l and fTT9l the following two scenari os w ere respectively introduced: 

1. Universal Trilinear Scenario (UTS)! I This case assumes universal soft 
terms, that is the soft breaking Lagrangian of Eq. (|12j) . It is realized in supcrgravity 
or gauge mediation (when the renormalization group running of the parameters is 
neglected) and corresponds to set 

Z a = Y a . (138) 
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Thus, in UTS the only flavour structure arises from the Yukawa couplings and 
both the flavoured total CP asymmetries e Q = ef + e^ + e J a and the flavoured 
washout terms W a arc proportional to the same flavour projectors P a . It follows 
that e Q oc W a , and moreover, as seen in Eq. (|132[) . the trilincar couplings have an 
unique phase (f>A a — 4>A — arg(AB*). 

2. Simplified Misaligned Scenario (SMS)! I To understand the possible effects 
of flavour dynamics, it is important to study also a case in which the flavoured CP 
asymmetries e a and the washouts W a are misaligned. This can be done without 
increasing th e nu mber of independent parameters with respect to UTS by imposing 
the conditio n I 

Z a = \ Y , , (139) 

where we have kept Z and Y explicitly as complex numbers. With this condition 
the CP asymmetries (except for the last term in Eq. (|135jl ) are equal for all flavours 
e a = e/3, while the washouts W a maintain their flavour dependence, so that an 
arbitrary misalignment can be realized. From Eq. (|132j) we see that there is again a 
unique phase (f>A a = 4>a — arg(AB*). Note that both Eq. (|138[) and Eq. (|139|) yield 
the same total asymmetry ^ Q e a = e, so that any difference between the UTS and 
SMS results can be ascribed directly to the differences in flavour configuration. Note 
also that in the case of flavour equipartition P e = P^ = P T = 1/3 the two scenarios 
are equivalent. 

Finally, it should be remarked that in both scenarios it is not possible to have 
flavour asymmetries of opposite signs, with |e Q | > |e| for some, or even for all 
flavours, as it is also excluded the possibility of having non-zero flavour asymmetries 
and a vanishing total CP asymmetry. This however, is simply due to the reduction in 
the number of physic al phases and, similarly to what happens in standard flavoured 
leptogenesisj ' ' ' in the most general case asymmetries of opposite signs (and 
possibly with a vanishing sum) are an open possibility. Given that these types of 
configurations are always characterized by larger CP asymmetries, the reader should 
keep in mind that enhancements of the final lepton asymmetry even larger than the 
ones found in UTS and SMS are certainly possible. 

5.3. Lepton flavour equilibration 

The possibility of having large enhancements of the baryon asymmetry yield of 
lcptogcncsis from flavour effects, relies on the fact that the density asymmetries 
stored in each flavour are independent from each other, and if for example a flavour 
that is weakly coupled to the washouts has a particularly large CP asymmetry, the 
final result will be essentially determined by the dynamics of this flavour. 

However, in the presence of lepton flavour violating (LFV) interactions, the 
density asymmetries in the different flavours are no more independent, and if LFV 
rates are sufficiently fast, they can bring the different flavours into equilibrium, with 
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the result that the amount of surviving asymmetry will be essentially determ ined 
by the flavour that is more strongly washed out: a potentially destructive effect! I 

In SL, LFV interactions are a natural possibility since they are directly related 
to off-diagonal entries in the soft mass matrices of the slcptons. For this reason a 
reliable study of flavoured SL must also include an analysis of LFV effects. 

In the basis where charged lcpton Yukawa couplings are diagonal, the soft slcpton 
masses read 

L soft D -mlptjp. (140) 

The off-diagonal soft slcpton masses fr? a ^a = ™^« (ct ^ 0) affect the flavour 
composition of the slcpton mass cigenstates so generically we can write 

1^ = Rapip, (141) 

where R a p is a unitary rotation matrix. In this basis the corresponding slepton- 
gaugino interactions in Eq. (|16|) become 

32 
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^=5 ab \iY iL P L t a R* a0 t£ + h.c. 



(142) 



The mixing matrix can be expressed in terms of the off-diagonal slepton masses as: 



Ra/3 ~ S a f 



hlT 2 



= 5 a f 



l a^0 



v 2 cos 2 /3 



m 2 M 2 



(143) 



where in the first line h a > hp is the relevant charged Yukawa coupling that deter- 
mines at leading order the thermal mass splittings of the sleptons, v in the second 



line is the EW symmetry breaking VEV with v 



vl + v 2 



174 GeV, 



_ M 



# and 



m a = m^ a (T = 0) is the zero temperature mass for the lepton £ a . For simplicity we 
parametrize the R a=j Lp entries in a way that they are independent of the particu- 
lar pair of leptons involved. Let us define fh eT = fh^ T = rh d and m eM = m d— -, 
where fh oe [ is a unique off-diagonal soft-mass parameter. We then obtain for the 
off-diagonal entries (a/3) = (er), (^j,t), (e/x) of the matrix R: 

m 2 od v 2 cos 2 /? 2 



R, 



a^/3 



where m T is the mass of the tau lepton. 



mlM 2 



(144) 
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Fig. 11. LFV lepton-slepton scatterings mediated by the SU(2)l and U(1)y gauginos Xq- 
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The terms in Eq. (|142[) induce LFV lepton-slepton scatterings through the ex- 
change of SU(2)l and U(l)y gauginosM^" There are two possible t-channel scat- 
terings £ a P <H- £pP*, £ a P -H- £pP and one s-channel scattering £ a £% <->• PP* (we 
denote P as fermions and P as scalars) as shown in Fig. QTJ For processes medi- 
ated by SU(2)r, gauginos P = £,Q, H u ^, while when mediated by U{l)y gaugino 
one must include the SU(2)l singlet states P = e,u, d as well. The corresponding 
reduced cross sections read: 



~a/3 
T tl,G 



-o(3 
J t2,G 



T s.G 



w 



(s) 



00 



- 


^G 


-Rq/3 


2 n£ 


p 


/ 


8tt 




r 


*Sg 


-Rq^J 


2 n£ 


p 


• 


8tt 




r 


*5g 


-R^l 


2 iEJ 



'2m 2 - 
Ag 



1 In 



m~ + s 

Ag 



Ag 



hi 



+ s 



P 



16tt 



s — m~ 

Ag . 



(145) 



where lip counts the numbers of degrees of freedom of the particle P (isospin. quark 
flavours and color) involved in the scatterings mediated by the SU{2)t, (G = 2) and 
U(l)y (G = Y) gauginos respectively. In this last case the hypercharges yu, and 
yp are also included in lip . If the flavour changing scatterin gs in Eq. f|145[) are fast 
enough, they will lead to lepton flavour equilibration (LFE) ! ' 

The values of fh d for which this occurs can be estimated by comparing the 
LFE scattering rates and the AL = 1 washout rates. Since the dominant AL = 1 
contribution arises from inverse decays, the terms to be compared are: 



■p , T > _ 7lfe(^) 
1lfe(J J = — TFFT- 



1id(J ) = 



(T)' 



(146) 



where n c L = T 3 /2 is the relevant density factor appearing in the washouts (see next 
section for more details) and 
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(148) 



In the first equality in Eq. (|147[l 7"^. with (x = tl,t2,s) is the thermally averaged 
LFE reactions for one degree of freedom of the P-particle, %\^{z) are the modified 
Bessel function of the second kind of order 1 and 2, T is the zero temperature width 
Eq. (|22j) . and n e ~(T) is the equilibrium number density for N at T. 

The LFE reaction densities, the Universe expansion rate H , and inverse decay 
rates all have a different T dependence. In particular Llfe ~ T~ 3 and H(T) ~ 
T 2 , so that once LFE reactions have attained thermal equilibrium (Flfe ~ H) 
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they will remain in thermal equilibrium also at lower temperatures. In contrast, 
for inverse decays the rate first increases, and then decreases exponentially Tid ~ 
e -M/T dropping ut of equilibrium at temperatures not much below T ~ M . The 
relevant temperature where we should compare the rates of these interactions is 
precisely when the condition Tid ~ H is reached, that is when the lepton asymmetry 
starts being generated from the out-of-equilibrium N± decays. Defining Zdec through 
the condition Tid(zdec) = H(zdec), LFE is expected to be quite relevant for flavoured 
leptogenesis when 

TlFE (Zdec) > TlD {Zdec) , (149) 

since when this condition is satisfied, LFE processes are already in thermal equilib- 
rium at the onset of the out-of-equilibrium decay era. 
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Fig. 12. The left panel shows the ratio Tlfe/^ at Zd ec as a function of m d for m e fj = 0.1 eV, 
tan/3 = 30 and three values of M. The right panel shows in the (P a m e ff, M) plane contours of 
constant values of ih ot i (in GeV) for which Tlfe { z dec) — ^"^id ( z de )• 



The left panel of Fig. [T2l shows the ratio V-LY^izdec) / H (zdec) as a function of fh d 
for ?7i ff = 0.1 eV, tan/3 = 30, and for different values of M. From the figure we can 
read the characteristic value of fh d for which LFE becomes relevant. Notice that 
the dominant dependence on tan/3 ~ l/cos/3 (tan/3 3> 1) arises from wcos/3 = Vd 
in Eq. (|144|) . Thus the results from other values of tan /3 can be easily read from 
the figure by rescaling m d(tan/3) = m o |/(30cos/3). 

Since we are interested in the dynamics of lepton flavours, to be more precise 
about LFE effects we should in fact consider the temperature at which the inverse 
decay rate for one specific flavour a goes out of equilibrium. We then define z c d y ec 
through the condition rj D (z^ ec ) = PcXiD{z dec ) = ^( z ^ ec ). For a hierarchical con- 
figuration of flavour projectors P 7 < Pp < P a we will have zj ec < z^ ec < z^ ec , 
that is the inverse decays involving the lepton doublet £ 7 which is the most weakly 
coupled to N± will go out of equilibrium first, and then r ID and r ID will follow. 
Hence, for given values of m c g and M, the minimum value fh" l d m for which LFE 
effects start being important is given by the following condition: 



Tlfe (z, 



deci ' 



T a 



ID 



\ z de, 



(150) 



September 7, 2011 0:24 Review: Lcptogcncsis from soft supcrsymmctry 
breaking 



50 Ghee Sheng Fong, M.C. Gonzalez- Garcia and Enrico Nardi 

For fh d <C in™™ LFE effects can be neglected, since they will attain thermal 
equilibrium only after leptogenesis is over. 

The right panel in Fig.[i"2"lshows in the plane of the flavoured effective decay mass 
P a m e g and of the RHN mass M, various contours corresponding to different values 
of fh d for which Tlfe ( z 2ec) = -^ib ( z dec)- For a gi yen value of M and m e s, and for 
a given set of flavour projections P 1 < Pp < P a , fn™ n is given by the value of the 
fhod curve for which the vertical line x = AI intersects the corresponding contour 
at y a = P a wieff ■ Since Tlfe has a rather strong dependence on fh d (Flfe oc fhg d ), 
one expects that the value fh^ d for which LFE effects completely equilibrate the 
asymmetries in the different lepton flavours will not be much larger than rh™ d n . 
Indeed our numerical results (see Section [57Tf show that fh eq d ~ 0(5 — 10) x rh™ d m . 
Clearly, larger values fh oc i 3> fh e od do not modify the final results with respect to 
what is obtained when m d = fh eq d . 

5.4. Flavoured Boltzmann equations 

In Rcfs.|6TJ[63l[64]the re l evan t equations including Savour effects associated to the 
charged lepton Yukawa couplings were derived in the density operator approach. 
One can define a density matrix for the difference of lepton and antilcptons such 
that p aa = Y&e a . As discussed in Refs. EU [63j EH as long as we are in the regime 
in which a given set of the chargcd-lepton Yukawa interactions are out of equilib- 
rium, one can restrict the general equation for the matrix density p to a subset of 
equations for the flavour diagonal directions p aa . In the transition regimes in which 
a given Yukawa interaction is approach ing equilib rium the off-diagonal entries of 
the density matrix cannot be neglected! I I I However, as we have discussed, 
SL always occurs in the three flavour regime, and thus there is no need to worry 
about this type of effects. In this Section we keep working, as in Section [3J under 
the assumption of superequilibration. However, we now include the effects of all the 
relevant fast reactions of type (i) (see Section 13- Xp with characteristic time scales 
much shorter than H~ l (z = 1), that is we include both flavour and spectator effects. 
This implies that we need to consider three flavoured lepton density asymmetries 
that also include the contributions of the SU(2)l singlet (s)leptons, defined as: 



Y ALa = 2 (Y A£q + Y Ala ) + Y Aea + Y A e Q , (151) 

where the factor of 2 comes from summing over the SU(2)l degrees of freedom. The 
contribution of e a and e a to the total flavour asymmetries is non- vanishing because 
scatterings with the Higgs induced by the charged lepton Yukawa couplings transfer 
part of the asymmetry generated in the SU(2)l lepton doublets to the singlets. 

In the BE for the evolution of the RHN and RHSN densities, a sum over flavour 
can be readily taken, and the resulting equations are not modified with respect to 
Eqs. (|106[) and (|107|) . The lepton charges AL Q arc anomalous, and thus are not 
only perturbatively violated in RHSN and RHN interactions, but also in nonper- 
turbativc sphalcrons transitions. This type of effects is however removed by writing 
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evolution equations for the anomaly free flavoured charges A Q = B/3 — Ln, w ith 
den sity asym metries defined as Y Aa = Vab/3 — Y A i, a . Including spectato r I I and 
j^p]jj l20 | 119 | e fr ec ^ S; th e corresponding BE is: 
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Notice 



ah u 



with F/ 9 
Y A n and we chose to express the 



that SE implies 2Y A £ 
asymmetries in term of the fermionic ones. The higgsino asymmetries enter in two 
ways: directly, when the Higgs(inos) are involved in the relevant process as external 
particles, as in RHSN decays and inverse decays; indirectly, when the top and stop 
quarks are involved in the relevant scatterings, and their chemical potentials are 
rewritten in terms of Y A ^ (see |Appendix A. 4 1 . The last line in Eq. (|152[) includes 
the reaction densities for the LFE processes in Eq. (J145I) . that as discussed above 
play the role of controlling the effectiveness of flavour effects. 

To bring Eq. (|152[) in the form of a closed differential equations for the charge 
density asymmetries Y Aa , the three asymmetries Y A i a as well as the higgsino asym- 
metry Y A jj must be expressed in terms of Y Aa according to: 



Y Ai 



Y ~ 

AH U 
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(153) 



^2-A. a pY Afj , 



The matrix A was first introduced in Ref.[6T]and the vector C in Ref. EU The values 
of their entries are determined by which interactions are in chemical equilibrium, 
and thus eventually by the specific range of temperature when leptogenesis takes 
place (see | Appendix B[ ). 



5.5. Results 

Unless differently stated, our results are obtained for M = 10 6 GeV and tan/3 = 30, 
although as long as LFE effects are negligible they are practically independent of 
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M. The dependence on tan/3 is also rather mild since it mainly arises from m e g as 
given in Eq. (J117I) through v\ ~ v 2 (l + 1/ tan/3 2 ) -1 . For tan/3 = 30 the <i-quark and 
electron Yukawa couplings arc sufficiently large to ensure that around T ~ 10 6 GeV 
the corresponding interactions arc, like all other Y ukaw a interactions, in thermal 
equilibrium. In this regime the A and C matrices are-^ 

1 /- 221 16 16 \ 4 

A= 16 -221 16 , C = (1 1 1). (154) 

9X237 V 16 16 -22l) 23? 

We anticipate that the impact of flavour and spectator corrections on the results 
will be sizable only in the strong washout regime. This can be easily understood by 
adding the equations for the three flavour asymmetries, Eq. (|152p . This results in 
an equation for Y& B _ L that can be recast in the form 

-Y AB - L =e(z)(^-2\^-WY,P a A a ^-W H Y,P a C a ^, (155) 

V N J " a/3 l a I 

where e(z) = J2 a e a{z), W is the washout term related to the lepton density asym- 
metries given in Eq. (|109[) . and Wh is an additional washout term related to the 
Higgs density asymmetry, whose expression can be easily worked out from Eq. (11521) . 
Of course, Eq. (|155p cannot be solved alone, since at least other two equations for 
two density asymmetries Ya q arc needed to get a closed system. However, since 
flavour and spectator effects arc encoded in the ^2 a g an( ^ 2 Q terms, this equation 
clearly shows that if W , Wh — > none of these effects will be important. 

The dependence of the efficiency factor on the flavour projectors P a and on 
m e g is shown in Fig. [13] The top panel shows the dependence of the efficiency on 
m e g in the flavour cquipartition case and for two other sets of flavour projections. 
As already mentioned, flavour effects become more relevant when the washouts 
get stronger. This is confirmed in this picture where it is seen that for the SMS 
scenario the possible enhancements quickly grow with m c g. Note that in SL this 
dependence is even stronger than in standard leptogenesis. This is due to the fact 
that the flavoured washout parameters P a m e g also determine the value of z°[ ec when 
the lepton asymmetry in the a flavour starts being generated, and since the CP 
asymmetry has a strong dependence on z, different values of P e , P M , and P T imply 
that the corresponding flavour asymmetries are generated with different values of 
the CP asymmetry even when, as in the SMS, the fundamental quantity e is flavour 
independent. In summary, what happens is that the flavour that suffers the weakest 
washout is also the one for which inverse-decays go out of equilibrium earlier, and 
thus also the one for which the lepton asymmetry starts being generated when 
e x Abf is larger. This realizes a very efficient scheme in which the flavour that is 
more weakly washed out has effectively the largest CP asymmetry, and this explains 
qualitatively the origin of the large enhancements that we have found. Furthermore, 
when P a m e g <C m* so that the inverse decay of flavour a never reaches equilibrium 
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Fig. 13. Top panel: the dependence of r\ on the washout parameter m e g for P e = P M and different 
values of P T . Bottom panel: the dependence of r\ for different flavour configurations P a normalized 
to rji / 3 that corresponds to the flavour equipartition case P a = 1/3. 



and the washout of the asymmetry Y& a is negligible, the maximum efficiency is 
reached. 

The bottom panel in Fig. [13] further illustrates how the departure from the 
equipartition flavour case results in an enhancement of the efficiency, and that par- 
ticularly large enhancements are possible for the SMS scenario. Note that the top 
line in the top panel of Fig. [relabeled P T = 0.99 represents the maximum enhance- 
ment that can be obtained in the SMS (relaxing the constraint in Eq. (|139p that 
defines our SMS, larger enhancements are however possible). This is because for 
P T = 0.99 both the asymmetries Y& e and Ya are generated in the weak washout 
regime, that is, approximately within the same temperature range, and in the SMS 
this implies e e (T e ) rj e At (r M ). The related combined efficiency is then simply deter- 
mined by (P e + P M ) m e ff ~ m* and is thus always maximal, independently of the 
individual values of P e and P„, as it is apparent from the figure. 
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We should however spend a word of caution for the reader about interpreting 
these results in the weak washout regime and, for the SMS, also in the limit of 
extreme flavour hierarchies (P a — > 0). At high temperatures (z < 1) the Higgs 
bosons (higgsinos) develop a sufficiently large thermal mass to decay into sleptons 
(leptons) and sneutrinos. The new CP asymmetries associated with these decays 
could be particularly largep^ and thus sizable lepton flavour asymmetries could be 
generated at high temperatures. This type of thermal effects are not included in the 
analysis here described. 

Concerning the flavour decoupling limit within the SMS, clearly when P a — > 
no asymmetry can be generated in the flavour a. However, in our SMS flavour 
asymmetries are defined to be independent of the projectors P and thus survive in 
the P — > limit. On physical grounds, one would expect for example that when one 
decay branching ratio is suppressed, say, as P < 10 -5 , the associated CP asymmetry 
will be at most of 0(10 -7 ) and thus irrelevant for lcptogcncsis. This means that for 
extreme flavour hierarchies, the SMS breaks down as a possible physical realization 
of SL, and thus in what follows we will restrict our considerations to a range of 
hierarchies P > 10~ 3 . 

As a result one finds that for the SMS scenario with hierarchical Yukawa cou- 
plings, successful leptogenesis is possible even for m e g ^> 0(eV). For example, as it 
is shown in the right panel of Fig. [13l for P e = Pfi = 5 x 10~ 3 and m e g ~ 5 cV, we 
obtain r\ ~ 10~ 3 , that yields the estimate 

Y£ B (SMS,P e = P M = 5x 10" 3 ,?7i off = 5eV)~ 10~ 6 x I. (156) 

Thus assuming a large, but still acceptable value of e ~ 10 -4 , SL can successfully 
generate the observed baryon asymmetry for values of m e g that are about two orders 
of magnitude larger than what is found in the unflavoured standard leptogenesis 
scenario. 

5.6. Natural values of B 

We next describe the impact that flavour enhancements can have in relaxing the 
constraints on the values of B (and M). Here we include only the leading CP 
asymmetry contribution from mixing Eq. (|134[) . From Eqs. (|119p . (|112[) . and (|54p . 
it is easy to derive the maximum value of B for given values of M and m e g: 



r (meg, M) |Inu4j y^ ??(m cff ) 
2 M Y™ B 




( M Y^ 



AB 



Vl Im ^l yNvim-eff) 



(157) 



where y$ = ^Y~ q0 , T(m cSl M) is given in Eq. <j22j) , ImA = Asm(j} A and Y^ IB is 
the observed baryon asymmetry Eq.(|6]). Consequently we obtain 

M < |M| ^" Kff) , (158) 

Y B obs 

3V3m off / |Im,4| y^ r]{m cS ) \ 2 
B ~^2^[ YCp ) ' (159) 
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where r](m e ff) = r](m e g, P a7 Z a ) and all residual dependence of -q on M has been 
neglected. As seen in the upper panel of Fig. [T3] assuming the SMS and for suf- 
ficiently hierarchical P a , r)(m c f{) decreases first very mildly with m e s and - once 
all the flavours have reached the strong washout regime- it decreases roughly as 
~ m~ s . Thus the product m e g x ?](m B g) 2 first grows with m c g till it reaches a 
maximum and then for sufficiently large m e g it decreases ~ m~g. Therefore, for a 
fixed value of the projectors, the upper bound on B does not corresponds simply 
to the maximum allowed value of m e g, but it has a more complicated dependence. 
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Fig. 14. Maximum values of B and M as a function of the flavour projections. To highlight the 
effects when P T is very large or very small, the curves are given as a function of P T and of 1 — P T . 



Fig. [14] shows the maximum values of B and M obtained for both the UTS 
and SMS cases as a function of the flavour projections. In order to have better 
resolution when either P T or 1 — P T is very small, we plot them both as a function 
of P T or 1 — P T . In the figure we set ImA = 1 TcV. The figure illustrates that within 
the UTS, the parameter space for successful leptogenesis is very little modified 
by departing from the flavour equipartition case (that corresponds to the point 
where the UTS and SMS curves join). On the contrary, in the SMS case with 
hierarchical flavour projections 1 — P T ~ few x 10~ 3 successful SL is allowed also 
with B ~ O(TcV), that is for quite natural values of the bilinear term. As mentioned 
above, even for hierarchical projections the maximum allowed values of B and M do 
not correspond to the maximum allowed value of m e g. In particular, for the range 
of flavour projections shown in the figure we obtain that the maximum values of B 
and M correspond to m e e ^ 2 eV. 



5.7. Lepton flavour equilibration and low energy constraints 

We finish this section by discussing the impact that the presence of the LFE scat- 
terings discussed in Section 15.31 can have on the enhancement of the efficiency due 
to flavour effects. We plot in Fig. [15] the dependence of the flavour enhancement of 
the efficiency as a function of the off-diagonal slcpton mass parameter m d- As can 
be seen in the figure (and as it was expected from the discussion in the previous sec- 
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tion) for any given value of M , LFE quickly becomes efficient damping completely 
the lepton flavours enhancements of the efficiency within a very narrow range of 



values ml 



< m od < m™ d ax . 



The figure corresponds to tan/3 = 30 however, as 



already said, the dependence on tan /3 arises from Vd — v cos /3 in Eq. (|144[) and is 
rather mild. Results from other values of tan/3 can be easily read from the figure 
by rescaling fh od = m o |/(30cos/3). 
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Fig. 15. The dependence of the efficiency (normalized to the flavour equipartition case P a = 1/3) 
on the off-diagonal soft slepton mass parameter m 0[ j, for different values of M and of the flavour 
projections (see text for details). 



Notice that while neglecting LFE the efficiency for a given value of m c g is 
practically insensitive to the particular value of M, this is not the case when the 
efficiency is evaluated by accounting for LFE effects. In fact, given the different 
scaling with the temperature of the Tlfe and Tid rates, the precise temperature 
at which lcptogcncsis occurs is crucial. For example, we see from Fig. [TS] that for 
reasonable values fh d ^S 200 GeV and for M > 10 6 GcV, LFE is not effective, and 
the large enhancements of the efficiency due to flavour effects can survive, while for 
M < 10 5 GeV all flavour enhancements disappear. 

It is interesting notice that the presence of a sizable m d would induce various 
LFV decays, like for example l a — ¥ lp~/ with rate 

>- 3 tan 2 /3 _ 4 



BR(l a -> l 0l ) 
BR{l a ->■ lpv a Vp) 



a 



-<2 



"SUSY 



, Q sin 2 /3 / 
2.9 x 10 -19 ^ - 



TeV 



cos 2 /3 



"od 



GeV 2 



(160) 



cos 6 /3 \msusY , 

where msuSY is a generic SUSY scale for the gauginos and sleptons masses running 
in the LFV loop. So it is possible to compare the values of fh d for which LFE occurs 
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(M SUSY /TeV)(cos(3) 3/4 /(sin(3) 1/4 



Fig. 16. Shaded in yellow: the region of fh 0( i cos /3 versus msj/sy- (cos /3) & /(sin/3) 4 excluded by 
the present bound BR(fi — > ej) < 1.2 X 10 -11 . The three bands corresponding to M = 10 5 GeV, 
M = 10 6 GeV and M = 10 7 GeV show the minimum value of fh od cos for which LFE effects 
start damping flavour effects. The width of the bands corresponds to P Q m e ff varying in the range 
0.003 eV— 10 eV, with P a the largest flavour projection. The vertical dashed line corresponds to 
the value of mg(/sy/(tan/3) 2 required to exp lain t he discrepancy between the SM prediction and 
the measured value of a M , assuming tan/3 = ll 



with the existing bounds imposed from non-observation of such flavour violation in 
leptonic decays. The result of such comparison in presented in Fig. 1161 The yellow 
shade gives the excluded region of fh d cos/3 versus msusY (cos /3)* /(sin (3)* arising 
from the present bound BR(/i — > ej) < 1.2 x 10 -11 , together with the minimum 
value of fh d cos /3 for which LFE effects start damping out flavour enhancements in 
SL. Three bands are shown respectively for M = 10 5 GeV, M = 10 6 GeV and M = 
10 7 GeV. The width of the bands represents the range associated with variations of 
the effective flavoured decay parameter P a m e s in the range 0.003 eV— 10 eV, where 
P a is the largest of the three flavour projections. For illustration we also show in the 
figure the characteristic SUSY scale that allows to explain the small discrepancy 
between the SM prediction and the measured value of th e mu on anomalous magnetic 
moment, a^. This values is nisusY K^&n-fi) 3 = 141 GeVJ I and the vertical dashed 
line in the picture corresponds to tan/3 = 1. As seen in the figure, in this case 
the off-diagonal slepton masses are bound to be small enough to allow for flavour 
enhancements in SL for M as low as 10 6 GeV. For larger values of tan (3, even lower 
values of M are allowed. 

In brief, LFE effects induced by off-diagonal soft slepton masses, when con- 
strained with the bounds imposed from the non-observation of flavour violation in 
leptonic decays, are ineffective for damping the flavour enhancements. 
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6. Soft Leptogenesis without Superequilibration: .R-genesis 

As mentioned in the previous Sections, early works on leptogenesis were carried out 
from the start within the unfiavourcd effective theory. Quite likely this happened 
because the corresponding Lagrangian is much simpler, given that the number of 
relevant parameters is reduced to a few. The main virtue of subsequent studies on 
lepton flavour effects was that of recognizing that for T < 10 12 GeV the unfiavoured 
theory breaks down, and the new theory brings in new fundamental parameters, 
which can give genuinely different answers for the amount of baryon asymmetry 
that is generated. 

In supersymmetric leptogenesis the opposite happened: the effective theory that 
was generally used assumed fast particle-sparticle equilibration reaction. But this 
assumption is only appropriate for temperatures much lower than the typical tem- 
peratures T 3> 10 8 GeV in which leptogenesis can be successful. In fact, only quite 
recently Ref. [71] clarified that in the relevant temperature range, a completely differ- 
ent effective theory holds instead. More specifically, in supersymmetric leptogenesis 
studies prior to Ref. 1711 it was always assumed (often implicitly) that lcpton-slcpton 
reactions like e.g. ££ <H> ££ (see Fig. [T7|) that are induced by soft gaugino masses 

L~ x = -\ (m 2 A^°P L A±'° + miAiPtAi + h.c.) , (161) 

as well as higgsino mixing transitions H u -h- Hd, that are induced by the superpo- 
tcntial term 

W H = fiH u H d , (162) 

are in thermal equilibrium. This implies that the lepton and slcpton chemical poten- 
tials equilibrate. However, in general, in supersymmetric leptogenesis superequili- 
bration (SE) between particles and sparticles chemical potentials does not occur. In 
fact, the rates of interactions induced by SUSY-breaking scale (A susy ) parameters, 
like soft gaugino masses rrig or the higgsino mixing parameter /^, are slower than 
the Universe expansion rate when 

%<2 5 £ * r>T SE ~5.10'(^) 2/3 G„V. ,163, 

Thus, when this condition is realized, these reactions should be classified as reactions 
of type (ii) (see Section I3.1|) and handled accordingly. Since leptogenesis occurs 
when the temperature is of the order of the heavy neutrino mass, in terms of M 
the assumption of SE breaks down when 



^susy \ 

500 



2/3 



M>5.10 7 (^g^J GeV. (164) 



Following the discussion in Section [3. 1[ the effective theory appropriate for studying 
supersymmetric leptogenesis, in which the heavy Majorana masses certainly satisfy 
the bound Eq. (|164|) . is thus obtained by setting mj, /i — > OJ ' 'In this limit, the 
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Fig. 17. Fcynman diagrams for lcpton-slcpton scatterings induced by soft gaugino masses m-g 
through the exchange of SU(2)l and U(1)y gauginos \q- The squared amplitudes of these pro- 



cesses arc proportional to rn~ and vanish in the m; 

9 y 



limit. 



supcrsymmetric Lagrangian acquires two additional anomalous global symmetries 
(respectively a i?-symmetry and a PQ-like symmetry) under which, besides the 
SU(2)l and SU(3) C SM fermions, also the gauginos and higgsinos transform non- 
trivially. As a consequence, the EW and QCD sphaleron equilibrium conditions are 
modified with respect to the usual ones: winos, binos and higgsinos are now cou pled 
to the SU(2)l sphalerons, while gluinos get coupled to the QCD sphaleronsM^ 
Therefore, besides the occurrence of non-supcrequilibration (NSE) effects, also the 
pattern of sphaleron induced lcpton-flavour mixing is different from what is ob- 
tained with a naive supcrsymmetrization of the SM case!^ Besides this, a new 
anomaly-free _R-symmctry can be defined and the corresponding charge, being ex- 
actly conserved, prov ides a constraint on the particle density asymmetries that is 
not present in the SMM-^J Nevertheless, in spite of all these modifications, in Ref. [7T1 
it was found that eventually the resulting baryon asymmetry would not differ much 
from what would be obtained with the (incorrect) assumption of SE. Basically, the 
reason for this is that by dropping the SE assumption and accounting for all the new 
effects only modifies spectator processes, while the overall amount of CP asymmetry 
that drives leptogenesis remains the same. 

In SL however, when M > 10 7 GeV the appropriate NSE effective theory not only 
implies profound theoretical modifications, but also results in very large quantitative 
differences. The reason for this is twofold: 

(I) In the NSE regime the leptonic density asymmetries for scalar and fcrmion 
evolve independently, and this implies that the corresponding efficiencies rj s j are 
different. When these different 'weights' are taken into account, the strong cancel- 
lation between the scalar and fcrmion contributions to the baryon asymmetry, that 
is characteristic of SL, gets spoiled, and a non-vanishing result is obtained even 
without the inclusion of thermal corrections!-^ 

(II) While the new symmetries R and PQ that arise in the high temperature ef- 
fective theory after setting m-g, fi — > arc anomalous, two anomaly free combinations 
involving R and PQ can be defined. These combinations, that have been denoted in 
Ref. 1721 as Rb and R x , are conserved in sphaleron transitions, and are only (slowly) 
violated by RHSN dynamics. Thus their evolutions must be followed by means of 
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Table 1. B, L, PQ and B charges for the particle supcrmultiplets. The labels in the 
top row refer to the supermultiplcts L-handed fermion components. The B charges for 
bosons are determined by B(b) = B(f) + 1. 
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two new BE and, since all density asymmetries get mixed by EW sphalerons, these 
equations couple to the BE that control the evolution of B — L. What is important 
is that the sources for Rb and R x are respectively the CP-violating asymmetries 
e s and e s — e/, which are not suppressed by any kind of cancellation. Thus the cor- 
responding density asymmetries Y&r b and Ya.r , remain large during lcptogcncsis, 
and act as source terms for the A Q = B/3 — L a density asymmetries Ya q that are 
thus driven to comparably large values. 

As regards the final values of Rb and R x at the end of leptogenesis, they are in- 
stead irrelevant for the computation of the baryon asymmetry since, well before the 
temperature when the EW sphalerons are switched off, soft SUSY-breaking effects 
attain in-cquilibrium rates, implying that R and PQ cease to be good symmetries 
also at the perturbative level implying that Ya_r b , Ya.r — ► 0. The baryon asym- 
metry is then determined only by the amount of B — L asymmetry at freeze-out of 
the EW sphalerons according to the usual relation B = Jj (B — L). 



6.1. Anomalous and non-anomalous symmetries above Tse 

As discussed above, when the temperature of the thermal bath satisfies the condi- 
tion Eq. (|163[) the appropriate effective supcrsymmctric Lagrangian is obtaine d by 
setting m-g, fj, — > 0, which results in the two new U(l) symmetries R and PQh^I 
The charges of the various states under these symmetries, together with the values 
of the other two global symmetries B and L are given in Table [TJ Note that to facili- 
tate the evaluation of the anomalies in the table we list the charges of the L-handed 
chiral multiplcts, and in particular those of u c , d c , e c that are opposite with respect 
to the charges of the R-handed states u, d, e whose chemical potentials will enter 
the chemical equilibrium equations below. 

Like L, also R and PQ are not symmetries of the seesaw supcrpotcntial terms 
MN C N C + XN C £H U , since it is not possible to find any charge assignment that would 
leave both terms invariant. In Table [1] the charges of the heavy N c supermultiplcts 
have been fixes so that RHSN do not carry any charge. This has the advantage of 
ensuring that all the RHSN bilinear terms, corresponding to the mass parameters 
M, M , B, are invariant, and thus RHSN mixing does not break any internal sym- 
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metry. However, since R(N C N C ) = 0, it follows that the mass term for the heavy 
RHN breaks R by two units.* 

All the four global symmetries B, L, PQ and R have mixed gauge anomalies 
with SU(2)l, and R and PQ have also mixed gauge anomalies with SU(3) C . Two 
linear combinations of R and PQ, having respectively only SU(2)l and SU(3) C 
mixed anomalies, have been identified in Ref. 11231 They axe:' 

R 2 = R-2PQ (165) 

R 3 = R-'iPQ. (166) 

The values of i?2,3 for the different states are given in Table [2] The authors of 
Ref. 11231 have also constructed the effective multi-fermions operators generated by 
the mixed anomalies: 



Oew = n Q (QQQ£ a ) H u il d W 4 , (167) 

QCD 



Oqcd = IT {QQu c d c ) i f . (168) 



Given that the three charges R2, B and L all have mixed SU(2)l anomalies, two 
anomaly free combinations can be defined. The most convenient are B — L ancM-^ 

Rb = \b + R 2i (169) 

whose values are also given in Table [2] The fact that Rb does not contain any 
B — L fragment, ensures that it will not enter in the final computation of the baryon 
asymmetry, and the fact that it is independent of L renders easier writing the BE 
for its evolution. The values of Rb in Table [2] imply that the superpotential term 
N c £H U has charge Rb = 2 and thus is invariant. It follows that RHSN decays into 
fcrmions conserve Rb- In contrast, the soft A term in Eq. (|T2|) responsible for RHSN 
decays into scalars violates Rb by 2 units, and more precisely for N± — > H u £ we 
have ARb = +2, while for N± — > H*£* we have ARb = —2. As regards the heavy 
neutrinos, their mass term violates Rb by two units. Note that this is precisely like 
the case when one chooses to assign a lepton number —1 to the singlet neutrinos 
N. Accordingly, the decays of the heavy Majorana neutrino violate Rb by one 
unit: for N —t £H U , £H U we have ARb = +1 while for the CP conjugate final 
states ARb = —1. All Rb violating reactions have, by assumption, rates that are 
comparable to the Universe expansion rate, and then a specific BE is needed to 
track the evolution of Y&r b . 



* Under R-symmetry the superspace Grassman parameter transforms as 9 — > e la d . Invariancc 
of J d9 8 = 1 then requires R(d8) = — 1. Then the chiral superspace integral of the superpotential 
J d8 2 W is invariant if R(W ) = 2. By expanding a chiral supcrmultiplet in powers of it follows that 
the supcrmultiplet R charge equals the charge of the bosonic scalar component R(b) = R(f) + 1, 
and thus for the fermion bilinear term R (A^./V£) = — 2. 

f With respect to Ref. 123, for definitcness we restrict ourselves to the case of three generations 
N g = 3 and one pair of Higgs doublets N^ = 1, and we also normalize i?2,3 in such a way that 
#2,3(6) =«2,3 (/) + !■ 
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Tabic 2. Charges for the fermionic and bosonic components of the SUSY multiplcts 
under the ij-symmetries defined in Eqs. 11651 1. 11661 1 and 11691 1. Supcrmultiplets are labeled 
in the top row by their L-handed fermion component. 
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At temperatures satisfying condition Eq. (|163[) there is at least one other anoma- 
lous global symmetry, that in Ref. 1721 has been denoted by X- It corresponds to U(l) 
phase rotations of the u c chiral multiplet that, for its fermionic component, can be 
readily identified with chiral symmetry for the right-handed up-quark. In fact, above 
T~2x 10 6 GeV, reactions mediated by h u do not occur and the condition h u — > 
must be imposed, resulting in a new anomalous 'chiral' symmetry. In the SU(3) C 
sector we then have two anomalous symmetries R3 and x, and one anomaly free 
combination can be constructed. Assigning to the L-handed u c L supermultiplet a 
chiral charge x = — 1 this combination has the forrrp-^ 



R\ 



Xu c T 



Ra 



(170) 



where k u c = 1/3. When the additional condition hd — > is imposed, a chiral 
symmetry arises also for the d c supermultiplet. A second anomaly free R x sym- 
metry can then be d efined in a way completely analogous to Eq. (|170p . with 
Kd c = Hu" = 1/3M-^ As regards perturbative violation of R x , this charge inher- 
its the same violation R3 suffers. The soft A term in Eq. (fl"2]) violates R3 by one 
unit, and so do RHSN decays into scalars. Moreover, since N c £H U has an overall 
charge R3 = 1, a violation by one unit occurs also for RHSN decays into fermions. 
Correspondingly, we have AR3 = +1 for the decays N, N* —> H u £, H u £ and 
Ai?3 = — 1 for N, N* —>■ H u £, H^£* . Of course, similarly to Rb, also the evolution 
of R x needs to be tracked by means of one BE. 



6.2. Chemical equilibrium conditions and conservation laws 

Because of the network of fast particle reactions occurring in the thermal bath, 
asymmetries generated in RHSN decays spread around among the various particle 
species, and this can affect directly or indirectly leptogenesis processes. In principle 
there is one asymmetry for each particle degree of freedom. There are however sev- 
eral conditions and constraints that reduce the number of independent asymmetries 
to a few. 

(i) Constraints imposed by reactions whose rates are much faster than the Universe 
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expansion have to be formulated in terms of chemical equilibrium conditions 
for the chemical potentials of incoming fij and final state particles fip: 



X> = X>^- (i7i) 



(ii) Conservation laws that arise when all the reactions that violate some specific 
charge are much slower than the the Universe expansion have to be formulated 
in terms of particle number densities An = n — h and. for a generic charge Q, 
read: 

Q = J2 Qi An * = const ' ( 172 ) 

i 

where Qi is the charge of the «-particle species. We will always assume as initial 
conditions for leptogenesis that all particle asymmetries vanish, and thus we 
will put the constant value of Eq. (|172[) equals to zero. 
(iii) Reactions with rates comparable with the Universe expansion have to be treated 
by means of appropriate dynamical equations. In this case, in order to reabsorb 
the dilution effects due to the Universe expansion, it is convenient to introduce 
as basic variables the number densities of particles normalized to the entropy 
density s so we define the density asymmetries per degree of freedom gf. 

Y Ai = ±^. (173) 

Clearly, u^, An; and Y& i are all related to particle asymmetries. In particular, 
the number density asymmetry of a particle for which a chemical potential can 
be defined is directly related with this chemical potential. For both bosons (b) and 
fermions (/) this relation acquires a particularly simple form in the rclativistic limit 
"nib.f <C T, and at first order in \xtj jT <C 1: 

An b = f T 2 /X6, An/ = |-T 2 u/ . (174) 

Eventually, to solve for the large set of conditions in a closed form one needs to use a 
single set of variables. Here we will take this to be the set {Ya;}, leaving understood 
that the solutions to the constraining conditions are obtained after expressing fit 
and An; in terms of this set through Eq. (|174j) and Eq. (|173|) . 



6.2.1. General constraints 

We first list in items 1, 2 and 3 below the conditions that hold in the whole tem- 
perature range M w <C T < 10 14 GcV. Conversely, some of the Yukawa coupling 
conditions given in items 4 and 5 will have to be dropped as the temperature is 
increased and the corresponding reactions go out of equilibrium. For simplicity of 
notations, in the following we denote the chemical potentials with the same notation 
that labels the corresponding field: <f>= \i^. 
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(1) At scales much higher than Mw, gauge fields have vanishing chemical potential 

IQQI 

W = B = g = 0. This also implies that all the particles belonging to the same 
SU(2)l or SU(3) C multiplcts have the same chemical potential. For example 
<f>(l3 = +i) = (f)(Is = —5) for a field <j> that is a doublet of weak isospin I, and 
similarly for color. 

(2) Denoting by Wr, Br and (Jr the right-handed winos, binos and gluinos chem- 
ical potentials, and by £, Q (£, Q) the chemical potentials of the (s)lepton 
and (s)quarks left-handed doublets, the following reactions: Q + c/r — s> Q, 
Q + Wr -> Q, I + Wr -> £, I + Br -> £, imply that all gauginos have the 
same chemical potential: 

-g = Q-Q = -W = l-i = -B, (175) 

where W, B and g denote the chemical potentials of left-handed gauginos. It 
follows that the chemical potentials of the SM particles arc related to those of 
their respective supcrpartners as 

Q,£ = Q,£ + g (176) 

H u4 = H u4 + g (177) 

u,d,e = u,d,e — g. (178) 

The last relation, in which u,d,e = UR,dR,eR denote the i?-handed SU(2)r, 
singlets, follows e.g. from u c L = u c L + g for the corresponding L-handed fields, 
together with u c L = —Ur, and from the analogous relation for the SU(2)r, singlet 
squarks. 

Eqs. (|176l) - (|178[) together with the vanishing of the chemical potentials of the 
gauge fields and the equality of the chemical potentials for all the gauginos, imply 
that we are left with 18 chemical potentials (or number density asymmetries) that 
we choose to be the ones of the fermionic states. They are 15 for the SM quarks and 
leptons, 2 for the up-type and down-type higgsinos, and 1 for the gauginos. These 
18 quantities are further constrained by additional conditions. 

(3) Before EW symmetry breaking hypercharge is an exactly conserved quantity, 
and we can assume a vanishing total hypercharge for the Universe: 

2/ to t = ^ An fo 2/ 6 + ^ An/ y/ = 0, (179) 

b f 

where ybj denotes the hypercharge of the 6-bosons or /-fermions. It is useful 
to rewrite explicitly this condition in terms of the rescaled density asymmetries 
per degree of freedom {Ya*} defined in Eq. (|173[) : 

Y, (Yaq s + 2Y Aui - Y Adt ) - J2 ( Y At a + Y A e a ) + Y Aku - Y AAd = 0. (180) 
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(4) When the reactions mediated by the lepton Yukawa couplings are faster than 
the Universe expansion rate, the following chemical equilibrium conditions are 
enforced: 

l a -e a + H d + g = 0, (a = e, jU, t). (181) 

For a = e the corresponding Yukawa condition holds only as long as 

T < 10 5 (l+tan 2 /3)GcV, 

when Yukawa reactions between the first generation left-handed SU(2) l lepton 
doublet £ e and the right-handed singlet e are faster than the expansion! I ' I 
Note also that, as discussed in Section 15.71 if the temperature is not too low 
lepton flavour equilibration induced by off-diagonal slepton soft masses will not 
occur. We assume that this is the case, and thus we take the three £ a to be 
independent quantities. 

(5) Reactions mediated by the quark Yukawa couplings enforce the following six 
chemical equilibrium conditions: 



Qi 


- m + H u +.9 = 0, 


(Ui = u, c, t), 


(182) 


Qi 


- di + H d +g = 0, 


(di = d, s, b) . 


(183) 



The up-quark Yukawa coupling maintains chemical equilibrium between the left- 
and right-handed up-type quarks up to T ~ 2 ■ 10 6 GcV. Note that when the 
Yukawa reactions of at least two families of quarks are in equilibrium, the mass 
basis is fixed for all the quarks and squarks. Intergencration mixing then implies 
that family-changing charged-current transitions are also in equilibrium: &l — > 
cl and t L -> s L imply Q 2 — Q3; s L ->• u L and c L -> d L imply Qi — Q 2 . Thus, 
up to temperatures T ^ 10 11 GeV, that are of the order of the equilibration 
temperature for the charm Yukawa coupling, the three quark doublets have the 
same chemical potential: 

Q = Q s = Q 2 = Qi- (184) 

At higher temperatures, when only the third family is in equilibrium, we have 
instead Q = Q3 = Q2 J^ Qi- Above T ~ 10 13 when (for moderate values 
of tan/3) also the r and 6-quark SU(2)l singlets decouple from their Yukawa 
reactions, all intergeneration mixing becomes negligible and Q3 ^ Q2 7^ Qi- 

6.2.2. Above the superequilibration temperature 

We now discuss the condition specific for ranges of temperatures satisfying Eq. (|1G3I) , 
for which the chemical potentials of particle (j) an d of its supcrpartner <f> are not 
equal (NSE) but are related through a (non-vanishing) gaugino chemical potential 
g, as in Eqs. (|176|) - (|178j) . For definiteness, we fix the relevant temperature around 
T ~ 10 8 GcV, and to emphasize that this condition applies only to the NSE regime 
we put the subscript 'NSE' on the numbering. 
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(6nse) Fast reactions induced by the generalized Q CD and EW sphalcron multi-fcrmion 
operators Eq. (fTFfj) and Eq. (|T68)) impljD^Sfc 

3^Qi + ^ Q +iJ u + ff d + 4s = 0, (185) 

i a 

2^g 4 -^(u J +d J )+6.g = 0. (186) 

i i 

At T ~ 10 8 GeV, Yukawa equilibrium for the u quark is never realized. For 
a = e and for the d-quark Yukawa, equilibrium holds as long as T < 10 5 (1 + 
tan 2 (3) GeV< 124 l 117 l and T < 4 ■ 10 6 (1 + tan 2 /3) GeV respectively. Then, for T ~ 
10 8 GeV both condition hold only if tan/3 > 35, while they both do not hold if 
tan/3 < 5. As we will discuss below, in the latter case the Yukawa equilibrium 
conditions get replaced by other two conditions, and thus the overall number of 
constraints does not change. Later in Section 16.2.31 we will present results for the 
large and small tan/3 cases, and since they do not differ much, we omit the corre- 
sponding results for the intermediate case 5 < tan/3 < 35. 

Counting the number of additional conditions listed in items 3-5 and 6nse, we 
have 1 from global hypcrcharge neutrality, 8 from Yukawa equilibrium plus 2 from 
intcrgencration quark mixing, and 2 from the EW and QCD sphalcron equilibrium. 
This adds to a total of 13 constraints for the initial 18 variables, meaning that 5 
quantities must be determined from dynamical evolution equations. These quanti- 
ties can be chosen, for example, as the density asymmetries of the three fermionic 
lepton flavours Y A i a , of the up- type higgsinos Y A ^ and of the gauginos Y A g, given 
that the last one allows to relate Y A i a and Y A ^ to the densities asymmetries of 
the corresponding supcrpartners. This choice would be a natural one since these 
are the density asymmetries that 'weight' the various interactions entering the BE 
for SL. However, the EW and QCD sphalcron reactions Eq. (|167j) and Eq. (|168j) 
imply fast changes of these asymmetries. A much more convenient choice will be 
to use appropriate linear combinations of the various asymmetries corresponding 
to anomaly free and quasi-conserved charges, where with 'quasi-conserved' we re- 
fer to charges that are not conserved only by the 'slow' RHSN-related reactions. 
These quantities can be identified with the three flavoured charges A Q and with the 
two charges Rb and R x discussed in the previous section. In terms of the rescaled 
density asymmetries (asymmetry abundances) per degree of freedom they read: 

Y Aa =GY AQ + Y, (Y Auz + Y Ach ) - 3 (2Y Ala + Y Aea )-2Y A - g , (187) 

i 

Y ARb - -6Y A q - J2 ( 13y A„ s -5Y Adi ) 

i 

+ ^(10F A ^ +7y Ae J + 68F Ag + 10Y Ai?d -27^, (188) 

a 
* These equations should be compared with the SE sphaleron conditions l)B.4[) and l)B.5[) . 
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and 



Y A r x = 3 ^ (3 Y Aui - 2 Y Ag _) + 1 F Ai? , 3 , (189) 

i 

where, in this last expression, 

Y ARs = -ISY AQ - 3 £ (11 Y Aui -AY Adi ) 

i 

+ J2 ( 16r Af„ + 13F Ae J + 82F A g + 161^ - 14i W ( 19 °) 

The asymmetry abundances of the five charges in Eqs. (|187l) - (|189p define the basis 
Y Aa = \Y Aoi ,Y A R B ,Y A n\ in terms of which the five fermionic asymmetry abun- 
dances Y A ^ a = {Y A £ a , Y Ag , Y A fi }, that are the relevant ones for the SL processes, 
have to be expressed. We will do this by introducing a 5 x 5 A-matrix defined 
according to: 

Y A ^ a = A ab Y Ab , (191) 

where the numerical values of A a b arc obtained from Eqs. (|187[) - (|189j) subjected to 
the constraining conditions listed in items 3-5 and 6nse- Let us note at this point 
that the 3x5 submatrix Ai a b for the lepton asymmetry abundances represents the 
generalization of the A matrix introduced in Ref. [61] A^ b generalizes the Higgs C- 
vector first introduced in Ref. |84j and Agb generalizes the C- vector for the gauginos 
first introduced in Ref. [7TJ As regards the asymmetry abundances for the bosonic 
partners of £ a and of H u , they are simply given by: Aj b — 2 {Ag a b + Agb) and 

A Hnb = 2[A s +Ag b ^ 



6.2.3. Additional conditions from Yukawa reactions 

(7-1) Case I: Electron and down-quark Yukawa reactions in equilibrium. 

If the down-type Higgs VEV is relatively small Vd <S v, the values of the electron 
and down-quark masses are obtained for correspondingly large values of the hd 
and h e Yukawa couplings. For v u jv c i = tan/3 ^ 35 we have a regime in which at 
T ~ 10 8 GeV, that is well above the NSE threshold Eq. ([TM]) . both h d and h e 
related reactions are in equilibrium. Since w-quark Yukawa reactions are never 
in equilibrium, in this case, only eight Yukawa conditions in Eqs. (|18ip - (|183[) 
hold. Solving for the density asymmetries Y At f, a = {Y A £ a , Y A g, Y^^ } in terms 
of the charge-asymmetries Y Aa = \Y Aa , Y A r b , Y A r x } subject to the constraints 
in items 3 to 5 and 6nse, yields 



A = 



t 

1 



9 x 827466 



788776 


38690 


38690 


-56295 


41931 \ 


38690 - 


-788776 


38690 


-56295 


41931 


38690 


38690 - 


-788776 


-56295 


41931 


41913 


41913 


41913 


124281 


12798 


102411 - 


-102411 - 


-102411 


108108 - 


-335907/ 



(192) 
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(7-II) Case II: Electron and down-quark Yukawa reactions out of equilibrium 

If Vd is not much smaller than v u , resulting in tan/3 ^ 5, then both h e and 
hd are sufficiently small that at T ~ 10 8 GeV the related Yukawa reactions 
do not occur. In this case we have to set hd, h e — > and the corresponding 
two Yukawa equilibrium conditions in Eqs. (|181l) - (|183p do not hold (on top 
of u-quark Yukawa reactions which are never in equilibrium). However, two 
conservation laws replace these conditions. h e — > implies that we gain a 'chiral' 
symmetry for the right-handed fcrmion and scalar electrons, ensuring that the 
total number density asymmetry An e + Ang is conserved. As usual, we assume 
that the constant value of this quantity vanishes, which in terms of the rcscalcd 
density asymmetries per degree of freedom implies: 

2 



Y 



Ae 



Y 



As 



0. 



(193) 



For the right-handed down quark we could define an anomaly-free charge com- 
pletely equivalent to Yar hi Eq. (|189[) but, given that in this regime all the 
dynamical equations are symmetric under the exchange u <-* d, it is equivalent, 
and much more simple, to impose the condition 



YAd = Yau ■ 



(194) 



The net result is that, with respect to the previous case, the total number 
of constraints is not changed, and again five quantities suffice to express the 
rescaled density asymmetries for all the fields. For the 5 x 5 A matrix defined 
in Eq. (|191j) we obtain: 



.4 



1 



9 x 162332 



/ -210531 21573 21573-12414 12483^ 
8676-165529 -3197-17958 29709 
8678 -3197 -165529 -17958 29709 
7497 7299 7299 23634 4833 

V -11322 -18477 -18477 23940-74385/ 



(195) 



6.3. Boltzmann equations for R-genesis 

In order to render clear the role played by the new charges ARb and AR X and 
by NSE effects, in this section we introduce a simplified set of BE including only 
decays and inverse decays of RHN and RHSN. In this approximation the evolutions 
of the number density of the heavy states normalized to the entropy density s is 
obtained from Eqs. (|106[) and (|107[) by retaining only the two reactions rates 7^ and 
7jv- In writing down the evolution equations for the five charges Y& a , Yar b , Yar x 
it is convenient to introduce a special notation for the scalar and fcrmionic asym- 
metry abundances (per degree of freedom) normalized to the respective equilibrium 
abundances Yf q 



2Y f eq 



15 . 

47r a g„ ' 




^As.A/ = 


Yas,A/ 
\req 



(196) 
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Using Eqs. ([T75|> and (|T74"|) together with Eqs. (fT75|) and (|T771) it is then easy to 
verify that 

VAf,AH lL = ^Ae,AH u "*" " A 9 ■ (197) 

Retaining only decays and inverse decays, the BE for the flavour density asymme- 
tries read: 



>a q =-4W M--2 f + (y^ + y A# J -§- + (W + w„ 



/ 4 



S,CK 



-4 (*) I ^ - 2 1 -f + [y Ala + w„ j -f- + \y Ala + y Ai? j -f , 

where j~*'' a denotes the rate of RHSN decays into scalars (fermions) of flavour 
a, while quantities without a flavour index are understood to be summed over all 
flavours. To an excellent approximation we have 7~ Q = 7~ , and furthermore the 
density asymmetries of the scalars can be expressed in terms of the ones of the 
fermions by means of Eq. (|197p . This yields: 



Y Aa =-e„ (z) U4- - 2 \-f + [y Aioi + y A5u + ^A ? j 2 JV , (199) 



N 4 



..a _|_ _.a 



lAH B =e s (4^- " 2 W "E 0^- +y A5 u +yA ? ) 2 ^^ - W^,(200) 



2 

*A*» = [^ (z) - e> (,)] (^" ~ 2 ) f ~ ^A 5 f ■ (201) 

It is possible, and formally straightforward, to add to these equations the ap- 
propriate terms that allow to extend their validity also in the SE regime, when the 
RHSN masses are below the bound Eq. (|164[) . In order to do this, one must add a 
7~ term characterizing the set of gaugino-mediated reactions with chirality flip on 
the gaugino line that are responsible for processes that equilibrate particle-sparticle 
chemical potentials. § Equivalently 7° characterizes the set of reactions induced by 
the higgsino mixing parameter jj, that enforce the chemical equilibrium condition 
H u + Hd = 0. The thermally averaged rates for these reactions can be written in 
an approximated form as: 



19 _ 9 WB _ M 



n f T n f 



(202) 



§Rcf, l56l included a similar term 7MSSM i n the BE for supersymmctric leptogenesis, corre- 
sponding to the thermally averaged cross section for the photino mediated process e + e ++ e + e 
that was computed in Ref ll25l However, the only contributions that do not vanish in the m^ — > 
limit are those that, like e.g. e^ + e^ <-> e^ + ejj, do not enforce SE. 
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where n e ? is the equilibrium number density for one fermionic degree of freedom, 
while m-g and fj, have to be understood as effective mass parameters in which all 
coupling constants as well as reaction multiplicities are reabsorbed. Extension of 
the validity of Eqs. (|199|) - (|201j) to the SE domain is then achieved by adding the 
following terms to the equations for Rb and R x : 

Yli B = {Y ARB }-y A ~ gl f, (203) 

Yii x = {y ARx } - ly A? 7 f + i (y ASu + y Ag J jf B , (204) 

where the {Iar} above stand for the r.h.s of the corresponding Eqs. (|200[) and (|20ip . 
Note that since the Rb charge of the [i term is RB(H u Hd) — 2, /j, conserves Rb 
and accordingly there is no term proportional to 7° in Eq. (|203|) . Since higgsino 

equilibration involves also the density asymmetry V A ^ we give below the C Hd 
vectors for the two cases: 

Case l:C Hd = — \ — (14237, 14237, 14237, 1260, -3915) , (205) 

827466 v ' y ' 

Case II : C" d = - ]— — (12469, 16768, 16768, 7056, -21924) . (206) 

3 x 162332 



It can be shown that the results of numerically solving Eq. (|199[) and Eqs. 

(|204[) with increasing values of m-g and \x converge to the solutions of the usual BE 

for the SE regime (sec | Appendix Al4 ). 



6.4. NSE regime: R-genesis in a simple case 

The role played by the asymmetries of the two R charges is easy to understand in 
a simple scenario, in which lepton flavour effects play basically no role and thus do 
not shadow the new effects. This scenario is defined by the following two conditions: 

• We assume equal branching fractions for the decays of N and of N± into the 
three lepton flavours, that is the P a defined in Eq. (|49| arc all equal to | 
implying e a = |c 1 and 7 « ~ = ± 7jv ^. 

• We assume the regime described in Case I in which the Yukawa equilibrium 
condition for the electron holds, and thus the three lepton flavours are all treated 
on equal footing (see the 3x3 upper- left corner in the A- matrix Eq. (|192[) ). 
Given the previous condition, it is then useful to define a 'flavour averaged' 
lepton asymmetry as: 



^A^\^ A ^ ( 207 ) 



"Here we assume the UTS scenario discussed in Section 15.21 
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With these conditions, the three equations for the flavour charges Eq. (|199[) can 
be resummed in closed form into a single equation for the B — L asymmetry: 

V^ + W ,+W) 2 ^^, (208) 




yielding a reduced set of just 3 BE. The 3x3 matrix relating {Y A £, Y A g, Yah } to 
the three charge- asymmetries {Y Ab _ l > Yar b > ^ar x } can be readily evaluated from 
Eq. (HMD: 

/ -26348 -6255 4659 \ 

A = 827466 465? 138 ° 9 U22 ' (2 ° 9) 

V -11379 12012-37323/ 

It is now easy to see that in the NSE regime we can rewrite the BE as 

Y Ab _ l =3Y ARx -Y ARb , (210) 

Yar b = e s (z)(^- - 2 J 7 ^ - (y Ae + y A5u + y A? ) 1E±M y A ~2^ 7 (2n ) 

Yar x = [e s (*) - ef (*)] [ 2jg- _ 2 W - y A ~M , (212) 



since the difference in the r.h.s. of Eq. (|210|) gives precisely Eq. (|208p . Eq. (|210|l 
makes apparent how Yah and Y A r b , that in the T — > limit keep having non 
vanishing CP asymmetries, are sources of the B — L asymmetry. Note that the only 
role of the two conditions listed above is simply that of allowing to collapse the three 
equations for A Q into a single one for Ab-l, while maintaining the BE in closed 
form. Therefore the above result is completely general, and in particular it holds 
also when scattering processes are included, and is independent of the particular 
NSE temperature regime (e.g. Case I and Case II) and flavour configuration. In 
short, in the NSE regime the evolution of Ab-l can be always obtained from the 
evolution of 3Ar x — Ar b , and the final value of Y Ag _ L can be equally well obtained 
from summing the values of the flavour charges asymmetries ^ Q Y Aa or from the 
final value of 3Y A r , — Y A r b . The reason why this happens is simple: by using 
the definitions Eqs. ()169|) - (|170|) together with Eqs. (|165|) - (|166j) one obtains that 
3R X — Rb = X u c — |i? — PQ. Of course, only the PQ fragment of this charge 
is violated in RHSN interactions, and from Table Q] we see that this violation is 
precisely the same as for B — L (e.g. for N — > IH U we have A(B — L) = —AL = 
—A(PQ) = —1). Thus, regardless of the fact that B — L, Rb and R x are all 
independent charges, in the NSE regime the BE for 3Y A r — Y A r b always coincides 
with the BE for Y Ab _ l = £ Q Y Aa . 

In this particularly simple case one can take a further step and rewrite the 
density asymmetry y A g and the combination (^ai + Va5 "^ ^ A s) m ^ ne r -h- s 
of Eqs. (|211j) - (|212jl in terms of Y Ab _ l , Y A r b , Y A r by means of the A matrix 
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Eq. (|209p . Replacing Y Ab _ l — > 3Y ARx — Y ARb and using jn = 7^ one obtains: 
3Y ARx = [e* W - e/ (,)] [ ^ - 2 ] if - 9152^+15393^ ^ (as) 




114424 y AflB - 245511 y Aflx 7j? 
827466 T- (2U) 

These two equations show that although 3R X and i?s have the same T = source 
term so that the difference of their asymmetries tends to cancel, their respective 
washouts arc quite different, and such a cancellation will never occur. With a general 
flavour configuration the set of BE cannot be collapsed to just two equations, but 
still the same mechanism is at work: because of the different washouts, the difference 
between 3Y ARx and Y ARb becomes of the same order of these density asymmetries, 
and so does Y Abl . Consequently, one expects that by increasing the washouts from 
weak strengths up to (not too) large strengths, the final value of B — L will increase. 
The numerical results in the next section confirm this picture. 

In the SE regime instead, things proceed in a different way. Eqs. (|203|) - ()204|) show 
that the BE for Y AR and Y ARb acquire new washout terms, that are proportional 
to the SE rates, while on the contrary no analogous terms enter the BE Eq. (|199[) 
for Y Aa or Eq. (|208| for Y Ab _ l . Thus, in the SE regime, Eq. (|2"T0f does not hold. 
One can argue instead that, because of the SE washouts, the roles of A R -l and of 
3Ai? x — ARb get reversed, since now we have 

3 Y ARx - Y ARb = Y Ab _ l + (y A5u + y A5 J if s . (215) 

In other words, since SE reactions conserve B — L but violate the R and PQ charges, 
the only source of asymmetry surviving SE is the Y Abl asymmetry generated 
by thermal corrections. Given that AR X and ARb both contain 'fragments' that 
carry B number, they do not vanish in the SE limit, but are driven to values 
that are proportional to Ab-l- The constants of proportionality are determined 
by the chemical equilibrium and conservation law conditions appropriate for the 
specific regime and, for example, in Case I are given by Y ARb = —\Y Abl and 
Y ARx = ~j § Y Ab _ l . 

6.5. Numerical analysis of R-genesis 

We summarize here some of the numerical results for SL in the NSE regime. They 
are obtained by integration of the BE given in |Appcndix A.5| that also include 
the various scattering processes. The comparative results for the SE case can be 
obtained in two formally different, but physically equivalent, ways. A first possibility 
is that of taking the limit nig, /j, — > 00 in the complete BE (given, for example, in 
their basic form in Eqs. p03j) - (|204p '). A second possibility, that corresponds to 
usual treatments, is to solve only the three equations for the flavour charge density 
asymmetries Y Aa with the corresponding A matrix and C vectors obtained under 
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^m.,,=0.2 eV 




10 10 



10 



Fig. 18. Evolution of Y& B _ L - The solid continuous (red) lines depict the complete results in the 
rrig = fi — ¥ limit. The dashed (blue) lines correspond to the same limit when thermal corrections 
to the CP asymmetries are neglected. The dotted (black) lines gives Y/± B _ L with thermal effects 
when SE is assumed. The picture on the left is for m e g = 0.05 eV and that on the right for 
m cS = 0.20 eV. 



the assumption of SE. For the two cases described in Scction r6.2.3l thc corresponding 
matrices assuming SE are given in |Appcndix B in Eqs. (|B.8|) and (|B.10|) . 

To single out the new NSE effects, for all the results a flavour equipartition 
configuration, with equal flavour branching fractions Eq. (|49[) P a = h is assumed, 
so that flavour effects are basically switched off. In all cases, the heavy RHSN 
mass is held fixed at M — 10 8 GeV, that is above the temperature threshold for 
the validity of the effective theory Eq. (|163|) . The values of the other relevant pa- 
rameters are: A = 1 TeV, 4>a — § and e = -A = 10~ 5 that correspond to a 
resonantly enhanced CP asymmetry in mixing e Eq. (|54[) . This is obtained for 
2 B ~ r ~ 2.6 ( o"i C cv ) GeV. As regards gaugino mass dependent contributions to 
the CP asymmetries from vertex corrections: i v Eq. (f55j) and e 1 Eq. (|56|) . they 
are suppressed by additional powers of A susy /M and thus have been neglected. 
Given the large value of M, they remain irrelevant even in the cases labeled as the 
"trig — > oo limit" , since in practice nig ~ 10 TeV is more than sufficient to enforce 
SE. The results are presented for Case I because, as shown Ref. [72j the differences 
between the situations in which the h eJ i Yukawa reactions are in equilibrium and 
when they are out of equilibrium are rather mild. 

Fig. [18] displays the evolution of Ya b _ l with increasing z = M/T. The solid 
(red) lines correspond to the full results obtained in the mjr, fi — > GeV limit, that is 
when particle-sparticlc supcrequilibrating processes are completely switched off. The 
dashed (blue) lines give the results obtained in the same limit, but when all thermal 
corrections to the CP asymmetries are neglected, and e s = —e/ = e/2. Both pictures 
display clearly that in the NSE regime neglecting thermal corrections in evaluating 
the final values of Ya b _ l is an excellent approximation. The dotted (black) lines 
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give Yis.i with thermal corrections included and under the assumption of SE, that 
in the BE fl203|) - (|204f) corresponds to taking the limit m-g,/! — > oo. The two panels 
are for two different washout strengths m c g = 0.05 eV (left) and m e g = 0.20 cV 
(right) and, as anticipated, they show that stronger washouts result in larger gain 
in the efficiency. 
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Fig. 19. Efficiency factor ijasa function of the washout parameter m e gf for Case I (h e <j Yukawa 
equilibrium) and different values of nig = fi. The red continuous line is for the NSE regime with 
nig = (j, = 100 GeV. The dashed blue line give the result for the same regime when thermal 



corrections are neglected. The red dash-dotted line corresponds to nig 
black dotted line to SE with m;, /i — > oo. 



500 GeV, and the 



Fig. 1191 shows the efficiency r\ as a function of the washout parameter m c ff. The 
red continuous line corresponds to m-g = \x = 100 GeV, and since it is practically 
indistinguishable from the m-g = \i — > case, the evolution occurs in the NSE regime 



in agreement with Eq. (|163|) . The red dash-dotted line corresponds to ■ 



/i = 500 



GeV. In this case we see that SE rates start suppressing the efficiency even without 
attaining full thermal equilibrium. The black dotted line corresponds to the m-g, fi — ¥ 
oo limit of complete SE. The figure illustrates that at T > 10 7 GeV the leptogenesis 
efficiency could be significantly underestimated if SE is incorrectly assumed. The 
size of this underestimation is a fast increasing function of the washouts, and for 
particularly large values of m e g can reach the two orders of magnitude level. Also, 
for m e g ^ 6x 10 -3 eV, the assumption of SE results in a baryon asymmetry of 
the wrong sign. Graphically, one can see this from the fact that at small values 
of m c ff the black dotted and red dash-dotted and continuous lines approximately 
overlap, and there is a change of sign in Y^ B _ L /e around m c g ~ 3x 10~ 4 eV 
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Fig. 20. The final value of Y& B _ L normalized to the SE result Y£ e 
fi for Case I (h e d Yukawa equilibrium) and m c g = 0.20 eV. The red continuous line corresponds 
to varying simultaneously both parameters holding rrig = fi. The blue dashed line corresponds to 
varying only rrig with fi — ¥ oo. The green dotted line corresponds to varying only fi with nig — > oo. 



But around m c s ~6x 10 -3 eV for the red dash-dotted and continuous lines there 
is another sign change. This marks the onset of i?-gcncsis domination; therefore, 
from this point onward, baryogenesis does not proceed through lcptogcncsis, but 
rather through i?-genesis. In the same figure the dashed blue continuous line shows 
the NSE results in the approximation of neglecting all thermal corrections to the 
CP asymmetries. By comparing with the full results (red continuous line) we see 
that for TO e ff ^ few x 10~ 2 cV thermal corrections give negligible effects. Thus, for 
.R-genesis, the zero temperature approximation yields quite reliable results. 

Fig. [20] displays the value of Y^ _ (labeled just as Ya b _ l for simplicity) as a 
function of different values of rrig and /j,, normalized to Y£ E _ that is the final value 
of the asymmetry obtained assuming SE. In order to enhance the impact of the new 
effects, the washout parameter has been fixed to a rather large value m e ff = 0.20 eV. 
The red continuous line corresponds to varying simultaneously both SE parameters 
keeping their values equal: rrig — ii. We see that for rrig = \x ^ 1 TeV the amount of 
B — L asymmetry produced by SL can be up to two orders of magnitude larger (and 
of the opposite sign) with respect to what would be obtained in the usual approach 
with SE. SE effects start suppressing the asymmetry around rrig = \x ~ 1 TeV. The 
asymmetry then changes sign around 3 TeV, that marks the transition from the R- 
gencsis to the leptogenesis regime, and eventually around 5 TeV SE reactions attain 
complete thermal equilibrium and Ya b _ l /Y£ e _ —$■ 1. 

The BE (|203p - ()204[) are general enough to allow to study what would happen 
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Fig. 21. Final values of the charge density asymmetries as a function of mgr = /x for Case I (h e ^d 
Yukawa equilibrium) and m e g = 0.20cV. Thick red line: Y& B _ L \ thick blue line: 3Yah x - ^Afl B ; 
thin dashed purple line: Y&r b ; thin dotted purple line: Y&r . 

if only one of the two anomalous symmetries U(1)r or U(1)pq were present. The 
corresponding results are also depicted in Fig. 1201 The blue dashed line corresponds 
to the £/(l)fl-theory where nig is varied while U{1)pq is broken." The green dotted 



line corresponds to the alternative f/(l)pQ-theory in which 



co and only /j, 



is varied. These results clearly show that the real responsible of the large effects is 
the i?-symmetry, while the effects of the PQ symmetry remains qualitatively more 
at the level of typical spectator effects. A theoretical justification of this behavior 
is not difficult to find, and we will discuss it in the following section. 

Some important aspects of the transition from i?-genesis (NSE regime) to lcpto- 
genesis (SE regime) are highlighted in Fig. [2T1 which displays the final value of the 
relevant charge density asymmetries as a function of rrig = /i, assuming Case I and 
m c ff = 0.20 eV. The thick solid red line corresponds to Y& B _ L , while the thin solid 
blue line corresponds to 3Ya.r , — Yar b ■ The thin dashed and dotted purple lines 
display respectively Y^r b and Yar x ■ We see that up to rrig = fx ~ 100 GcV we have 
Y\ B _ L ~ 3Yar x — Yar b that is in agreement with Eq. (|210|l . and thus implies that 
baryogenesis occurs almost only via i?-genesis. As the soft SUSY-breaking param- 
eters are increased, SE reactions begin to wash out efficiently Yar b and Y&r but 
the difference 3Ya_r — Yar b still remains of the order of Y/\ B _ L , and i?-genesis still 
gives the dominant contribution to baryogenesis. Around wig — fi ~ 3 TeV all the 



II Note that since fj, breaks both symmetries, the case of the £/(l)ij-theory is somewhat academic. 
We include it to put in evidence the fundamental role of U(1)r in enhancing the baryon asymmetry. 
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charge asymmetries change simultaneously their sign. This is the benchmark of the 
onset of the regime in which leptogenesis dominates. The only relevant source for 
generating the density-asymmetries is now the (opposite-sign) thermally induced 
B — L asymmetry, that is not affected by SE washouts, and that is feeding (small) 
asymmetries into all the other charges. In this regime Y Ar and Y Ar do not have 
anymore an independent dynamics, and can be simply computed in terms of Y Ab _ l 
yielding Y ARb = -\Y Ab _ l and Y ARx = —^ Y Ab _ l . 

6.6. Discussion 

In the temperature regime quantified by Eq. (|163|) all reactions that depend on the 
soft gaugino masses do not occur. In this regime the early Universe effective theory 
includes a new i?-symmetry. In SL, this i?-symmetry is violated in the out of equi- 
librium interactions of RHN and RHSN. In particular, i?-number CP asymmetries 
in heavy RHSN decays can be defined, and constitute important quantities. In fact, 
given that it'-symmetries do not commute with SUSY transformations, it is hardly 
surprising that no cancellation occurs between the i?-number CP asymmetries for 
scalars and fermions. For this reason, a sizable density asymmetry for the R charge 
can develop in the thermal bath, and this asymmetry turns out to play the main 
role for the generation of the baryon asymmetry. 

To keep higgsinos sufficiently light, in SUSY one needs to assume [i ~ rrig, and 
thus when the gaugino masses are set to zero, one must set a — > as well. In this 
limit the effective theory acquires another quasi-conserved global symmetry, that is 
a U(1)pq symmetry of the Peccei-Quinn type. PQ is also violated in RHSN inter- 
actions and thus it also has an associated CP asymmetry. However, since U(1)pq 
is a bosonic symmetry that commutes with SUSY, the same cancellation between 
fermion/boson CP asymmetries occurring for lepton number also occurs for PQ. 
Accordingly, PQ does not play an equivalently important role in the generation of 
the baryon asymmetry. 

In order to make more understandable the previous two remarks, let us start 
from the beginning, by listing the relevant global symmetries of the effective the- 
ory. For simplicity we concentrate on Case I (/i e ,d Yukawa equilibrium). Neglecting 
lepton flavour, that is irrelevant for the present discussion, these symmetries are: 
L, R, PQ, B and X„c . The first three L, R, PQ are violated perturbatively in the 
interactions of the RHSN, and all five symmetries are violated by non-perturbative 
sphaleron processes. In this review, in carrying out our analysis, we have first iden- 
tified the anomaly free combinations of the five charges, that are B — L, Rp and 
R x , and then we have written down the BE to describe their evolutions. Here, we 
want to sketch a different procedure. We first write a set of evolution equations for 
the five anomalous charges, that have the form: 

Y Aa = §a q + Sa q + S££ ■ ( 216 ) 

In this equation § represents the source term for Y Al S is the (s)neutrino- related 
washouts with all density-asymmetries and signs absorbed, and S represents 
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the non-perturbative EW and/or QCD sphaleron reactions that violate Aq. The 
latter are reactions of type (i) discussed in the introduction in Section 13.11 and 
in this Section in 16.21 that is fast processes, that eventually will be convenient 
to eliminate in favour of chemical equilibrium conditions. Now, given that B and 
X u c are good symmetries at the perturbative level, they have no CP-violating source 
terms: Sa b , §a x = (they also do not have perturbative washouts, and 9a b Sa x = 
too). The only source terms thus are §Ax,, §a pq and §a h - However, as we already 
know, in the T — > limit, for S& L we have a cancellation between the fermion and 
scalar contributions: S A + S A — > 0. This straightforwardly implies that S A + 
S A —} too, since the RHSN processes contributing to the CP asymmetry for 
PQ are the same as for L: they are simply multiplied by the appropriate PQ 
charge that is, however, the same for fermion and scalar final states. For the R 
charge we have instead § Ah oc Rf ■ § A + R s ■ § A , where Rf iS are respectively 
the overall i?-chargcs of the fermion and boson two particle final states, and thus 
satisfy R s = Rf + 2. We then straightforwardly obtain that in the T — >• limit 
the i?-charge source term does not vanish, and is given by §a r —* 2§ A . Fast 
in-equilibrium sphaleron processes enforce equilibrium conditions between particle 
densities carrying R charge, and those carrying a B and L numbers and, as a 
result, baryon and lepton asymmetries roughly of the same order than the R charge- 
asymmetry develop. Eventually, with the decreasing of the temperature, gaugino 
mass related reactions will start occurring with in-equilibrium rates erasing any 
asymmetry in the R charge. It is important to notice that when the _R-symmetry gets 
explicitly broken, generalized EW sphalerons reduce to the standard EW sphalerons, 
and sphaleron induced multi- fermion operators decouple from gauginos, and reduce 
to their standard B + L violating form.** Since gaugino mass reactions as well as all 
other MSSM processes conserve B — L, the asymmetry initially generated through 
-R-genesis will remain unaffected. 

Now that we have identified where the large density asymmetries come from, 
we can complete our procedure by constructing suitable linear combinations of the 
five equations (|216|) for which the sphaleron terms $ NP cancel out. Since there are 
only two such terms, 9^/ and Sqcd, we can construct three linear combinations 
in which only processes of type (iii) enter. These are the BE for the three anomaly 
free charges that have been discussed at length in Section 16.11 The equilibrium 
conditions enforced by 9^^/ and 9ncr) have to be imposed on the system, and to 
obtain the BE in closed form, the various density-asymmetries appearing in the 
washout terms 9 must be rotated into the densities of the anomaly free charges by 
means of an appropriate A matrix. 



**Here we concentrate on the role and fate of the ij-symmetry. However, given that eventually 
also the PQ symmetry gets explicitly broken, higgsinos decouple from sphalerons as well. 
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7. Soft Leptogenesis Testability and Variations 

As is well-known, leptogenesis models are plagued with the undesirable feature that 
their experimental verification is very difficult, and in minimal scenarios it appears 
to be impossible, at least in the light of foreseeable experimental tests. This is 
because, to establish leptogenesis experimentally, we would need to produce the 
heavy states responsible for the generation of the lepton asymmetry, and measure 
the CP asymmetry in their decays. With the dawn of the LHC, this issue has been 
more pressing than ever, since new states with masses of the order of the TeV could 
become for the first time accessible. However, in the most natural scenarios, the new 
states relevant for leptogenesis lie at a scale that is several orders of magnitude above 
the TeV. With some severe fine tuning, masses light enough to fall within the energy 
range accessible at the LHC could be accommodated. However, in this case to keep 
the light neutrino mass scale within the experimental limits the Yukawa couplings 
of the heavy states must be extremely tiny, preventing again any possibility of direct 
production. The possibility of indirect verifications of leptogenesis, for example by 
pinning down the whole set of the eighteen parameters of the seesaw, and from this 
deriving a prediction for the baryon asymmetry, is also not viable. This is because 
only half of the seesaw parameters are (in principle) accessible at low energy, while 
the values of some of the remaining (unmcasurable) high energy parameters are 
crucial for leptogenesis predictions. 

As regards SL, it is clear that the discovery of SUSY at the LHC can be regarded 
as a basic condition to keep considering this scenario as a possible explanation of the 
cosmic baryon asymmetry. However, in spite of an energy scale that is intrinsically 
much lower than that of standard leptogenesis, and that could even fall within the 
range of energies accessible at the LHC, with respect to the possibility of direct ex- 
perimental verifications SL is in no better shape than standard (or supcrsymmetric) 
leptogenesis. Even if the RHSN mass is low enough, direct production of RHSN faces 
the same no-go issue of extremely suppressed couplings. Indirect evidences could in 
principle come from measurements unrelated to the neutrino sector because, as it 
has been thoroughly discussed in Section [21 SL depends also on soft SUSY-breaking 
parameters that are not directly related with the seesaw. These parameters could 
in principle be measurable through their effects on other low energy observables 
like 1 LFV lepton decays, or the electric- dipole moments of the charged leptons, that 
receive contributions from the complex phases of the soft SUSY-breaking terms. 
However, Ref. 11261 that addressed specifically this issue, found that all the related 
effects arc much smaller than other MSSM contributions, and therefore unobserv- 
able. One can thus conclude that, much alike standard leptogenesis and supcrsym- 
metric leptogenesis, the simplest SL scenario based on the supersymmetric type I 
seesaw, and on the related soft SUSY-breaking terms, also escapes the possibility 
of experimental verification. 
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7.1. Variations of soft leptogenesis 

The exceedingly strong suppression of the production rates of relatively light RHSN 
states is clearly the direct consequence of their gauge-singlet nature, that leaves the 
tiny Yukawa interactions as the only mechanism for their production. In order to 
obviate this problem one can assume that these states are non-singlet under SU{2)t j1 
or under some new gauge symmetry, so that they could be produced at colliders 
through the corresponding gauge interactions. However, in following this approach, 
one has to be very careful because fast gauge scatterings could potentially keep the 
RHSN in complete thermal equilibrium, and/or RHSN annihilation through gauge 
boson channels could leave a too small fraction of out-of-equilibrium decaying states. 

In Refs. QUI ES the MSSM with the addition of an SU(2) L triplet of scalars 
with non-zero hypercharge (as required for the type II seesaw) was considered, 
and the possibility of SL at a low scale was investigated. Fast annihilations of the 
scalar triplets through gauge interactions keep the triplet abundance very close to 
its equilibrium value, and strongly suppress the final lepton asymmetry. However, 
it is argu ed th at for a triplet mass scale 10 3 — 10 4 GeV leptogenesis could still be 
successful! I Due to the small neutrino Yukawa coupling ~ 10~ 6 all low -energy 
LFV processes like fi — > ej or \i — > 3e remain strongly suppressed! ' ' On the 
other hand, Tevatron and LHC have the potential to produce these triplets, and a 
marked signature will be the decay o f the doubly charged component of the triplet 
to lepton and Higgsino pairs! I I 

A supersymmctric seesaw model associated with an extra U(l)' gauge symmetry 
spontaneously broken at the TeV scale has been studied in Ref. 11331 and shown 
to be a viable option for the generation of the cosmological baryon asymmetry 
via the SL mechanism. Such a scenario leads to testable predictions in colliders, 
through the production of Z' and their subsequent decays into RHSN. The RHSN 
will further decay to final states with pairs of same-sign leptons (slcptons) and 
charginos (charged Higgses) through N — N* mixing and CP violation. 

SL in the inverse seesaw model was considered in Rcf. 11341 This is another 
interesting possibility with the potential of being verified experimentally. In this 
scenario the lightness of the active neutrinos is not associated with tiny neutrino 
Yukawa couplings, but with a small dimensional parameter that breaks U(1)l- The 
unsuppresscd Yukawa couplings together with a low mass scale for the RHSN can 
result in a relatively large mixing with the SM leptons, and through this mixing 
direct production and detection of the RHSN at the LHC become possible. 

Ref. ll35l put forth the more speculative idea of implementing SL in a warped five 
dimensional scenario. It was shown that, within the context of extra dimensions, the 
condition of out-of-equilibrium decay and the phenomcnological constraints on the 
neutrino mass can be both satisfied in a natural way, and that all necessary elements 
needed for SL to predict a correct value for the baryon asymmetry can be obtained. 
While the specific SL mechanism of this model does not seem to be easily verifiable 
experimentally, the general idea of extra dimensions could potentially be probed 
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at LHC through the production of the Kaluza-Klein excitations. An experimental 
confirmation of this scenario would certainly increase the phcnomenological interest 
of SL in the context of extra dimensions. 

Other interesting alternative models which utilize soft-SUSY breaking terms 
to realize leptogenesis at a low scale have been considered in Rcfs. 11361 11371 11381 
It remains to be seen if any of these models can yield some clear experimental 
signature. 

8. Final Remarks and Conclusions 

The matter-antimatter asymmetry of the Universe and the experimental confirma- 
tion of tiny but non- vanishing neutrino masses are two among the very few evidences 
of physics beyond the Standard Model. The type I seesaw can elegantly explain the 
strong suppression of the neutrino mass scale, and through the leptogenesis mecha- 
nism can also provide a natural explanation of the cosmic matter-antimatter asym- 
metry. Leptogenesis can be quantitatively successful without any fine-tuning of the 
seesaw parameters, yet, in the non-supcrsymmctric seesaw framework, a fine-tuning 
problem arises due to the large corrections to the mass-squared parameter of the 
Higgs potential, that are proportional to the heavy Majorana neutrino masses. The 
supersymmetric version of the seesaw has the virtue of stabilizing the Higgs mass- 
squared parameter under radiative corrections, but at the same time also introduces 
a serious tension between the lower limit on the seesaw scale that follows from the 
requirement of successful baryogencsis, and the upper bound on the reheating tem- 
perature that must be satisfied to avoid an overproduction of gravitinos. 

However, supcrsymmctry (SUSY) has to be broken, and this yields the possibility 
that leptogenesis could proceed through right-handed sneutrino decays, thanks to 
the new sources of CP violation from complex phases in the SUSY-brcaking sector. 
In this scenario the lcptogcncsis scale is naturally lowered and successful baryogen- 
esis can be obtained anywhere in the temperature range 10 4 GeV < T < 10 9 GeV. 
Accordingly, the tension with the gravitino problem gets generically relaxed and, 
in the lowe r temp erature window, is completely avoided. This scenario, termed soft 
leptogenesi a ' ' (SL) has been the subject of this review. 

As discussed in Section O SL is plagued by the problem of a congenital low 
efficiency, that is related to the cancellation between the asymmetries produced in 
fcrmions and bosons carrying lepton number. It should be stressed that it is the fact 
that lepton number commutes with supersymmetric transformations (that is that 
scalar and fermionic members of the lepton supcrmultiplcts have the same lepton 
charge) that plays the crucial role in enforcing this cancellation. Finite temperature 
corrections break SUSY spoiling the cancellation between the scalar and fermionic 
CP asymmetries, and can eventually rescue SL from a complete failure. 

The basic mechanism of SL was reviewed in Section |3l To highlight the role of 
thermal factors and of the different types of CP asymmetries, in this section the 
simplifying assumptions of a single lepton flavour and of equal density asymmetries 
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for the lepton and slepton states (supcrequilibration) were adopted. 

When the SL CP asymmetries are dominated by the contribution from mixing, 
the RHSN have to be highly degenerated, and in this situation quantum effects 
can become important. These issues have been reviewed in Section [4j In the strong 
washout regime quantum effects can enhance the absolute value of the final asymme- 
try and also induce a change of its sign. However, altogether, the allowed parameter 
space for SL to be successful is not modified substantially. 

Issues related with lepton flavour effects have been addressed in Section [5) Given 
that SL can only occur at temperatures low enough that all the three lepton flavours 
are resolved by their fast charged lepton Yukawa interactions, the inclusion of flavour 
effects in SL studies is mandatory. We have seen that enhancements of the produced 
asymmetry up to factors ~ 10 3 are possible when flavour effects are accounted for, 
and this is sufficient to avoid the need of additional enhancements from resonant 
conditions. Thus, the RHSN do not need to be highly degenerated, and a natural 
scale for the sneutrino mixing parameter B ~ msuSY is allowed. 

In Section |6] we discussed the recently discovered possibility that baryogenesis 
could proceed through _R-gencsis!-^ If the RHSN mass lies above M ~ 10 7 GeV, a 
new scenario different from SL must be considered. In this scenario the asymmetry 
is first generated in a new R charge that is conserved in the effective theory when 
all gaugino-related reactions can be neglected, which is the case when T > 10 7 GeV. 
This asymmetry is then transferred in part to baryons via generalized clcctroweak 
sphaleron interactions. Given that i?-symmctries do not commute with SUSY trans- 
formations, the scalar and fcrmionic members of the lepton supermultiplets have 
different i?-charges, and no cancellation between the i?-asymmetries produced in 
fcrmions and bosons occurs. Thus, in the high temperature window where SL is 
replaced by i?-genesis, a sizable baryon asymmetry can be generated regardless of 
thermal effects. 

In conclusion, SUSY allows for baryogenesis via leptogenesis to occur at any 
temperature from somewhat below the GUT scale down to the TeV. Above T ~ 
10 9 GeV baryogenesis can occur through the supersymmetric version of standard 
leptogenesis, although this possibility is disfavoured by considerations related to 
gravitino overproduction, whose decays would affect Big Bang nucleosynthesis. In 
the intermediate temperature range 10 7 GeV < T < 10 9 GeV, where the gravitino 
constraint gets relaxed, baryogenesis can occur through i?-genesis, which is a highly 
efficient mechanism in which the production of asymmetries does not rely on thermal 
effects. In the lower temperature window T < 10 7 GeV where the gravitino problem 
is usually evaded, the usual SL mechanism can take place, with an efficiency that 
is suppressed by the scalar-boson asymmetries cancellation, but with possible large 
enhancements from flavour effects. Finally, carefully constructed variations of SL 
can allow for a scale as low as the TeV and, as was briefly reviewed in Section I7~T| 
in some of these cases experimental verifications are also possible. 



September 7, 2011 0:24 Review: Lcptogcncsis from soft supcrsymmctry 
breaking 



Soft Leptogenesis 83 



Acknowledgements 



This work is supported by USA-NSF grant PHY-0653342, by consolider- 
ingenio 2010 program CSD-2008-0037, by CUR Gcncralitat dc Catalunya grant 
2009SGR502 and by MICINN grant FPA2010-20807. CSF would like to thank Juan 
Racker for delightful collaboration and discussion in leptogenesis. He is also grate- 
ful to CNYITP, Stony Brook for the generous support. EN would like to thank 
Diego Aristizabal, Enrico Bertuzzo, Pasqualc Di Bari, Sacha Davidson, Guy En- 
gelhard, Ferruccio Feruglio, Yuval Grossman, Marta Losada, Luis Alfredo Muhoz, 
Yosef Mr, Jorge Noreha, Marco Peloso, Juan Racker and Esteban Roulet for fruitful 
collaborations in leptogenesis, and for the pleasure of working with them. 



September 7, 2011 0:24 Review: Lcptogcncsis from soft supcrsymmctry 
breaking 



84 Ghee Sheng Fong, M.C. Gonzalez- Garcia and Enrico Nardi 

Appendix A. Boltzmann Equations 

Appendix A.l. General Boltzmann equations 

Our Universe is very well described by a spatially homogeneous and isotropic metric 
known as Robertson- Walker (RW) metric 



ds 2 = dt 2 -R 2 (t) 



dr 2 

i 2 i Q 2 i 2-2 

'-r at) + r sin 



1 -k 



(A.l) 



where (i, r, 9, (f>) are comoving coordinates, R(t) is the cosmic scale factor, and k = 
0, +1, —1 describe spaces of zero, positive, or negative spatial curvature, respectively. 
For a general process a+b+... -<->• i+j + .... in the RW space, the Boltzmann equation 
(BE) for the phase-space distribution of the particle species a can be written as: 

t- H Mmi=-^k clM - (A - 2) 

where 

|2 



|2 



c[f a ] = -i^ K 3 ::: |m (ab... -> ij...)\ 2 u b ... (i + m u) (i + ,&/,•) ... 

5a 



- |M (*j...->o6...)r/i/i-(l + »/«/«) (1 + W6) -J • (A.3) 

In the above p a is the number of spin degrees of freedom of particle a and 

Aj?;;- = / dn & ...dn,dn J ... (2^) 4 5^ (p a + Pt + ... - Pt - Pj - ...) , 

d n x = ^ . (A.4) 

(2irf 2E X 

In Eq. (|A.3|) . |M a 6..._>i_ 7 -...| is the squared amplitude summed over initial and final 
quantum numbers (spin states and gauge multiplicity) and f x is the distribution 
function of x with r/x — ± if x is a boson or fermion respectively. The factors 
(1 ± f x ) are known as Pauli-blocking (for x being fermion) and Bose-enhancement 
or stimulated emission (for x being boson) factors, respectively. In Eq. IA.21 the 
Hubble expansion rate of the Universe H in the radiation-dominated era is given 

by 



R 3 V 5 M p i' y ' 

where M p i = 1.22x 10 19 GeV is the Planck mass, 5* is the total number of relativistic 
degrees of freedom (#* = 228.75 for MSSM). 

Using the definition of the number density in terms of the phase space distribu- 
tion 

d 3 p 

{2nf- 



9<> I 7^33"/^ ( A - 6 ) 
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and upon integration by parts, the BE (|A.2[) can be rewritten in the form 



—^ + 3Hn a = - ^ [ab... «■ ij...} , 

b,...i,j,... 


^ab... 


\M(ab...^i 3 ...)\ 2 f a f b ...(l + ni .f l )(l 


\M(ij. 


.^ab...)\ 2 fJ 3 ...(l + Va f a )(l + Vb f b 



(A.7) 



where 

[ab... o ij...] = A^;;. |M (aft... -> »j...)f /a/6- (i + mh) (i + »&/,■) ■ 



(A. 



Notice that in writing the BE (|A.7|) , we have implicitly assumed that the right hand 
side of this equation can also be written in term of n a . However, without certain 
approximations, this cannot be done in general and we have to resort to the BE 
(|A.2[) and solve it in term of phase space distribution. In |Appcndix A.2.1 we list 



the approximations that allow us to write the right hand side of Eq. (|A.7[) in terms 
of number densities, and then to use the BE (|A.7[> . 

In order to scale out the effect of the expansion of the Universe, one defines the 
particle abundance i.e. the particle density n a normalized to the entropy density s 
as: 

n = — , (a.9) 

s 
where the entropy density in the radiation dominated era is given by 

s = ^9«T 3 . (A.10) 

Using the conservation of entropy per comoving volume (i.e. s R 3 = constant), 
replacing the time t with the temperature T (in the radiation dominated era t = 
jjj ~ T~ 2 ) and defining the dimensionless parameter 

* = Y' (A ' n) 

where M is any convenient mass scale, the left hand side of Eq. (|A.7[) becomes: 

dn dY dY 

££ + 3Hn a = s^- = sHz^. (A.12) 

at at az 

Regarding the distribution functions in Eqs. (|A.8I) . for particles for which the clastic 

scatterings are much faster than the inelastic scatterings, one can assume that they 

are in kinetic equilibrium and have either Fermi-Dirac distribution (for fcrmions /) 

or Bose-Einstein distribution (for scalars s) given respectively by 

where fj,'s are the chemical potentials and the "bar" or "star" refers to the corre- 
sponding antiparticles. The equilibrium distributions /J 9 are defined as those with 

f eq - = f eq , = (A 141 

•'/,/ e Ef/T + 1 > J",s* e E s /T_l y^-^) 
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So for #r <C 1 one can expand the kinetic equilibrium distribution function in #r as 



J fJ - J f ±J f e T 



/. .. = f? ± ft 



T 

•,2 Es/T^s 

T 






w\ 





m 2 } 



(A.15) 



It follows that 



//-/ 7 =2(l-/F)/ff 



Js /s* 

where we have used the identities 



2(l + /| 9 )/ s e9 y 







T 



3 (B S TM»)/T 



/,■ 



(A.16) 

(A.17) 



l-f fl = e^*W 1 f f j, 1 + /.,.. 

Using Eq. (|A.6[I one gets that the difference between number densities of massless 
particles and antiparticles at leading order in chemical potentials is 



n Af 



9_j_ T zV± 



riAs 



*ZiT 3 — 



(A.18) 



if — n-T 
i f 6 t ' 3 T 

Defining the density asymmetries per degree of freedom Fa/,s = ^A/,s/(ff/,«s) as 
the the number density asymmetries per degree of freedom «a/,s/<7/,s normalized to 
the entropy density s, we can rewrite the chemical potentials for massless fcrmions 
and bosons as: 



fx f _ 8n 2 g t 



T 



15 



Y 



Y 



A/ 



A/ 
yeq ; 



, Ms 

'T 



15 



r 



As 



>As 



(A.19) 



where Y f e<? = ^§- and F/? = r ^-. 



Let us introduce the following shorthand notation: 

^,.,,..(...)EA^|MK..^y...)| 2 (...), 
F^— .{...) = Al: |M (ij... -+ ab...)\ 2 (...) , 



(A.20) 



where (...) denotes some function to be integrated over. Note that CPT invariancc 
implies F a b...ij... (...) = F^ — —■ — (...). The thermally averaged reaction densities can 
be defined as 

7 (oft... -> ij...) = F ab ... ij ... x f?ff... (1 + w /?«) (1 + mff) ..., (A.21) 

where we have used the equilibrium distribution functions with vanishing chemical 
potentials Eqs. (|A.14|) . 

Neglecting Pauli-blocking and Bose-enhancement factors and assuming that 
all the particles follow Maxwell-Boltzmann distribution / = e~ E / T , and that 
|M (ab... — > ij.--)\ does not depend on the relative motion of particles with respect 
to the plasma, Eq. (|A.21[) for the decay N — > ij... reduces to 



7(0 -> ij...) = j(ij... -> a) 



a X 2 (z) 



r„, 



(A.22) 
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where T a is the decay width in the rest frame of a, % q is the modified Bcsscl function 
of the second kind of order q, and n^ q is the equilibrium number density of a: 



9a 



d 3 Pa 

(2tt) 3 



-E a /T _ 9a 



rp3 



(A.23) 



For a two-body scattering ab — > ij, Eq. ([A. 211) reduces to 

T 



j(ab -> ij) 



64tt 4 



j°° dsVi*{s)Xi(£\, (A.24) 



where s is the center of mass energy squared, s min = max[(?7i a + m^) 2 , (rrii + rrij) 2 ] 
and a(s) is the reduced cross section which is related to the total cross section a(s) 
(summing over initial and final spin states) by 

1 



2X 2 (s,m 2 a ,m 2 ) 
a(s) = 2 CT(S) 8ns 



t + 



dt |M(a6-M.7')|' 



with 



A(a, b, c) — y(a — b — c) 5 



46c, 



(A.25) 
(A.26) 



4.s- 



2\2 



(s + m 2 a - m 2 ) 
4s 



ml =F 



(s + m 2 — m 2 ) 2 
As 



(A.27) 



Because of the thermal-statistical nature of the CP asymmetry in SL, a rigorous 
treatment would require the use of the BE for the particle distribution functions 
Eq. (|A.2[) rather than the integrated BE for the number densities. We will describe 
in the next section the derivation of the BE for SL and the approximations required 
to write them in integrated form, while keeping the relevant thermal statistical 
factors. 



Appendix A. 2. Boltzmann equations for soft leptogenesis 

Appendix A. 2.1. Unflavoured Boltzmann equations 

In the rest of this section, unless otherwise stated, we will ignore thermal masses. 
The BE for the RHN abundance Y/v can be written down as: 



Y N = - N o H u £ - [N «-► HJ], - N£ o 
J+ + L 

NQ «-> tu*\ - \Nu ^ I*q] - \Nu f* £*Q 
J+ L J+ L 

- [M ^Qu] + - [Nu <+ 1Q] + -- [NQ <r*lu] + 
2-Fjv tso - --, Jen + 2-F/v 



N£^ Qu 

NQ* <-> tu 



req 



N 



1 



£ eq 

■' J 



N 



req 



N 



1 



Jn 



AP} ' + 4F (1) + AP t 



(2) 



2f; 



(3) 



AR 



(4) 



AT 



i-/ 



A' 



/•eg 

./a 



i-.r 



IN 



(A.28) 



September 7, 2011 0:24 Review: Lcptogcncsis from soft supcrsymmctry 
breaking 



88 Chee Sheng Fong, M.C. Gonzalez- Garcia and Enrico Nardi 

where the time derivative is defined as Y = sHzj^, [ab -H- ij], = [ab «-»• ij] + 
[ab <-> TJ\ , and [ab <-> ij]_ = [ab •(-> ij] — [ab -H> TJ\ and we have further defined the 
following shorthand notations: 

F N {...) = F NSJ r N «(i + f?) (i -/£)(...) , 
f n (...) = F NHut f% (i - n q ) (i + /i?j (...) , 
^ (0) (■■•) - Fsifrftf? (i - /a 9 ) (i + #0 (■■•) 

F« (...) = F^vftft (l + /f ) (1 + #) (.-) 
= ^^Q/^/„ e9 (l + /f)(l + /| )(...), 
if > (...) ee F^fPf? (l + /f ) (l - #) (...) 



i 7 , 



(3) 



(...) = F NtQtt f%f? (i - /a 9 ) (i - /«*) (-) 

f t (4) (...) ee F Nv3Q f%f? (1 - jf) (l - /§*) (...) 

= F Nm f N q r Q q (i-fr)(i-fu q ) (...). 



(A.29) 



The BE for the RHSN abundances Y~ are: 



N± 



N ± <+ #^ 

N±Q o £u 



W± <+ H.J 

N ± u^I*Q 



N ± ^£uQ* - N ± £^u*Q 



N±£ o Qu* 



N±u o £Q - N ± Q* ^iu - N ± Q o £u* 



N ± e o Qm 



N ± t <-> Qu 



7V ± Q^£u 



,(3) 



F k +F h +2F h:+ QF ^ 



JN 



N ± u o £Q 



f~ q i + /-* 



-2 (2F^ + 2F™ + 2F« + F® + 2f£> 



^#V (A.30) 



req 



N = 



J N ± 



where terms of order (e^), have been dropped. The shorthand notations are 

*s ± (•■•) - (*W+ *W*) /£ (i + /f ) (i + /£) (■■■) . 

*g (■••) - ^ ± &q./h (i + /f ) (i + m (i + /a 9 ) (■■■) . 
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F 



F 



F 22± (...) e, F^^fgf^ (1 + /-) (l + /-) (...) 

f£> (...) ^ F s±tQS .f% ± f? (i - /-) (i + #) (...) 
= ^ ± ^/£/r(i + /?)(i-/« 9 ) (••■). 

^ (-) - ^ ± ^q/£/~ 9 (i - /d (i - /a 9 ) (■•■) 

F t ( p (...) = F N±iQ ^r N q jr (i - /^) (i - fu q ) (-) , 
*& } (■■•) - F Ni: *QfS ± f? a - /f ) (i - /q ) (•••) 

= F N± QguJ]y ± jQ (1 — Jl )(! — /«) (•••) 7 

F t ( ^ (...) ee F N±ulQ f N yt? (1 - /f ) (l - /*) (...) 
= ^N ± QeuJN±JQ (1 ~ h ){^~ fu ) (•■•) ■ 
The BE for the asymmetry in the lepton doublets Y&e = (X« — Yj) /2 is 

2Fa^ = ^ [JVi ^ #„*] _ - ^ [H u £ <-> y 

i=± ij 

- [Nl ^Qu}_- [Nu ^lQ]_- [NQ «-► lu] 

»=± ^ 

= 2^ 1 I Niii u i 



(A.31) 






[JV ++ H u £]_ 

*]- 






i=± 



"% 






1 + / 



A', 



1 + /S? 



1 + / 



J _Nj_feq, ■ J Ni f-, _ req\ 

req J I T 1 , j. e g I, 1 ./£ j 



A', 



-2F, 



A 



J Ni 

In ,eq 1 ~ In 

req J I + 1 /.i 

■/JV X -'AT 



<?</ 



(1-/D 






r 



- (// -> -/j?„ > W -> Mff J 



+4^«T (f " f ) /?*" (l + /f) 2 + ^ + W AL = 2 , 



(A.32) 



*Here Y^£ refers to lepton asymmetry abundance in single SU(2)l gauge degree of freedom. 
However, since the amplitude on the r.h.s is summed over gauge multiplicity, we have to multiply 
by a factor of two in the l.h.s. of the BE. 
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where (// -> / 



peg 



/i£ — > fijj ) etc. refers to the term obtained by replacing the cor- 

sub 

refers 



responding / and [i in the preceding term. In Eq. (|A.32[) . y^- 



H u £ «• ij 



to the sum of all possible AL = 2 scatterings H u £ o ij and, if N± exchange in the 
s-channel is involved, the on-shell contributions are subtracted out to avoid double 
counting. The AL = 2 scatterings H u £ -f-> ij with t- and u-channel exchange of N±, 
and the leftover off-shell contribution for s-channel exchange of N± are all collected 
in Wal=2- The details of the subtraction procedure is given in |Appcndix A.2~3 



In the numerical calculation, we neglect Wal=2 since it is subdominant in the SL 
temperature range T < 10 9 GeV. The top and stop scattering term St in Eq. (|A.32|) 
is given by 



Sr- 



-2R 
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(A.33) 



where in the last line (Q — > u) and (Q — 
is replaced by u and Q is replaced by u. 



u) denote respectively terms in which Q 
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The BE for the slepton asymmetry Y A j = [Yj — Yj,) /2 can be written as: 



2Y. 



j = ^2 [n+ *+ hJ\ -J2 [hj o ij 

i=± ij 

Y^ [ [Ni «+ luQ*] _ + [Nil* *+ uQ* 



N o B.,1 
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(A.34) 



where the AL = 2 scatterings -ff M ^ O ij with t- and u-channel exchange of N± and 
the off-shell contribution for s-channel exchange of N± are all collected in Wal=2 
that, as already said, can be neglected in the SL temperature range. 
The term S22 hi Eq. IA.34I from scalar potential terms is given by 



S 22 =£ 2J? 
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(A.35) 
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while the top and stop scatterings term St reads: 



St 



.4p(o) 



f eq 

JN 



1 + ./V 



i-.r 



■' /:■ 



<<l 



-4\F<n+FP 



/jV ,eg 
j?eg.n i 

J/v 



JV 

I -In 
1 



T 



/■eg 



JV 



l + /~' 



T 



/jV -eg 1-/AT /-. _ 
fCq-'Q T -i _ r eg I 1 ./ 1 



eg 



+2F« 

+2(f(°)+F (1) ) 



JV x JN 

fN /% c eq\ , : - /JV ,eg 

•/JV 



*(! /qJ + ^^/q 



T 
T 



/jV req - - In 



fN*'® 



i-r 



/V 



1 



-2F& 



I*L( 1+ fi 



Q 



l-f- 



N 



1 f e « 
'J-'o 



T 



£ H 



(8) 



»=± I 
(9) 



Wjv 

JNj /-, reg\ _ ^ + / JVj ,e 

,eg ^ + /J J 1 req I J 



JVi 



JV,- 



fS) 



IJq 

T 



T 



+4F 



JJVi ,eg _ + ^ JV, 

~ req J J + -I , jeg 

•fjv, i + / JV ! 



1 + /7 



T 




(A.36) 



Appendix A. 2. 2. Approximations: integrated Boltzmann equations 

In order to write the BE in the integrated form of equations for the number densities, 
the following approximations are needed: 

1 + Vafa 



1 + Vafa 



>■<! 



1, 



Vil, 
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T 



0. 



(A.37) 



where a refers to N or N± and i for any other particles. The approximations above 
are equivalent to neglect the chemical potentials in the Fermi-blocking and Bosc- 
cnhanccmcnt factors. In addition, wc need to assume that N and N± are in kinetic 
equilibrium, that is 



f_N 
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(A.38) 



The BE in their non-integrated form have been considered in Refs. 11391 11401 1551 
and it was found that in the strong washout regime the distributions of N and 
N± are close to kinetic equilibrium, and thus Eqs. (|A.38j) represent a very good 
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approximation. However, it should be mentioned that, as discussed in Rcf. 1851 in 
the weak washout regime numerical differences when using the integrated form 
of the BE can be up to one order of magnitude. In any case, by means of the 
approximations (|A.37|) and (|A.38[) on can write a set of integrated BE as 
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Y N ± 



veq 

1 N 
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\ 'n 



N± + 7 JV± 




4 7 f) + 2 7 f)+4 7 f) 



(A.39) 



(27&+27&+27£+7ff+27i 




(A.40) 



where we have defined the reaction densities for RHN decays as jn = ^jv(l) + 
Fjv(l), RHSN decays as 7~ s = F~ s (l), interactions from scalar potential as 7- 
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By using the following approximations 
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that are justified by the fact that the N± mass splitting is small B <C M, we can 
sum up the BE for iV + and iV_ (|A.40|1 by defining Y^ 
up with: 
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(A.41) 



With the same approximations we can write the BE for Y&e and Y.j as follows: 
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(A.42) 
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where we define the CP asymmetries for N± decays as follows: 

K F NiH u t ~ F nMV ^nS 1 ~ ^ H 1 ~ ^, 
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(A.43) 



(A.44) 



(A.45) 



In order to write the BE (|A.42[) and (|A.43[) in a closed form, all chemical potentials 
have to be expressed in terms of \x t j, and then these quantities must be rewritten 
in terms of the lepton and slepton density asymmetries by means of Eqs. (|A.19[) . 



Appendix A. 2. 3. Subtracted 2-^2 scatterings 

Although 2 -<->• 2 scatterings are of higher order 0(F 4 ) with respect to decays and 
inverse decays which are 0(Y 2 ), the CP asy mmetries of the 2 •<-> 2 subtracted rates 
are of the same order than that of the decays] I I and hence cannot be ignored. 



The term 



H u £ «• ij 



sub 



in the BE (|A.32[) consists of the following two terms 



sub 



F S« 



F z 



H u £ -n- ij j. H^Hj 

ub / \Xi + fJ,j fit + &H, 



Sill) 



l + /««e*£ + /Se 



r 



H„ 



T 



-F 



H„ 



T 



T 



f?f%Jl + Tkf?)(l + Tl j f?) 



+ { F H u u 3 -FlJ~) SUh feV e i (1+THf?) (1+%/f) (^/r| + ifc/r f ) . (A.46) 



September 7, 2011 0:24 Review: Lcptogencsis from soft supcrsymmctry 
breaking 



Soft Leptogenesis 95 



and, for the CP conjugate states 

; — i sub 



H u £ <-> ij 



F=-_ -F~-- 



eq ^A fit _ f eq „ZMm. V-R^ 



i-/r^f -/£ e 



T 



x/r/f u (i+r?i/r)(i+^/r) 

-sub / /Uj + /Uj W + Mff. 



-F- 



ff u < 






r 



T 



f! q f e l(i+vifn(i+vjf; q ) 



F-= 



/f /£ (i+^/D (i+^/f ) (%/f | +%■/;« f 



(A.47) 



where 77; = +1, —1 for boson and fermion respectively. From Eqs. (|A.46|) and (|A.47[> 
and using CPT invariance, we obtain 

Ah = T,{ 2 fam ^T n q f e l (i + thJ?) (1 + mfT) 



J2 [h u £ <-> ij 



F H U 



F- 



H u Uj 



sub 



fJ-i 1 jJij „eq neq 



T 



F~ + F-=-_ 

H v^J H u £ij 



Sill.) 



T 



^fPfi (l + r?i/D(l + %7f) 



"^ ' H u tij ^H„iij 



SU V/£(l+^/D(l+%/D(^/rf +%/ff)}-(A-48) 



For i7^^ -f> ij with AT^ exchanged in s-channcl we have 

sub 

(i + ^/f)(i + ^f) 






-Fl 



(././■• 



i'ij 



F H u wA l + fN k )^[ N ^i3 



- F ku 3 + F ^N k (1 + /ivj Br (^ -»• yjj (i + %/D (1 + ^/f ) 



where Br I A^ — > ij \ is the branching ratio for the corresponding process, that sat- 
isfies the unitarity condition 

]T Br (N k -> ij) (1 + mf?) (1 + %/f) = 1. 

y 

In the last equality in Eq. (|A.49|) we have neglected terms proportional to the CP 
asymmetry I F~ — F4^-_ I which arc of higher order ~ 0(Y 6 ). Substituting 

V Hu.£ij H u Hj J 

Eq. (|A.49j) into Eq. (JA.48JI and ignoring the term of order O(e^), we have 



£[*. 



I <H> ij 



Sl.ll".> 



2 V(f~ - -F — )±±1*± 

xfg (1 ~ if ) (l - /f J - ^AL=2 , (A.50) 
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req 

where we have used the identity f^ q f e ~ = 1+ ^^ (1 — ff> q ) ( 1 — /~ 9 ) and 

w AL=2 = -J:(f^ uj +f w Jf^fTf* ^^ (i+^)(i+%/; ? ). (A.51) 



T 



Following the same procedure we also obtain 

ij k J N k 

Xfe l i 1 + fl) C 1 + fZ) " ^al =2 , (A.52) 

where 



i j 



Appendix A. 3. Lepton flavours and lepton flavour equilibration 

The unflavoured BE (|A.42j) and (JA.43J) can be easily generalized to the flavoured 
case. By denoting (s)lcpton flavours by a = e, /i, r we simply have to replace Y A £ ~ ► 
Y A £ a , e - ► £q, 7 - > la (except for the normalization factor of e Q that is still 7jy), 
Wal=2 —> W / al=2 c ^ c - However, before writing the flavoured BE we now take some 
further step. Chemical equilibrium enforced by the top Yukawa interactions implies 

- m + M« = Vh u , -MQ + M« = L l h u ' -Vq + M« = M5 u , (A.54) 

which also yields 

- Mq + Ms = 2 /i^ u - /iff„ ■ (A. 55) 

This equation and the first relation in Eq. (|A.54j) allow to eliminate in the BE the 
chemical potentials of the (s)quarks for those of the Higgs(inos). It is also convenient 
to trade chemical potentials u^ for the corresponding density asymmetries Y Ar p. This 
can be done by using Eqs. (|A.19[) . For / = l a ,H u and s = £ a ,H u we have 

2/i/ Y Af _ 2^ s Y As _ 

-jr = yw =VAf, — = yeq=yAs, (A.56) 

where 2Yf 9 = Y^ q = . \ 5 . With these conventions the flavoured BE can be written 



as: 



2Y Ata =E a +W£ L=2 , 

2Y Al =E a +W£ L=2 , (A.57) 
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E - f> (z)-^- 




f,a 



-J(^Jn f^v i o 1 _ 2 ] _22 



(W + ^gj - \l% (^Ai a + VahJ 



(3)q *iV 9 «)a , 9 „ ; (6)a (7)a „,(5)a Y N tot \ y 

7t ^g + 2 7t + 2 7t + 2 7t + It -TFW 9A£ Q 



Y 



N 



y: 



' N 



I?'" + 7, (4)Q ( 1 + ^) + I?"" + l { t 6)a + ll^-y^-) W u 



(5)a.(7)a. 1 (6)a Y N m 



- I 7t + 7* "I" o7t 

+ 2 7f (y A i-„-y 



9 'i v e, J 



?A£ Q , 



(2y A 5 u - ^ 



(A.58) 



and 



7tv / Y N tot 



E a = eS, (*) -f I -^ - 2 J " ^f" (^~ + ^„) \l a N [ 



3)q 

'at 



1 „a N tot 



2 ,zz r~ q 



2 72 Q 2 ) (^-+2^-^^) 



27}° )a ^ + 2 74 (1)a + 2 7t (2) " + ^7 t (8)a ^f + 2 7t (9)fe I y 
V )Q + 7 t (1) "^ + lf )a + 7 t (9)Q (l + y^\ ) ¥a*. 



(0)a . (1 

7t + 7* 



(l)a . (2) a r N \ L, ,. 

t + 7t y5g ) (^ A g„ - W, 

- 2 7f (y A ?„-VA, 



A£ Q 



(A.59) 



Appendix A. 3.1. Lepton flavour equilibrating interactions 

The off-diagonal soft slepton masses induce lepton flavour violating (LFV) inter- 
actions through the exchange of SU(2) L gauginos A| and U(l)y gaugino Ai (see 
the Lagrangian (|142|l ). and this can result in lepton flavour equilibration (LFE). 
There are two t-channel scatterings £ a P -f-s- IpP* , £ a P *H> £pP and one s-channel 
scattering £ a £% O PP as shown in Fig. QTJ We denote fermions as P and scalars 
as P. For processes mediated by SU(2)l gauginos P = £,Q, H UJ i, while for U(l)y 
gaugino one must include the SU(2)l singlet states P — e, u, d as well. We have in 
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general 



t P O (.fjP 



'l a P <+ £bP 



&OL*~-R ~^r ±± 



£ a P£ B P I "^1 



2 (^lp +^P Hi a + HP 



T T 

ft q jp ( 1+ C)( 1+ /f)' 



= -2F, 



,pe B p I ^^l 



T 



T 



ntf$(}+fg)v-fP): 



T 



£ £* PP l^ a /3 



/;:/£(i-/; 9 )(i+/f). 



2 (ixp + np H a V>e f 



T 



(A.60) 



where the factors of two come from the CP conjugate processes. Each l a l$— gaugino 
vertex involves one element R a p of a unitary matrix. In Eq. (|A.60[) we have explic- 
itly factored out the entries \R a p\ and hence, if we ignore the zero temperature 



lepton and slepton masses, F f 



PigP 



(...), F. 



(...) and F e j, PP (■■•) are flavour 



independent. With the same approximation the distributions f^ q and f~ q are also 
flavour independent thus, from now on, we will drop the flavour index whenever 
possible. For simplicity, we only keep the thermal masses mr and mr of the 
SU(2)l and U(l)y gauginos. With this approximations, we can define the flavour 
independent LFE reaction densities as follows 



7ti,G - *irfp(fle)£7? (1 + ff) (1 + f~ q ) . 
7e,G = F^{9G)f?f? (l + ff) (1 ~ l P q ) , 



7s, G 



F^ PP (ge)ft q ff(l-f P q )(l + f P q 



(AM) 



where G — 2, Y for the scatterings mediated respectively by the SU(2)l and U(l)y 
gauginos and correspondingly g& = 92 or gy. For example, to rack the evolution of 
the abundance of l a for P = £ P = t we need the following terms: 



Y. \^< *» ^] _= - 2I1 ^ ^ 3 E l^a/3 I 



9lX,P,V 



T 



E ^4|,(A.62) 



T 



E [Wc«¥, 

geX-.P,n 



-2U n g 3^|i? a/ 



I j. "T" ^^ ji 



, (A.63) 



[*J 



gn.X,P,n 



- 2n<7f [ 3 £|i^^ + £^< 



(A.64) 
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where we have used the following properties of unitary matrices: 



z2 \ R <*p\ 





(no sum over a), 



(A.65) 



22\R a p\ -3. 

In Eqs. (|A.62[) - (|A.64[) IT is a factor arising from summing over isospin degrees of 
freedom of leptons and sleptons, and for example for the scatterings mediated by 
Aj we have He = 3. Note that since £ a — t a = /j.g, for j3 = a the sum of chemical 
potentials always cancel (soft slepton masses can only induce LFV interactions but 
not superequilibration) . Hence Eqs. (|A.62[) - (|A.64[) simply become 

|2 »!, ~ W- 



£,P,V 

Y, \u<^~^u 



9e£,P,V 



ge£,P,v 



-m n ZY\R a p\ 



-m n gY\R 



a/3 



-6U e ^Y\ R ^ 



T 

2 /% ~ M° 
T 

2 v-i* - w« 



(A.66) 
(A.67) 
(A.68) 



9t£,P,v P¥=<* 

Similarly, for processes mediated by SU(2)l gauginos we have the scatterings 
with P = Q,H u _d for processes mediated by U(l)y gauginos P — Q,H Ut d,e,u,d. 
The corresponding term that has to be added to the BE for Y&.i a is 

(*k„) = - J2 n G {lg + lg + 7?) E l^^l 2 0^ o - «A?J • ( A - 69 ) 

G=2.Y /3^a 

with ni = 42 and riy = 38. The analogous term for Y A j can be obtained from 
Eq. (|A.69[) by exchanging the asymmetry labels £ a — > £ a , £p — > I p. Given that the 
LFV factor \R a p\ is explicitly accounted for in the terms above, the reduced cross 
sections for the LFE interactions can be defined in a flavour independent way: 



a t\ K s ) = ^z 



°% (*) 



*. 00 = ^ 




1 In 



m~ + s 



Ac; 



Ag 



(A.70) 



where the gaugino thermal mass is m~ = (9/2) g^T 2 , and quantum statistical 
factors have been neglected. 



Appendix A. 4. The superequilibration regime 

Superequilibration (SE) is defined by chemical potentials condition 



(A.71) 
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where (0, <fi) are the fermion and boson components of a generic supermultiplet. In 
the SE regime, the BE for the RHN and RHSN abundances are still given respec- 
tively by Eqs. (JA.39I) and (|A.4L|) . However, Eqs. (|A.54|) combined with Eq. (|A.71j) 
yields 

MQ - Mu = Mq - /"« = -Mff„ = -M-ff«) (A. 72) 

and, as we will now see, this allows us to sum up the two BE (|A.57[) into a single BE 
since with SE the relation between scalars (s) and fermions (/) density asymmetries 
is Yas = 21a/ which implies Y A j = 2Ya£„ and Y A ^ = Yah u /2- 

Summing up the two equations (|A.57|) and including LFE effects we obtain: 



Yu~= [E 



E a ) +[ Y Ai 

SE 



LFE 



(A.73) 



(y A( 



Y a j ) and we have dropped the AL = 2 off-shell scatter- 



where Ya£ q 

ing terms Eqs. (|A.51|) and (|A.53|) which in SL are always negligible. The LFE term 

is obtained by summing Eq. (|A.69|) to the analogous term (Y A j )lfe which gives: 



\A h ), (A.74) 
while summing up Eq. (|A.58j) and (JA.59J) under the assumption of SE yields: 



[E a + E a \ 



SE 2 



iVtot 

Y~ q 

, N 



1 



N 



TjV , (3)a 

2 'n 



lY, 



N t , 



(7 t (1) 



2 y* 



li 



(2)a 



(4)c 



'o (0) , (3)a 

27* + It 



-li 
Y N 

1 N 



li 



(6)a 



7< 



722 



(7)a 



+ 7* 



„(5)c 



(8)a 



(27 



(0)a 



li 



(l)a 



7i 



(3)a 



(4)a 



,7 



(l)a 



•7< 



(2)a 



7t 



(4)a 



"7t 



27 



(6)a 



(6)a 



77 



(7)a 



+ 7t vy + 



7 



(6)a 



li 



(7)a 



'7i 



(9) a 



(W 


+ VaS u ) 


if )a ) ^ 


\ *ATtot 


^At a 


-7 t (8)a +7f )Q )y A H u 


iVtot 

N . 


^AH U ■ ( A - 75 ) 



We note at this point that Eq. (|A.73|) is in fact incomplete, since A£ a is also 
violated by EW sphalerons, however, no term accounting for these processes has 
been included. This can be corrected by writing instead BE for the flavour charge 
asymmetries Ya q = Yab/3 — Yal that are violated only by the RHN and RHSN 
interactions appearing in the r.h.s of Eq. (|A.75[) . Here Yal q = 2 \Y&.i a + Y A j J + 
Y&e a +Y/\e a and the contributions of e a and e a to the total flavour asymmetries have 
to be included because scatterings with the Higgs induced by the charged lepton 
Yukawa couplings can transfer part of the asymmetry generated in the SU{2)l 
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lepton doublets to the singlets. To take into account the EW sphalerons, we write 
down the complete BE for Yal„ and Yab- 

Y ALa = (Y ALa ) + (Y ALa ) , (A.76) 

V / pert \ / non— pert 



Yab = (Yab) , (A.77) 

V / non— pert 

where 'pert' refers to the violation of AL from perturbative interactions, i.e. the 
r.h.s. of Eq. (|A.73[) . while 'non-pert' refers to the violation of AL and AB from 
the non-perturbative EW sphaleron processes. Since the EW sphalerons conserve 
B/3 — L a we have 

\(Yab) -(Ya Lo ) =0, (A.78) 

O \ / non— pert \ / non— pert 

and then by taking the difference between Eqs. (|A.76|) and (|A.77[) with the proper 
factor of 1/3, we arrive at 



Y Aa = - [Yal J ■ (A.79) 

V / pert 

To get the BE for Y Aa in closed form, the density asymmetries Y A i a and Y A ^ 
multiplying the washout reactions on the r.h.s. of Eq. (|A.79[) must then be expressed 
in terms of Y Aa according to 

^Ai a = ^-Ai^Ae, ))ah u = Y^ C p u ^ a p ■ 

P P 

The values of the entries in the matrices A and C H ™ depend on the particular set 
of reactions that are in equilibrium when leptogenesis is taking place, and are given 
in|Appcndix B| 



Appendix A. 5. The non-superequilibration regime 

The BE describing the RHN and RHSN abundances in the non-superequilibration 
(NSE) regime (when ^ ^ fir) are still given by Eqs. (JA.39J) and (|A.41|) . The 
evolution of the flavour charges Y Aa is given by: 

Y Aa =-(E a + E^ , (A.80) 

where E a and E a are respectively given in Eqs. (|A.58j) - (|A.59[) . and the AL = 2 
W-terms as well as LFE effects (that are irrelevant at T J^ 10 7 GeV) have been 
neglected. To derive the BE for the evolution of Rb and R x defined in Eqs. (| 16911 - 
(|170[) we need to know by which amount these charges are violated in the different 
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scattering processes. This information is collected in the following table: 



Reaction 


Ai? B 


Ai? 3 


72 « 2 = 7 (iV ± ^ o Qu*) =i{N ± Q* o t a u*) =~f(N ± u <-> l* a Q) 





1 


7 f a = 7 (N ± ^u*I* a Q) 





1 


7t (0)Q = 7 (M a <-> Qu*\ = 7 (./Vl Q o Qu) 


-1 





lt 1)a = 7 f NQ <-> £*u*J = 7 [Nu <-> £* Qj 


-1 





7 t (2)Q = 7 (iVu ^> CQ) = 7 (#Q* ^> C«) 


-1 





7t (3)Q = 7 (iV£ Q ^> QU) 


-1 





7t (4)Q = 7 (TVw «• T a Q) = 7 (iVQ «■ T a u) 
7t (5)Q = 7 fAr±* a o Qu*) = 7 (iV±£ Q o Qu) 


-1 





1 


7t (6)Q ee 7 In ± u <+ 4Q) = 7 (iV±Q* o M 





1 


7t (7) " = 7 MV±Q ^> £«*) = 7 (#fcU <-> ^QJ 





1 


7f )Q ^ 7 f #fc% o Qu) 


2 


1 


7t (9)Q = 7 (/V ± Q o £ Q u) = 7 (jv±b o £ q q) 


2 


1 



The evolution equation for Yar b and Yar x are: 
Y ARb = ]T (2F a + F a ) - 7 f y A g, 






'AH a 



(A.81) 
(A.82) 



where the SE rates 7~ ff and 7° ff have been also included. F a and -F Q are given by: 



Fa = -f OJm. + Va*J - (t^"^ + 2 7t (4)Q ) 



and 



^ = 4 (*) 



7/ 



7w / y w t , 



(3)a 



. . Yn \ (4)q , 
1 + yeg ) 7* 



(A.83) 



9 I V e 9 



-2 - 



'JV 



(y A?a + Ha*. 



7^ 



(^A4, + ^Afl, 



X0)a VJV 



„(l)c 



„(2)c 



7* p? + 7* 1-7* +" o7t 



U 0)a 



1 (l)q yjv 

o 7* v eg 



li 



(8)a 



(8)a y AT tot 
ye« 
AT 

1^ 



^// 



(9)c 



7 Af Q 



2 y 



Wtot. ] -,(9)a 
eq I r* 



?Ai? u 



'AT 



1 / (0)a . (l)a . (2)a ^V < 

-o I 7* +7* +7* ^ ) l 2 ^Ai5„ - ^ 



V, 



iY 



(A.84) 
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For G a and G a we have: 



(5)q . (6)a . 1 (7)a Y Nu± 



Tt +Tt + oTt 




2 lz Y~ q 

N 




Vah u . (A.86) 

The density asymmetries of the five charges in the BE (JA.80I) . (|A.81|) and (|A.82|) 
define the basis Ya = {Ya„, YAfl B , Ya,r x } in terms of which one needs to express 
the five fermionic density-asymmetries Ya-j/v — {Yai , Ya§, Y^h }■ The relation is 
given by a 5 x 5 matrix defined according to: 

YAtp a = A a b YA b , 

and the numerical values of A a i, for different cases are given in Section 16.2.21 

With the inclusion of jf and jf 6 the BE fOO]) . (|X8l|) and (jX82|) are also 
valid in the SE regime. To verify this, one can compare the results obtained with 
the complete BE given above, assuming large in-equilibrium SE reactions 7~ ff and 
7® ff , with what is obtained from the set of BE specific for the SE regime (basically 
Eq. (|A.75[) with dY^a /dz — > —dY^a/dz). Of course, one also has to use the A ma- 
trix Eq. (|A.87|I and the corresponding A £ and C Hu matrices Eq. (|A.80j) appropriate 
for the specific temperature regime. For Case I (ft. e ,d Yukawa equilibrium) A is given 
in Eq. (fT92|) while A e and C" u are given in | Appendix B| in Eqs. (|B.8|) . For Case II 
(h e ,d Yukawa non-equilibrium) A is given in Eq. (|195j) and the corresponding SE 
matrices are given in Eqs. (|B.8[) and (|B.10[) . 
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Appendix B. Chemical equilibrium conditions in the SE regime 

In Section 16.2.11 items (l)-(5) a set of general constraints on particles/sparticles 
chemical potentials were given. At relatively low temperatures additional conditions 
hold, that are listed here. To simplify notations, in the following we denote chemical 
potentials with the same symbol that labels the corresponding fields: <f> = jjla,. 

6se-I- Equilibration of the particle-sparticle chemical potentials \i$ = fj,r (superequi- 

libration (SEJf^) is ensured when reactions like It -H- U are faster than the 
Universe expansion rate. These reactions are mediated by gaugino exchange 
but also require a chirality flip on the gaugino line, and thus are proportional 
to the soft mass m~. 

Fast reactions induced by the superpotential higgsino mixing term [iH u Hd im- 
ply that the sum of the up- and down-higgsino chemical potentials vanishes. 
Since /j, is expected to be of the same order than m§ it is reasonable to assume 
that both these reactions are in equilibrium in the same temperature range. The 
corresponding rates are given approximately by Tg ~ w?/T and T^ ~ fi 2 /T, 
and are faster than the Universe expansion rate up to temperatures 

r ~ 5 - io7 (5i^r &v - ^ 

When the temperature is sufficiently low that the limit rrig — > is not valid, 
then gauginos must be considered as Majorana fermions with an associated 
vanishing chemical potential: 

9 = 0. (B.2) 

Then, fast reactions £ 4-t £ + g, Q 4-> Q + g, H u ,d ++ H u .d + g etc. imply SE. 
Similarly, when the limit fi — > is not valid, fast H u <->• Hd reactions yield: 

H u + H d = 0. (B.3) 

6se-H. For temperatures satisfying Eq. (|B.1[) the MSSM has the same global anomalies 
than the SM: the EW SU(2) l -U(1)b+l mixed anomaly and the QCD chiral 
anomaly. EW and QCD sphalcron effects can be described by the effective 
operators Oew = H a (QQQ£ a ) and Oqcd = Tli{QQu c u d c Li ). Above the EWPT 
reactions induced by these operators are in thermal equilibrium and yield the 
conditions (compare with the corresponding NSE conditions (|185|) and (|186p ): 

CI. 

6Q-J2( u * + d *) = 0, (B.5) 

i 

where the same chemical potential has been assumed for the three quark dou- 
blets (Eq. (|184| l). which is always appropriate below the limit (|B.1|) . 
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Eqs. (|Wj) and (fTg2"j) - ([Tg4"|) . together with the SE conditions (|B~2J) - (|B~3)) . the two 
anomaly conditions (|B.4|) - (|B.5|) and the hypercharge neutrality condition f|180[) . give 
11 + 2 + 2 + 1 = 16 constraints for the 18 chemical potentials. Note however that 
there is one redundant constraint, since by summing up Eqs. (|182[) and (|183|) and 
taking into account conditions (J184I) . (JB.2I) . and (|B.3|) we obtain precisely the QCD 
sphalcron condition Eq. (|B.5|) . Therefore, like in the SM, we have three independent 
chemical potentials, which can be conveniently taken to be Y Aa = Yab/3 — Yal o : 



Fa 



\ Y, ( 2y AQ, + Y Aui + Y Ad J - {2Y Ata + Y Aea ) 



(B.6) 



The density asymmetries of the leptons and higgsino doublets, that weight the 
washout terms in the BE, can be expressed in terms of the dens ities of the anomaly 
free charges Eq. (|B.6[) by means of an A matrix!^ and C vector^ as: 

Y Aia = A l ap Y Af! , Y A ^ d = Cf- Y Aa . (B.7) 

In the following we give the results for A e and C Hud that refer to fermion states, 
and we recall that in the SE regime the density asymmetries of the corresponding 
scalar partners are given by Y As = 2 Y A f with the factor of 2 from statistics. 



Appendix B.l. Yukawa reactions in chemical equilibrium 

(7se-I) All Yukawa interactions in equilibrium. 

At temperatures below the limit in Eq. (|182[) all Yukawa interactions are in 
equilibrium and we have 

x /-221 16 16 

A ' = 97237 16 " 221 16 
^ Xl6< \ 16 16 -221 

C^=-C*« = ^(l, 1, 1). (B.8) 

In the SE regime Ya£ q = Y A t a + Y^j = 31a^„ and then the relation between 
Ya.£° and Y Aa is obtained by simply multiplying the A matrix in Eq. (|B.8[) by 
a factor of 3. This gives the same A matrix obtained in the non-supersymmetric 
case in the same regime (see e.g. Eq. (4.13) in Ref . IM|) . The C matrix (multiplied 
by the same factor of 3) differs from the C matrix of the non-supersymmetric 
result (that is of course given for the scalar density- asymmetry Y A h) by a factor 
1/2. This is simply because Y A ^ U = Y A ^ + Y A h u — ^Y A h u . These results 
agree with Rcf.l89[ and hold in general for supcrsymmctry within the SE regime. 
(7se-II) Electron and up-quark Yukawa reactions out of equilibrium. 

For temperatures above 10 5 (1 + tan 2 /3) GeV the interactions mediated by the 
electron Yukawa h e drop out of equilibrium, and one of the conditions Eq. (|18ip 
is lost. However, in the effective theory one can then set h e — » and one 
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global symmetry is gained. This corresponds to chiral symmetry for the R- 
handed electron, that here translates into a symmetry under phase rotations of 
the e chiral multiplct. Conservation of the corresponding charge ensures that 
An e + Ang = 3An, e is constant, and since leptogenesis aims to explain dy- 
namically the generation of a lepton asymmetry we set this constant to zero, 
implying a vanishing chemical potential for the R-handcd electron e = 0. In 
this way the chemical equilibrium condition that is lost is replaced by a new 
condition corresponding to the conservation of a global charge, and the three 
non-anomalous charges (|B.6[) are again sufficient to describe all the density 
asymmetries. At temperatures above T ~ 2 ■ 10 6 GeV interactions mediated by 
the up-quark Yukawa coupling h u drop out of equilibrium. In this case how- 
ever, by setting h u — > no new symmetry is obtained, since chiral symmetry 
for the R-handcd quarks is anomalous, and the corresponding charge is violated 
by QCD sphalerons. However, after dropping the first condition in Eq. (|182|) 
for m = u, the QCD sphaleron condition Eq. (|B.5j) ceases to be a redundant 
constraint, with the result that also in this case no new chemical potentials are 
needed to determine all the particle density asymmetries. In this case, the A 
and C are given by 

/ 1221 156 156' 

A e = 111 -910 52 

3X2886 V HI 52 -910, 

C H U = _ c H d = __ / 3? 52 52 ) /g g x 

2886 *-'''' y ' 

(7seTII) First generation Yukawa reactions out of equilibrium. 

At temperatures T > 4 • 10 6 (1 + tan 2 j3) GeV, also the d-quark Yukawa coupling 
can be set to zero (to remain within the SE regime we assume tan/3 ~ 1). In this 
case the equilibrium dynamics is symmetric under the exchange u -f-> d (both 
chemical potentials enter only the QCD sphaleron condition Eq. (|B.5[) with 
equal weights) and so must be any physical solution of the set of constraints. 
Thus, the first condition in Eq. (|183[) can be replaced by the condition d = u, 
and again three independent quantities suffice to determine all the particle 
density asymmetries. The A and C matrices in this case are: 

120 120' 
-688 28 
28 -688, 

C H - -C H * = ^ (37, 52, 52) , (B.10) 

that agrees with what is obtained in the non-supcrsymmctric case (see Eq. (4.12) 
of Rcf. l64| after the factor of 1/2 relative to the higgsinos discussed above in 
(7se-I) is taken into account. 
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